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1  INTRODUCTION 


Manual  flight  control  system  design  for  fighter  aircraft  continues  to  be  one  of  the 
most  demanding  problems  in  the  world  of  automatic  control.  It  is  the  job  of  the 
manual  flight  control  system  to  provide  a  satisfactory  dynamic  response  to  pilot 
inputs.  The  problem  is  inherently  multivariable,  that  is,  a  controller  must  drive 
multiple  effectors  based  on  information  from  multiple  sensor  and  command  inputs. 
Modem  fighter  aircraft  dynamics  generally  have  highly  coupled  dynamics  that  are 
both  uncertain  and  nonlinear. 

In  the  case  of  such  multivariable  problems,  the  only  efficient  means  of  obtaining  a 
solution  is  to  use  multivariable  control  design  techniques.  Classical  control 
techniques  that  close  a  .single  feedback,  loop  at  a  time  are  inefficient  for  designing 
control  systems  for  plants  with  multiple  inputs  and  outputs  and  strong  coupling 
between  the  loops.  Further,  classical  techniques  do  not  address  robustness  in 
multivariable  systems  adequately,  as  they  only  address  stabilit’  margins  in  single 
loops  and  neglect  simultaneous  perturbations  in  .several  loops.  Multivariable  control 
techniques  can  provide  an  efficient  means  of  finding  control  laws  for  complex  systems 
and  make  it  straightforv/ard  to  incorporate  specific  performance  and  robustness 
requirements. 

The  magnitude  of  the  manual  flight  control  problem  is  diiven  by  the  nonlinear  and 
uncertain  nature  of  aircraft  dynamics.  Linear  models  of  these  system.s  are  only  valid 
for  small  regions  about  trim  conditions.  The  conventional  solution  to  this  problem  is 
to  perform  point  design.s  for  a  large  set  of  trim  conditions  and  then  construct  a  gain 
schedule  by  interpolating  gains  with  respect  to  flight  condition.  This  ad  hoc  procedure 
is  time  consuming  and  expensive,  but  is  well  accepted  and  has  yielded  satisfactory 
results  for  dozens  of  aircraft. 

One  advantage  of  classical  controllers  is  that  they  are  straightforward  to  schedule 
with  flight  condition  becau.se  there  is  a  clear  meaning  for  each  one  of  the  gains  or 
dynamic  filter  elements.  The  use  of  muitivari;ible  control  in  full  envelope  flight  control 
law  design  is  difficult  because  the  resultant  controllers  are  usually  dynamic 
compensators  in  state  space  form  and  ofte  i  are  of  high  order.  Huch  control  laws  are 
extremely  difficult  to  schedule  with  varying  flight  condition,  as  the  elements  in  the 


state  space  description  of  the  controller  do  not  vary  smoothly  with  flight  condition, 
and  their  relation  to  any  particular  parameter  is  not  clear. 

It  is  tempting  to  represent  the  variation  in  aircraft  dynamics  with  flight 
condition  as  a  model  uncertainty  and  design  a  single,  fixed  control  law  for  a  nominal 
aircraft  model.  This  would  eliminate  the  need  to  schedule  the  controller  with  flight 
condition.  However,  making  the  control  law  robust  to  all  changes  in  the  dynamics  Hue 
to  changes  in  flight  condition  would  mean  sacrificing  performance.  No  single,  fixed 
control  law  will  give  maximum  performance  across  a  wide  flight  envelope.  In  addition, 
although  the  dynamics  at  a  particular  flight  condition  are  usually  uncertain  to  some 
degree,  there  is  usually  a  reasonable  amount  of  information  about  how  the  flight 
dynamics  change  as  the  operating  condition  of  the  vehicle  changes.  Some  alternative 
to  representing  the  plant  parameter  variation  as  model  uncertainty  is  required. 

An  obvious  solution  to  this  problem  is  to  design  constant  gain  multivariable 
controllers  whose  elements  can  be  scheduled  with  flight  condition.  This  is  a  useful 
alternative  in  some  cases,  although  the  number  of  techniques  available  for  this  are 
somewhat  limited,  and  in  some  cases  these  have  restrictive  assumptions  on  the 
design  model  such  as  requiring  state  feedback.  Furthermore,  many  of  the  most 
powerful  robust  control  design  approaches,  which  account  for  modeling  uncertainty 
using  frequency  dependent  bounds  and  guarantee  robustness  via  the  small  gain 
theorem,  produce  dynamic  controllers.  While  it  is  becoming  more  straightforward  to 
obtain  multivariable  controllers  for  linear  time  invariant  plants  that  meet  specific 
performance  and  robustness  requirements,  using  multivariable  control  as  a  tool  in 
the  overall  design  of  a  nonlinear  system  by  individual  point  designs  is  still  an  issue. 


Recently  more  attention  has  been  given  to  applying  a  more  mathematical  basis  to 
the  issue  of  gain  scheduling  flight  control  systems.  In  [1.1],  the  author  transforms 
plant  dynamics  into  a  quasi-linear  parameter  varying  form  in  which  dynamics  depend 
on  exogenous  variables  that  are  unknown  but  can  be  measured.  This  transformation 
allows  explicit  gain  scheduling  relationships  to  be  derived  which  preserve  stability  in 
the  presence  of  time  varying  plant  parameters.  In  [1.2],  the  author  describes  a 
systematic  approach  for  automating  gain  schedule  calculations.  The  approach 
guarantees  both  stability  and  performance  using  structured  .singular'  value  theory. 
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There  are  a  few  examples  ot’ attempts  to  gain  schedule  multivariable  controllers  in 
the  literature.  In  [1.3],  the  authors  designed  control  laws  at  four  widely  spaced 
operating  points  for  a  pitch  axis  autopilot  for  a  highly  maneuverable  missile.  The 
plant  dynamics  are  defined  as  an  explicit  function  of  some  parameters  such  as  angle 
of  attack  and  dynamic  pressure.  The  linear  controllers  at  the  different  equilibrium 
points  are  interpolated  dependi.ng  upon  the  operating  condition.  Unfortunately, 
describing  the  aircraft  dynamics  as  an  explicit  function  of  all  its  parameters  is  not  a 
straightforward  approach.  In  [1.41,  the  authors  design  controllers  for  a  V/STOL 
aircraft  at  discrete  points  and  then  switch  between  them  based  upon  operating 
condition.  Finally,  in  [1.5]  and  [1.6],  the  authors  design  flight  control  laws  for  a  test 
aircraft  using  integral  LQ  regulators.  The  controllers  found  were  constant  gains  plus 
integrators,  and  were  scheduled  as  a  function  of  dynamic  pressure. 

Ail  alternative  to  gain  scheduling  is  to  use  control  design  methods  which  directly 
consider  the  nonlinear  nature  of  the  problem.  Adaptive  control  has  received  much 
attention  recently  and  shows  promise.  Nonlinear  dynamic  inversion  has  been  applied 
successfully  to  a  number  of  flight  control  problems  [1.7,  1.8,  1.9], 

In  this  report,  the  problem  of  multivaricihle  gain  scheduling  is  addressed  by 
using  an  irmcr/outer  loop  control  law  structure  that  separates  the  issues  of  scheduling 
the  controller  and  designing  control  laws  to  meet  specific  perfonnanco  and  robustness 
requirements.  The  inner  loop  consists  of  a  control  selector  that  transforms 
generalized  control  commands  into  actual  control  effector  commands  and  feedback 
compensation  that  minimizes  the  relative  error  between  the  inner  loop  dynamics  at 
different  flight  conditions.  Both  the  control  selector  and  the  inner  loop  feedback  are 
functions  of  flight  condition.  The  idea  is  to  use  the  inner  loop  to  make  changes  in  the 
aircraft  dynamics  transparent  to  the  outer  loop.  The  outer  loop  consists  of  a  fixed 
dynamic  compensator  to  meet  the  performance  and  robustness  requirements  of  the 
system. 

Two  design  examples  are  presented  to  illustrate  the  effectiveness  of  the  inner/outer 
loop  approach.  Manual  flight  cor  -rol  systems  are  designed  for  the  lateral/directional 
axes  of  the  VISTA  F-16  test  vehicle  and  for  the  longitudinal  and  lateral/directional 
axes  of  a  supermaneuverable  lighter  aircraft.  These  two  vehicles  bring  out  a  number 
of  interesting  flight  control  problems  including  effector  limiting,  thrust  vectoring,  and 
redundant  controls.  Three  different  inner  loop  de.sign  approaches  are  used  to  explore 


equalization  methods.  Dynamic  inversion  is  used  for  the  F-16  inner  loop  design,  a 
reduced  order  method  is  used  for  the  longitudinal  supermaneuverable  vehicle  inner 
loop  design,  and  eigenstructure  assignment  is  used  for  the  iateral/directional 
supennaneuverable  vehicle  inner  loop  design.  Flying  qualities  are  built  into  each 
design  using  an  implicit  model  following  p-synthesis  outer  loop  formulation.  An  ideal 
model  of  the  aircraft  response  which  represents  tfie  desired  flying  qualities  is  included 
in  the  S3mthesis  model.  The  infinity  norm  of  the  transfer  function  from  the  command 
input  and  the  frequency  weighted  error  between  the  ideal  model  response  and  the 
actual  model  response  is  minimized.  Performance  and  robustness  analysis  is 
performed  for  wade  ranges  of  operating  conditions  for  each  of  the  designs.  Nonlinear 
simulation  results  are  shown  to  demonstrate  that  the  inner/outer  loop  design 
approach  yields  controllers  that  perform  well  in  a  highly  dynamic  and  nonlinear 
environment. 

Chapter  2  of  this  document  presents  the  basic  technical  background  behind  the 
manual  flight  control  problem  and  the  theory  used  in  the  applications.  Chapter  3 
introduces  the  inner/outer  loop  robust  control  design  methodology.  Chapter  4  is  a 
detailed  application  of  this  methodology  to  the  design  of  a  lateral/directional  manual 
flight  control  system  for  the  VISTA  F-16  aircraft.  Chapter  5  is  a  detailed  application 
of  this  methodology  to  a  supermaneuverable  F-18  manual  flight  control  problem. 
Chapter  6  is  a  conclusion  which  sums  up  the  results  and  recommends  future 
research  directions. 
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2  TECHNICAL  PRELIMINARIES 


2.1  AIRCRAFT  DYNAMICS 

The  nonlinear  equations  of  motion  for  an  aircraft  are  derived  using  Newton's 
Second  Law  of  motion.  That  is,  the  total  sum  of  all  external  forces  acting  on  a  body 
must  equal  the  time  rate  of  change  of  linear  momentum  and  the  total  sum  of  all 
external  moments  acting  on  a  body  must  equal  the  time  rate  of  change  of  angular 
momentum.  Some  critical  assumptions  necessary  to  simplify  the  derivation  of  the 
aircraft  equations  of  motion  arc:  the  aircraft  is  a  rigid  body,  tlie  mass  of  the  aircraft 
remains  constant  with  time,  and  the  earth  provides  a  fixed  inertial  reference  frame 
[2.1].  Fig.  2.1  shows  the  body  axis  forces,  moments,  angles,  velocities,  and  rotational 
rates  necessary  to  describe  the  motion  of  an  aircraft. 


The  nonlinear  equations  of  motion  are  made  up  of  three  translational  and  three 
rotational  equations.  The  translational  equations  of  motion  are: 

m[  U  +  QW  -  RV  4-  gsinG]  =  X 
m[V  +  KU  -  PW  -  geos  0  sin  fid  =  Y 
m[W  +  PV  -  QU  -  geos  0  CO.S  fl>1  =  Z 


(i 


(2.1) 


The  rotational  equations  of  motion  are: 


Plx  -  RIxz  +  QRdz  -  ly)  -  PQIxz  =  L 

Qly  +  PKdx  -  Iz)  -  =  M  (2.2) 

RIz  -  PIxz  +  PQCIy  -  Ix)  +  QRIxz  =  N 

U,  V,  and  W  are  the  translational  velocities,  P,  Q,  and  R  ai-e  the  rotational  rates,  m  is 
the  aircraft  mass,  Ix,  ly,  Iz,  ^n^d  Ixz  are  the  moments  of  inertia,  g  is  gravity,  and  X,  Y, 
Z,  L,  M,  and  N  are  the  external  forces  and  moments  due  to  the  aerodynamics  and 
propulsion.  Eqs.  (2.1)  and  (2.2)  completely  describe  the  motion  of  an  aircraft.  The 
Euler  angles,  0,  <!),  and  T,  describe  the  orientation  of  the  aircraft  with  respect  to  the 
Earth. 


* 

(h  =  P  +  Q  tan  0  sin  <h  +  R  tan  0  cos 

0  =  Q  cos  <h  -  R  sin  «!>  (2,3) 

.  R  cos  (h  Q  sin  (h 
cos  0  cos  0 

O  is  the  roll  angle,  0  is  the  pitch  angle,  and  T  is  the  yaw  angle  [2.2J, 

2.1.1  TRIMMED  EQUATIONS 

A  trimmed  condition  is  a  local  equilibrium  condition  at  which  all  of  the  linear  and 
rotational  accelerations  in  eqs.  (2.1)  and  (2.2)  aie  zero.  Let  the  trimmed  flight 
condition  be  described  as  variables  with  zero  subscripts.  The  general  nonlinear  trim 
equations  of  motion  are  [2.2] 


m[QoW()  -  RoVo  +  g.sin©,)!  = 

m[RDUo  -  P()Wo  -  geos  00  .sin 

m[PoVo  -  QqUo  '  geos  0()  cos  (l)o]  =  Zq 
Qoltodz  -  ly)  -  PqQo^xz  =  Lq 
P()Podx  -  Iz)  -  l^o^Ixz  +  Pfl'^lxz  =  Mq 
P()Qody  ■  Ix)  +  Qol^lxz  =  Nq 
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While  these  equations  allow  for  solutions  with  non-zero  rotational  rates  and  side 
velocity,  additional  assumptions  are  usually  applied  to  further  simplify  the  trim 
solution.  The  most  common  of.  these  is  straight  and  level,  coordinated  flight.  In  this 
ca.se  eq.  (2.4)  simplifies  to 


mg  sin  ©0  =  Xo 

-  mg  cos  ©0  sin  <I)o  =  Yq 

-  mg  cos  ©0  cos  cl>o  ^  Z^)  (2.5) 

0  =  Lo 
0  =  Mo 
0  =  No 

That  is,  the  available  aerodynamic  control  surfaces,  throttle  controls,  etc.  must  be 
set  to  cancel  the  gravitational  and  aerodynamic  forces  and  moments  on  the  aircraft. 

2.1.2  LONGITUDINAI.  LINEAR  EQUATIONS  OF  MOTION 

The  d3Tiamics  of  a  rigid  aircraft  are  described  by  the  six  simultaneous  nonlinear 
equations  as  shown  in  eqs.  (2.1)  and  (2.2).  These  equations  can  be  programmed  and 
digitally  integrated  by  a  computer  to  simulate  aircraft  motion.  Because  most 
analysis  and  design  tools  require  a  linear  representation  of  a  system,  it  is  useful  to 
make  some  assumptions  to  linearize  the  aircraft  equations  of  motion.  The  first  step 
is  to  break  the  six  equations  into  two  sets  of  simultaneous  equations,  three 
longitudinal  and  three  lateral/directional  I2.1J. 

To  develop  the  longitudinal  linear  equations  of  motion,  it  is  assumed  that  the 
aircraft  is  in  straight  and  level  unaccelerated  flight.  The  only  disturbances  on  the 
system  considered  are  external  forces  X  and  Z,  and  external  moments  M.  These 
disturbances  on  the  equations  of  motion  do  not  create  any  sideforce,  Y,  or  any  rolling 
moment,  L,  or  yawing  moment,  N.  Itoll  rate,  yaw  rate,  and  aide  velocity  remain 
undisturbed  so  three  of  the  equations  can  be  neglected.  'I'he  remaining  equations  are 
simplified  because  V  -  P  =  R  =  =  0. 


m[U  QW  +  gsinGJ  =  X 

m[W-  QU  -  geos©]  =  Z  <2.6) 

QIy=  M 


H 


Fig.  2.2  shows  the  orientation  of  the  longitudinal  variables  with  respect  to  the  aircraft 

W 

and  its  total  velocity  vector,  V^.  The  total  angle  of  attack,  wi,  is  equal  to  sin'^C^-), 


Fig.  2.2  Longitudinal  Axis  System 

By  assuming  that  the  aircraft  is  in  equilibrium,  the  total  linear  and  angular  velocities, 
Fiuler  angles,  and  the  total  external  forces  and  moments  can  be  written  as  the  sum  of 
their  equilibrium  values  and  pcrturbational  values. 

U  =  LJo  +  u  W  =  W()  +  w 

Q  -  Qo  M  =  M()  +  dM  (2.7) 

X  =  Xy  +  dX  0  0Q  +  H 


If  only  small  disturbances  about  the  equilibrium  point  arc  considered,  the  product  of 
the  perLui’batioaal  values  cun  be  assumed  to  be  very  small  and  neglected.  Tlie 
angular  values  between  the  equilibrium  and  disturbed  conditions  are  also  assumed  to 
be  small.  By  writing  the  equations  in  the  stabihy  axis,  Wy  can  be  set  to  zero.  The 
longitudinal  force  and  moment  eciuations  in  (2.6)  can  now  be  represented  as 


mfu  +  (gc()s0())Ol  n:  dX 

m[w  -  qU()  +  (gsin0o)0J  -  dZ  (^*8) 

qly=  dM 
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The  equations  in  (2.8)  are  linear  with  respect  to  the  perturhalional  variables. 
Continued  development  of  the  linear  longitudinal  ecjuations  of  motion  involves 
expressing  expanded  representations  of  the  external  forces  and  moments  in  tenns  of 
the  changes  in  them  resulting  from  the  perlurhations  in  the  linear  and  angular 
velocities.  In  other  words  the  partial  derivatives  of  (he  forces  and  moments  are  taken 
with  respect  to  the  perturbational  variables.  A  detailed  development,  of  those 
equations  can  be  found  in  [2.1]  ,  12. 2|,  or  [2.8J.  'fhe  resulting  longitudinal  linear 
equations  of  motion  can  be  expressc'd  in  state  space  form: 
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where  a  =  -  exo-  X^,  'Aa,  Mq,  etc.  are  longitudinal  stability  derivatives  tind  Xsi,  Xs». 

Mfini)  etc.  are  longitudinal  control  derivatives.  The  variables  8),  ,  ^in  ^*^0 

changes  in  external  control  inputs  such  as  elevator,  Haps,  and  thrust  vectoring, 

Typically,  an  open  loop  aircraft  with  a  classical  configuration  operating  in  a 
trimmed  condition  at  a  conventional  flight  condition  will  exhibit  two  longitudinal  modes 
of  motion:  the  short  period  and  phugoid.  The  short,  period  mode  is  normally  fa.st  and 
oscillatory  and  takes  place  at  nearly  constant  speed.  It  is  dominated  by  angle  of 
attack  and  pitch  rate  response.  The  stability  derivative  drives  the  natural 
frequency  of  this  mode.  The  phugoid  mo<k*  is  normally  slow,  o.scillatory,  and  lightly 
dEunped  and  takes  place  at  nearly  constant  angle  of  attack  |2..‘}j  b'or  a  manual  flight 
control  system  design  problem,  the  short  period  is  the  primary  mode  of  interest 
because  it  dominates  the  aircraft's  response  to  pilot  inputs.  The  phugoid  mode  is  the 
most  important  in  autopilot  designs. 

Because  the  short  period  mode  is  virtually  decoupled  from  the  speed  and  pitcli 
angle  response  of  the  aircraft,  it  is  possible  to  reduce  the  model  in  eq,  (2.9)  to  a  second 
order  short  period  approximation. 
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This  model  provides  an  accurate  measure  of  the  aircraft’s  transient  response  to  small 
amplitude  inputs  over  a  short  time  frame  (~  5  seconds).  Because  the  primary  focus 
of  this  document  is  manual  flight  control  system  design,  the  phugoid  approximation  is 
not  described  in  detail  here.  The  reader  is  referred  to  [2.1],  [2.2]  or  [2.3]. 

2.1.3  LATERAL/DIRECTIONAL  LINEAR  EQUATIONS  OF  MOTION 

To  develop  the  lateral/directional  equations  of  motion,  it  is  again  assumed  that  the 
aircraft  is  in  straight  and  level  unaccelerated  flight.  Only  disturbances  on  the 
external  force  Y  and  external  moments  L  and  N  are  considered  [2.1].  The  three 
equations  used  to  derive  the  decoupled  lateral  equations  of  motion  are: 

m[V  +  RU  -  PW  -  geos  0  sin  <I>]  =  Y 

PIx  -  Rlxz  +  QRdz  -  ly)  -  PQlxz  =  L  (2.11) 

RIz  -  PIxz  +  PQdy  -  lx)  +  QRIxz  =  N 


Fig.  2.3  shows  the  orientation  of  the  lateral/directional  variables  with  respect  to  the 

V 

aircraft  and  its  total  velocity  vector,  V^j'.  The  angle  of  sideslip,  p,  is  equal  to  sin'l(^). 


Fig.  2.3  Lateral  Axis  vSystem 
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The  aircraft  is  assumed  to  be  in  a  straight  and  level  equilibrium  condition,  so  the  total 
linear  and  angular  velocities,  Euler  angles,  and  the  total  external  forces  and  moments 
are  represented  as  the  sum  of  their  equilibrium  value  and  perturbational  values. 


P  =  Pq  +  p 
V  =  Vq  +  V 
L  —  Ijq  *4"  dh 


R  -  Rq  +  r 
y  =  Yo+dY 
N  =  No  +  dN 
4^  =  M>0  +-  Vj/ 


(2.12) 


Zero  total  pitch  rate,  Q,  and  zero  equilibrium  side  velocity,  Vq,  can  be  assumed  as  a 
consequence  of  a  decoupled  trim  condition.  The  unaccelerated  flight  condition  dictates 
that  equilibrium  roll  and  yaw  rates  must  be  zero,  Pq  =  Rq  =  0.  Small  perturbations 
about  the  trim  condition  are  assumed,  so  the  product  of  these  perturbations  is 
neglected  [2.1],  In  the  design  of  manual  flight  control  systems  for  fighter  aircraft, 
angles  of  attack  beyond  what  may  be  considered  small  mu.st  often  be  addressed.  In 
order  to  more  clearly  express  an  angle  of  attack  dependency,  the  translational 
equation  can  be  rewritten  in  tem.s  of  the  stability  axis  variables  sideslip  and  angle  of 
attack.  To  maintain  decoupling  in  the  equations  of  motion,  angle  of  attack  is 
assumed  constant,  =  ao-  lateral  force  and  moment  equations  in  (2.12)  can 
now  be  expressed  as 


mUofP  +  r  cos  exq  -  p  sin  oq  -  g(!)/Uo  cos  Oq]  =  dY 

pix  ■  rlxz  =  dL  (2.13) 

•  « 

rlz  ■  plxz  =  dN 


The  lateral  equations  in  (2.13)  are  linear  with  respect  to  the  perturbational  variables. 
They  are  shown  here  in  a  body  axis  system  to  show  the  (Xq  dependence  in  the  (3 
equation.  Changes  in  expanded  representations  of  the  external  forces  and  moments 
resulting  from  perturbations  in  the  linear  anu  angular  lateral  velocities  must  be 
derived.  Detailed  derivations  are  given  in  (2.1]  ,  [2.2],  and  [2,3].  The  resulting 
lateral/directional  linear  equations  of  motion  can  be  expressed  in  state  space  form: 
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Yp,  Lp,  N,-,  etc.  are  the  lateral  stability  derivatives  and  Ygi,  L52,  Msm,  etc.  are  the 
lateral  control  derivatives.  Note  that  Yp  and  Y,.  are  assumed  to  be  neglig^ible. 

The  lateral/directional  equations  of  motion  are  usually  characterized  by  three 
distinct  modes  of  motion;  a  first  order  spiral  mode,  a  first  order  roll  mode,  and  a  second 
order  Dutch  roll  mode.  The  spiral  mode  is  dominated  by  bank  and  heading  angle,  <])  and 
v;/,  while  sideslip,  p,  is  small.  It  noi*mally  has  a  long  lime  constant  and  may  be  slowly 
divergent.  In  manual  flight  control  system  designs,  the  spiral  mode  is  often  neglected. 
Bank  and  heading  angle  are  trajectory  states  which  should  naturally  evolve  during 
manuevering  flight  and  therefore  should  not  be  regulated  by  the  controller.  These 
states  are  important  in  autopilot  designs.  The  roll  mode  is  usually  fast  and  stable.  It 
is  dominated  by  rolling  motion  and  the  stability  derivative  Lp,  The  Dutch  roll  mode  is 
dominated  by  sideslip  and  yawing  motion  if  the  derivative  L|^,  corresponding  to 
dihedral  effect,  is  small.  It  is  usually  oscillatory  and  lightly  damped.  The  term  Np 
drives  the  frequency  of  this  mode,  wliile  Yp  and  N,.  drive  the  damping  [2.3J. 


2.2  FLYING  QUALITIES 


Fl3ring  qualities  include  everything  that  is  involved  in  the  safe  flight  and  effective 
performance  of  an  aircraft,  from  the  pilot’s  point  of  view.  Military  standards  for  fljdng 
qualities  provide  guidelines  for  analytical  parameters  which  have  been  correlated  to 
safety  and  mission  performance.  MIL-STD-1797A,  Flying  Qualities  of  Piloted 
Vehicles  [2.4],  provides  guidance  for  U.S.  military,  fixed  wing  aircraft. 

In  order  to  interpret  requirements,  it  is  necessary  to  classify  an  aircraft's  mission. 
Requirements  for  a  light  trainer  aircraft  obviously  should  not  bo  the  same  as  those 
for  a  heavy  strategic  bomber.  An  airplane  can  bo  put  into  one  of  the  following 
categories: 

Class  I  Small  Light  Aircraft 

Class  II  Medium  Weight,  Low-to-Medium  Maneuverability  Aircral’t 

Class  III  Large,  Heavy,  Low-to-Medium  Maneuverability  Aircraft 

Class  rV  High-Maneuverability  Aircraft 

Another  important  discriminator  in  flying  qualities  requirements  is  flight  phase. 
Categories  of  flight  phases  separate  flight  conditions  and  tasks  into  groups  which 
require  similar  aircraft  responses  and  pilot  workloads.  These  categories  are  further 
divided  into  nontenminal  and  terminal; 

Nonterminal 

Category  A:  Flight  phases  requiring  rapid  maneuvering,  precision 
tracking,  or  preense  flight  path  control. 

Category  B;  Flight  phases  that  are  normally  accomplished  using 

gradual  maneuvers  without  proension  tracking.  Accurate  flight 
path  control  may  be  required. 

Terminal 

Category  C;  Flight  phases  that  are  normally  accomplished  using 

gradual  maneuvers  and  requiring  accuram  flight  path  control. 


The  qualitative  degree  of  acceptability  for  hajidiing  qualities  is  given  in  tenns  of  levels. 
These  levels  specify  the  adequacy  of  the  aircraft  response  in  meeting  mission 
requirements  and  ctre  based  on  the  Cooper-Harper  pilot  opinion  rating  scale  [2.4]. 

Level  1  ■  Satisfactory:  Flying  qualities  clearly  adequate  for  the  mission  flight 
phase.  Desired  performance  achievable  with  no  more  than  minimal  pilot 
compensation. 

Level  2  -  Acceptable:  Flying  qualities  adequate  to  accomplish  the  mission 
flight  phase  with  some  increase  in  pilot  workload  or  degradation  in 
mission  effectiveness. 


Level  3  -  'Controllable:  Fl34ng  qualities  such  that  the  aircraft  can  be 

controlled  in  the  context  of  the  mission  flight  phase,  but  with  excessive 
pilot  workload  or  inadequate  mission  effectiveness. 


Most  flying  quality  specifleations  have  been  correlated  to  modal  parameters. 
These  modal  parameters  come  from  specific  low  order  forms  for  the  transfer 
functions  that  describe  the  aircraft  response  to  pilot  inputs.  In  general,  the  transfer 
functions  that  describe  the  dynamics  of  an  aircraft  include  not  only  the  open  loop 
aerodynamics,  but  also  actuator  dynamics,  controller  dynamics,  structural  modes, 
sensor  dynamics,  and  so  on.  These  high  order  transfer  functions  must  be  translated 
into  equivalent  low  order  forms  so  that  the  modal  parameters  which  define  the  flying 
qualities  may  be  identified.  These  low  order  forms  are  called  low-order-equivalent- 
systems  (LOFS).  The  LOKS  is  derived  by  matching  the  high-order-system  (HOS) 
over  a  specific  frequency  range,  usually  0.1  lad/.s  <  co  <  10  rad/s.  The  fit  is  measured 
by  a  weighted  sum  of  squares  function  of  the  differences  in  magnitude  and  phase 
between  the  LOES  and  HOS  at  n  discrete  frequencies.  This  function  is  given  by 
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(2.15) 


The  relative  weighting  between  the  gain  and  phase  errors  dictate.s  that  1  db  of  gain 
mismatch  is  roughly  equal  to  7  deg  of  phase  mismatch.  Any  value  of  the  mismatch 
function  less  than  10  is  acceptable.  Because  values  greater  than  10  may  or  may  not 
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be  acceptable,  an  additional  measure  of  the  IjOES  fit  may  be  used.  Fig.  2.4  shows  an 
approximate  bound  on  the  maximum  unnoticeable  dynamics  between  the  LOES  and 
HOS.  These  bounds  represent  the  level  of  mismatch  that  will  result  in  a  1  pilot  rating 
change  on  the  Cooper-Harper  scale.  A  pilot  is  most  sensitive  to  dynamics  in  the 
region  from  1  to  4  rad/s  [2.4]. 


Frequency  (rad/s)  Frequency  (rad/s) 


Fig.  2.4  Bounds  on  Maximum  Unnoticeable  Added  Dynamics 

It  should  be  noted  that  the  military  standard  allows  for  a  good  amount  of  tailoring. 
Due  to  the  bi’oad  possibilities  for  interpretation,  only  one  possible  set  of  requirements 
for  Class  IV  aircraft  in  Categoiy  A  flight  phases  is  presented  here. 

2.2.1  LONGITUDINAL  RESPONSE  REQUIREMENTS 


The  fl5dng  qualities  requirements  for  the  pitch  axis  are  determined  by  the  response 
of  the  aircraft  to  a  longitudinal  input  from  the  pilot,  usually  a  stick  force  or  deflection. 
The  low  order  modal  form.s  required  for  the  LOES  fits  arc 


_ ^0^ 

Sp(s)  (s2  4-  2^pa)pS  +  0)2)  (s2  +  2^,;pO)spS  +  o)  ) 


(2.1B) 


n^ls)  (s  +  lAl'np 

5p(s)  (s2  +  2^p0)pS  +  0)2)  (s2  +  2^j^pO)j,pS  +  0)2,1 


(2.17) 
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q  is  pitch  rate,  is  the  normal  acceleration  at  the  instantaneous  pitch  center  of 
rotation,  and  8p  is  a  pilot  pitch  stick  input. 

The  primary  flying  quality  measure  for  the  short  period  mode  is  the  control 
anticipation  parameter  (CAP).  A  value  for  CAP  may  be  estimated  by  the  expression 

Oi? 

CAP  =  ^-  where  (n/a)  -  (V-yg)(lAro.d  (2.18) 

(n/cx) 

is  the  airspeed  along  the  tliglit  path  (ft/s)  and  g  is  gravitational  acceleration  (ftys^). 
The  Level  1  cmd  2  requirements  for  the  estimated  CAP  parameters  are  shown  in  Fig. 
2.5. 


Fig.  2.5  Short  Period  Frequency  Requirements 

The  requirement  specifies  that  oigp  >1.0  rad/s  for  Level  1  and  cOgp  >  0.6  rad/s  for  Level 
2.  Short  period  damping,  ^^p,  must  be  between  0.35  and  1.3  for  Level  I  and  0.25  and 
2.0  for  Level  2.  Another  measure  of  short  period  flying  qualities  is  the  product  of  oigp 
and  T02.  Specifications  for  this  product  can  be  used  in  place  of  or  as  a  supplement  to 
CAP  requirements.  Fig.  2.6  shows  the  values  of  oi^pTy^  necessary  for  Level  1  and  2 
flying  qualities. 
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Fig.  2.6  Short  Period  Damping  Requirements 


The  LOES  parameter  is  the  equivalent  pitch  time  delay.  This  value  repi’osents 
lags  and  time  delays  in  the  control  system  and  the  phase  loss  due  to  the  liigh  order 
dynamics  that  are  not  represented  in  the  low  order  model.  Excessive  time  delay  can 
lead  to  pilot  induced  oscillation  problems.  Table  2.1  shows  the  requirements  for 
maximum  pitch  time  delay  . 


Table  2. 1  Pitch  Time  Delay  Requirements 


Level 

Maximum  'fime  Delay 

1 

0.10  sec 

2 

0.20  sec 

3 

0,25  sec 

TVio  lontT 

^  AA\y  a\Jaa^  i/V'AAAA  &  V  • « K.r  ^  V'C^tk  ijl  V/X  C'AAV' 

parameters  in  the  LOES.  The  stability  requirements  in  Table  2.2  are  recommended 
for  any  longitudinal  mode  with  a  period  greater  than  15  seconds. 
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Table  2.2  Phugoid  Stability  Roquiremenl.s 


Level 

1 

>  0.04 

2 

Cn  >  0.0 

6 

T,i  >  55  sec 

Tj  is  the  time  to  double  for  an  unstable  oscillation.  No  aperiodic  longitudinal 
instabilities  are  acceptable.  Often  a  simple  inspection  of  time  response  data  can  be 
made  to  check  the  criteria  in  'i’able  2.2.  It  is  then  not  necessary  to  find  the  phugoid 
parameters  through  a  LOPS  fit.  If  this  is  done,  it  suffices  to  use  second  order  LOES 
to  find  the  short  period  parameters.  The  first  order  numerator  and  second  order 
denominator  terms  corresponding  to  the  phugoid  mode  in  eqs.  (2. 16)  and  (2. 17)  can  be 
neglected . 

2.2.2  LATERAL/DIRECTIONAL  RESPONSE  REQUIREMENTS 

The  flying  qualities  for  the  luteral/directional  axes  arc  determined  by  the  response 
of  the  aircraft  to  a  lateral  stick  input  from  the  pilot.  The  military  standard  suggests 
using  roll  rate  or  roll  angle  as  the  output  for  a  LOEH  fit.  Because  we  are  concerned 
with  fighter  aircraft  which  may  maneuver  at  elevated  angles  of  attack,  an  extension 
to  the  standard  is  made  to  use  stability  axis  roll  rate  as  the  output.  The  low  order 
modal  form  required  for  the  LOES  fix  is  given  by 


^(s)  _  I^s(s^  + 

6i;,t,(.s)  (s  +  l/T'rlis  -h  lAl'y)  +  2QcO(jS  +  oi,^) 

where  p  is  the  stability  axis  roll  rate,  p  =  p  cos  a  4-  r  sin  (x,  and  6i.j(,  is  a  lateral  stick 
input.  A  three  step  process  may  be  taken  in  finding  a  low  order  fit  for  eq.  (2.19). 
First,  the  high  order  transfer  function  from  lateral  stick  input  lo  stability  axis  roll  rate 
is  used  to  find  a  first  order  roll  mode  approximation 


p(s)  _  Kj,  e-^iri 
Si.,i(.s) 


(2.20) 
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Next  the  high  order  transfer  function  from  rudder  pedal  input  to  sideslip  is  used 
to  find  &  second  order  Dutch  roll  approximation 


■  ■’! 


JKsJ_.  ^  (2.21) 

The  roll  mode  time  constant  found  from  the  fit  of  eq.  (2.20)  and  the  Dutch  roll 
frequency  and  damping  found  from  the  fit  of  eq.  (2.21)  are  put  into  eq.  (2.19)  and  fixed. 
The  fourth  order  DOES  in  eq.  (2. 19)  is  then  fit  to  the  HOS  to  find  the  equivalent  time 
delay,  numerator  poles,  and  spiral  mode  time  constant. 

The  general  requirement  for  equivalent  roll  mode  time  constant  must  be  less  than 
or  equal  to  1.0  second  for  Level  1  and  1.4  seconds  for  Level  2.  These  values  are 
conservative  for  modem  fighter  aircraft  which  typically  have  values  of  T,.  ranging 
from  1/3  to  1/2  second  for  conventional  flight  conditions.  At  high  angles  of  attack,  roll 
requirements  are  highly  dependent  on  control  power  limitations,  but  typically  time 
constants  are  slower. 

Requirements  for  Dutch  roll  frequency  and  damping  are  given  in  Table  2.3.  The 
high  Level  1  value  for  Dutch  roll  frequency  and  damping  i.s  driven  by  f  r.  '  and 
pointing  requirements  in  the  air-to-air  and  air-to-ground  missions. 

Table  2.3  Dutch  Roll  Frequency  and  Damping  Requirements 


Level 

Min  Cd 

Min  QjtOd  (rad/s) 

Min  Oj  (ra 

1 

0.4 

0.4 

1.0 

2 

0.02 

0.03 

0.4 

3 

0 

... 

0.4 

Equivalent  roll  time  delay  requirements  are  identical  to  the  pitch  requirements  in 
Table  2.1.  The  spiral  mode  typically  has  a  slov/  first  order  response.  For  a  manual 
flight  control  problem,  the  spiral  mode  is  generally  not  a  problem,  even  if  it  is 
unstable.  The  requirements  in  Table  2.4  define  the  speed  of  divergence  that  is 
tolerable  in  an  unstable  spiral  mode. 
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Table  2.4  Spiral  Mode  Time  to  Double  Amplitude  Requirements 


A  number  of  specific  requirements  have  been  described  for  Class  IV  aircraft  in 
Category  A  flight  phases.  These  specifications  are  derived  from  the  guidelines  in 
MIL-STD-1797A  and  by  no  means  constitute  the  complete  story  on  flying  qualities. 
For  more  detailed  information,  consult  the  standard  and  the  references  therein  [2.4]. 

2,3  UNCERTAINTY  REPRESENTATION 

It  is  important  to  remember  that  any  model,  no  matter  how  accurate,  can  only 
approximate  the  true  behavior  of  a  sy.stem.  In  fact  the  primary  purpose  of  feedback 
is  to  reduce  the  effect  of  uncertainty.  Robustness  is  a  measure  of  how  tolerant  a 
system  is  to  some  level  of  uncertainty,  either  structured  or  unstructured. 

Unstructured  uncertainty  is  characterized  as  additional  dynamics  that  are  not 
represented  in  the  analysis  model.  Classical  single  loop  gain  and  phase  margins  arc 
one  way  of  quantifying  robustness  to  unstructured  uncertainty.  For  multi-loop 
systems  a  more  general  framework  is  necessary.  Unstructured  uncertainty  may  be 
captured  as  either  an  additive  or  a  multiplicative  perturbation.  Let  Go(s)  be  the 
nominal  plant  and  G(s)  be  the  true  plant.  A.,  represents  additive  uncertainty  and 
represents  multiplicative  unc(*rtainty.  1'hen  for  an  additive  uncertainty 

G(s)  =  Go(.s)  +  AJs)  .  (2.22) 

For  a  multiplicative  uncertainty  at  the  plant  input 

G(s)  =  Go(s)(i  +  A„,(s}|  ,  (2.23) 

and  for  a  multiplicative  uncertainty  at  the  plant  output 


G(s)  =  [I  +  A,„(s)]Go(s) . 


(2.24) 


U  can  be  said  that  an  uncertainty  is  structured  if  the  perturbation  corresponds  to 
a  parameter  or  set  of  parameters  whose  variation  is  independent  of  the  rest  of  the 
system.  If  this  structure  can  be  captured  in  an  uncertaintv  model,  robustness 
results  will  be  less  conservative.  Structured  uncertainty  can  also  be  captured  as 
either  an  additive  or  a  multiplicative  perturbation.  Let  Iq  be  the  nominal  parameter 
value  or  set  of  parameter  values  and  /  be  the  true  value  or  set  of  values.  Then  using 
the  same  nomenclature  as  above 


Z  —  /q  + 

(2.25) 

1  —  Zq(  1  +  Ajj^) 

(2.26) 

The  structured  singular  value  (j-O  framework  described  in  section  2.4  provides  the 
analysis  tools  necessary  to  capture  the  robustness  of  a  system  to  any  of  the  above 
types  of  uncertainty,  structured  or  unstructured,  additive  or  multiplicative. 
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2.4  STRUCTURED  SINGULAR  VALUE  ANALYSIS 


Structured  singular  value  or  (i-analysis  is  a  framework,  based  on  the  small  gain 
theorem,  in  which  the  robustness  of  a  system  can  be  quantified.  Consider  the  closed 
loop  uncertain  system  in  Fig.  2.7, 


Fig.  2.7  M-A  Formulation 


Any  closed  loop  system  with  plant,  l^(s),  controller,  K(s),  and  uncertainty,  A{s),  can  be 
rearranged  in  what  is  called  M-A  form.  An  equivalent  representation  of  the 
uncertainty  block  is; 

A(s)  =  Wo  A(s)  Wi  ,  where  II  Afs)  lU  <  1.  (2.27) 

The  small  gain  theorem  guarantee.^?  that  if  M(s)  and  A(.s)  are  stable,  then  the 
uncertain  system  will  remain  stable  if  lor  all  frequencies,  0  <  O)  <  oo. 


o(  M(ja))  A(j{o) )  <  1  , 


(2.28) 


which  can  be  equivalently  expressed  as 


IIM(s)A(s)IL  <  1.  (2^9) 

An  inequality  can  now  be  invoked, 

II  M(s)  A(s)  IL  <  II  M(s)  IL  II  A(s)  IL. .  (2.30) 

Since  it  is  known  that  II  A(s)  IL  <  1,  a  sullicient  condition  for  stability  is; 

IIM(s)lL<  1.  (2.31) 

Eq.  (2.31)  can  be  shown  to  he  a  necessary  as  well  as  a  sullicient  condition  for  stability 
L2.51. 


When  the  uncertainty  model  is  highly  structured,  this  one  block  small  gain 
theorem  analysis  is  potentially  conservative.  The  general  perturbation  matrix  A(s) 
may  not  be  an  accurate  representation  of  the  true  uncertainty,  since  it  may 
overbound  the  true  uncertainty  and  represent  perturbations  to  the  nominal  plant 
that  are  unrealistic.  In  such  a  case,  eq.  (2.31)  is  only  a  sufficient  condition  for 
stability.  The  analysis  can  bo  made  less  conservative  by  considering  a  certain 
structure  for  the  uncertainty. 

Structured  robust  stability  assumes  that  the  perturbation  block,  A(.s).  ha.s  a 
particular  structure  based  on  specific  knowledge  of  how  the  uncertainty  enters  the 
plant.  A  general  model  of  the  uncertainty  is 

A=  {  diag(5ilrj,62lr^,...,hmIr,„-^l-^\2.--Ai)  >  bit-  C7,  Aj  G  CKr><.il  (2.32) 

=  {  A  e  A  I  a(A)  <  1  }  (2,33) 

where  A  is  the  set  of  all  possible  perturbations.  The  uncertainty  model  consi.sts  of  a 
diagonal  matrix  of  repeated  scalar  blocks  and  full  complex  uncertainty  blocks,  A 
single  scalar  perturbation  is  a  special  case  of  the  repeated  block  hmlr,,,-  where  r,„  =  1. 

Unstructured  uncertainty  is  represented  a.s  a  full  complex  block. 
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The  structured  singular  value  [2.6]  a  measure  of  robustness  to  complex 
perturbations  that  have  a  given  structure.  The  structured  singular  value,  p,  of  a 
complex  matrix,  M,  is  defined  as  the  inverse  of  the  maximum  singular  value  of  the 
sm2dlest  destabilizing  perturbation  that  has  the  specified  structure. 

-i—  =  min  {  a  (A)  I  det.(l  -  MA)  =  0  }  (2.34) 

|a(M)  Ac  BA 

If  M(s)  is  a  stable  closed  loop  transfer  matrix,  and  p(M(ja)))  is  evaluated  along  the 
imaginary  axis,  then  |i(M(ja)))  is  a  function  of  frequency  that  gives  the  size  of  the 
smallest  allowable  A  which  inovc.s  a  closed  loop  polo  to  the  imaginary  axis.  Since  the 
uncertainty  representation  is  scaled  to  be  less  than  one,  if  pdVKjo)))  is  less  than  one 
over  all  frequencies,  then  the  system  is  stable  for  all  possible  uncertainties  in  the 
allowed  set,  RA. 

While  in  the  general  case  |i(M)  cannot  be  calculated  exactly,  its  value  can  be 
placed  between  lower  and  upper  bounds.  We  can  define  matrices  U  and  D 

U  =  ldiag(Ui,  U2, !  UtUi  =  1|  (2.35) 

D  =  (diag(Di,  02-  dp  d2,  d,,)  I  D,  e  C^ixri,  1)^  e  >  0,  dj  €  94)  (2.36) 

In  [2,6],  the  following  properties  of  the  structured  singular  value  are  defined: 

(a)  p(aM)=  lalp(M) 

(b)  p(I)  -  1 

(c)  p(AB)  <  d(A'4i(B) 

(d)  |i(A)=d(A)  for  all  A  G  A 

(e)  UA  6  A  and  AU  g  A  for  all  A  g  A  and  U  e  U 

(f)  DAU'^  =  A  for  all  D  e  D  and  A  e  A 

(g)  |i(UM)  =  m(MU)  =  ,4(M)  for  all  U  6  U 

(h)  p(DMI)'M  =  pdVl)  for  all  I)  e  D 

(i)  p(M)  <  pdVI)  <  diM)  where  p(M)  denotes  the  spectral  radius  of  M 
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By  using  the  properties  (g),  (h),  and  fi),  the  following  bounds  on  |I(M)  are  derived. 


max  p(UM)  <  1J.(M)  <  inf  o(DMD-l)  (2.37) 

UeU  DeD 

It  can  be  proven  that  the  lower  bound  is  always  equal  to  pdVI),  but  the  maximization 
of  p(UM)  is  not  convex.  Local  maxima  can  occur,  making  a  global  solution  difficult. 

The  minimization  ofa(DMD-^)  is  convex.  While  this  upper  bound  is  an  equality  only 
for  n  ^  3,  computational  experience  has  shown  that  this  bound  provides  a  close 
estimate  to  the  actual  value  of  |i(M)  [2.6,2.71. 


In  analyzing  the  robustness  of  a  closed  loop  system,  it  is  useful  to  form  an  analysis 
model  of  the  form  in  Fig.  2.8. 


W, 


A 


w 


Fig.  2.8  Robust  Performance  Diagram 


is  the  weighted  plant  transfer  function.  M(s)  is  a  set  of  transfer  functions; 


■  Za  ■ 

[Mn 

Ml2l 

-  Wa 

z 

-  •^21 

M22  - 

w 

(2..38) 


It  has  been  shown  that  robust  stability  exists  if  and  only  if 

sup  |i(Mii(j(o))  <  1 .  (2.39) 

0) 
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It  is  often  possible  to  characterize  the  performance  of  a  system  in  terms  of  a 
weighted  closed  loop  sensitivitj'  function  or  model  following  error.  In  such  a  case,  a 
condition  for  closed  loop  performance  may  be  written  in  terms  of  an  infinity  norm 
bound.  The  transfer  function  between  w  and  z  in  Fig.  2.8  represents  this  weighted 
sensitivity  or  model  following  error.  The  matrices  are  scaled  such  that  nominal 
performance  exists  if 


I!  M22(s)  IL  <  1  (2.40) 

It  was  shown  in  [2.7]  that  by  creating  a  fictitious  uncertainty  block  between  w  and  z, 
the  performance  problem  can  be  put  into  the  framework  of  p-analysis.  This  fictitious 
uncertainty  is  represented  by  a  full  complex  block. 

It  can  be  said  that  nominal  performance  exists  if  and  only  if 

sup)I(M22(j<i>))^  1.  (2.41) 

(0 

By  combining  the  fictitious  performance  uncertainty  block  and  the  uncertainty 
blocks  for  stability  robustness  analysis,  a  robust  performance  analysis  model  can  bo 
formed.  If  the  structured  singular  value  of  this  combined  model  is  less  than  unity, 
then  the  performance  specification  of  cq.  (2.41)  will  be  satisfied  in  the  presence  of  the 
bounded  uncertainties. 

It  can  be  said  that  robust  performance  exists  if  and  only  if  [2.7] 

sup  ladVKjo)))  <  1 .  (2.42) 

(0 
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2.5  DYNAMIC  INVERSION 


The  purpose  of  dynamic  inversion  is  tx)  develop  a  feedback  control  law  that 
linearizes  the  plant  response  to  commands.  In  {.tenoral  the  nonlinear  aircraft 
dynamics  can  take  the  form 

X  =  flx,u),  y  =  Cx  (2.43) 

where  x  is  an  n-dimensional  state  vector,  u  is  a  m-dimensional  input  vector,  C  is  a 
pxn  matrix,  and  y  is  a  p-dimensional  vector  of  output  variables,  A  transformation  is 
necessary  to  put  the  equations  in  a  form  from  which  the  inverse  dynamics  can  be 
constructed.  Each  controlled  output,  yj,  is  dilferentiated  until  an  input  term  from  u 
appears  [2.8].  Only  m  outputs  can  be  controlled  independently  by  the  rn  available 
inputs,  therefore  p  must  equal  m.  As  shown  in  [2,9],  the  output  equations  may  now 
be  written  in  the  form.. 


yildil 

y,ld2l 


=  h(x^  +  G(x)u 


ypl^^pi  J 


(2.44) 


where  yjl^il  represents  the  djth  derivative  of  the  output  yj.  The  inverse  dynamics 
control  law  can  be  written  as 


u  =  G(x)  -h  V  -  h(x) ) 


(2.45) 


h(x.)  represents  the  nonlinear  output  dynamics  and  G(x)  represents  the  nonlinear 
t  ontrol  distribution.  The  parameter  v  represents  the  desired  linear  dynamics  of  the 
closed  loop  system.  With  the  inverse  dynamics  control  law  implemented,  the  closed 
loop  system  has  the  form. 


ylcll  V 


(2.46) 


If  the  system  is  observable  and 
p 


1’ 

Z  dj  =  n, then 
1=1 


placed.  If  Z  dj  <  n,  then  closed  loop  stability 
i=l 


all  of  the  closed  loop  poles  may  be 
cannot  bo  proven.  In  this  case  the 


28 


unobserved  dynamics  or  the  internal  dynamics  of  dynamic  inversion  must  be  checked 
at  local  operating  points  to  insure  stability  f2.8J. 

2.6  ROBUST  EIGENSTRUCTURE  ASSIGNMENT 

Eigenstructure  assignment  [2.10]  is  a  technique  that  uses  constant  output 
feedback  gains  to  arbitrarily  place  the  eigenvalues  and  eigenvectors  of  the  closed  loop 
system.  The  number  of  eigenvalues  and  eigenvectors  that  can  be  placed  arbitrarily 
depends  upon  the  number  of  plant  inputs  and  outputs.  Assume  a  linear  time- 
invariant  system  in  state  space  form 


X  =  Ax  +  Hu  (2.47) 

y  =  Cx 

fiere,  the  number  of  states  is  n,  the  number  of  inputs  is  ni,  and  the  number  of 
outputs  is  r.  The  number  of  eigenvalues  that  can  arbitrarily  be  placed  is  equal  to  the 
number  of  outputs,  r.  The  number  of  elements  of  the  eigenvectors,  corresponding  to 
the  r  eigenvalues  that  can  be  placed,  is  equal  to  the  number  of  inputs,  m. 

The  eigenstructure  assignment  problem  s(.atement  is  as  follows.  Given  a  set  of 
desired  eigenvalues  and  eigenvectors  vi*^,  find  a  real  feedback  matrix  F  of 
dimension  m  by  r  such  that  the  closed  loop  eigenvalues  and  eigenvectors  of  (A  -r  BFC) 
are  close  to  the  desired  ones.  The  gain  matrix  is  found  by  first  reordering  each  of  d  u;  r 
desired  eigenvectors  into  m  specified  and  n-in  unspecified  components  using  a 
reordering  operator,  Ri 


jvjdjRi  ^ 


ro  AQ\ 

Viy*‘xo/ 


where  li*^  are  the  ni  components  of  the  desired  eigenvector  that  are  specified,  and  di 
are  the  n-m  components  that  are  unspecified.  Gompute  the  achievable  eigenvector 
for  each  defined  eigenvalue  by  defining 

-  (^i^l  -  A)-lB  (2.49) 

and  reorder  the  achievable  eigenvector  using  the  corresponding  reordering  operator 
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(2.60) 


where  Li  is  the  specified  part  of  the  achievable  eigenvector  and  Di  is  the  unspecified 
part.  Now,  the  projection  of  the  desired  eigenvector  onto  the  achievable  subspace,  zi, 
is  found  by  minimizing 


J  =  11  lid  -  Liziir^ 

The  value  of  zj  is  computed  as 

Zi  =  (Li'^UylL^lid 

and  then  the  achievable  eigenvectors  are  found  as 

vi^  =  Li'/.i 


The  system  is  transformed  so  the  B  matrix  is  of  the  form 


A  = 


=  Q-lB 


Q-lAQ 


(~)  CQ 

X  =  CJx 


X\  -  Xi 
vi‘i  =  Qvi^ 


(2.51) 


(2.52) 


(2.53) 


(2.54) 


where  Q  is  a  similarity  transfonnation.  The  rows  of  vi^  and  A  are  partitioned  into  the 
iirst  M  rows  and  tlie  last  n-tn  rows 


Si 

^  1 

Ai' 

’  ^  = 

-Wi, 

-A2. 

(2.55) 


;i() 


Define 


S  =  [Xi^si  X2^S2  ...•  ?^r^Sr]  (2.56) 

V  =  [vi^  v2^  ....  vr®]  (2.67) 

Finally,  the  feedback  gain  matrix  is  computed  using: 

F  =  (S  -  AiV)(CV)-^  (2.58) 

2.7  FULL  ORDER  H«.  DESIGN 

The  general  Hoo  output  feedback  problem  can  be  characterized  simply  as  finding 

a  controller  K,  if  one  exists,  that  internally  stabilizes  the  closed  loop  system  and 
satisfies  the  condition  llT.^^li„o  <  y.  Fig.  2.9  illustrates  the  problem. 


G{s)  represents  the  design  model  including  the  plant  dynamics  and  weighting  filters. 


G(s)  = 


'_A I  B ' 
Gilf 


-  A 

Bv  1 

7^  - 

LC2 

^21 

(2.69) 


The  well  known  state  .space  lU,  theory  [2.1 1,2.12]  gives  the  equations  to  produce  a 

stabilizing  feedback  controller  I'rom  output  vector  y  to  input  vector  u  that  minimizes 
the  induced  Hoo  norm  between  w  and  z.  (Fs)  must  satisfy  the  follovdng  conditions: 


31 


(1)  (A,  B2)  is  stabilizable  aiid  (A,C2)  is  detectable 

(2)  Di2  is  fi-dl  column  rank  and  D21  is  full  row  rank 

(3)  for  all  0),  has  full  column  rank 

L  '^1  *Jl2. 

(4)  for  all  0),  has  full  row  rank 

L  ^2  b)2i 

The  first  condition  ensures  the  existence  of  a  stabilizing  controller.  Condition  (2) 
can  be  considered  analogous  to  the  nonsingular  control  penalty  and  nonsingular 
sensor  noise  weight  requirements  in  linear  quadratic  Gaussian  (LQG)  control.  The 
column  and  row  rank  conditions  of  (3)  and  (4)  are  made  to  assure  that  the  transfer 
matrices  from  w  to  y  and  from  u  to  z  have  no  invariant  zeros  on  the  imaginary  axis. 
Also  assume  for  convenience  that  Du  anci  zero,  and  that  the  system  is 

scaled  so: 

'Hi]  [O' 

Di2'r[Ci  D12I  =  [0  I];  D2i'1’=:  (2.60) 

LB21J  LI.I 

These  last  assumptions  arc  made  only  to  simplify  the  solution,  and  can  be 
removed.  The  plant  can  be  manipulated  using  the  approach  of  L2.13J  to  tran.sforin  an 
arbitrary  plant  to  one  that  meets  these  conditions.  Dcnne  two  Hamiltonian  matrices 
and  two  Riccati  equations; 

A  A’l'  -  C2'‘'G2' 

H«=  (2.61) 

.-Ci'i'Ci  -AT  J  LBjBiT  -a  J 

ATXoo  +  X<«A  +  X^(y-2BiBiT  -  B2B2T):X.<,  +  CiTCi  =  0  (2.62) 

AY.,.  +  Y.J^T  +  Y^{y  -20 1'l’C  1  -  C2TC2 )Y„,  +  B 1 B  iT  =  0  (2.63) 

Under  the  conditions  that  the  Hamiltonian  matrices  H,,,  and  Jc«  do  not  have  any 
eigenvalues  on  the  imaginary  axis,  X...  and  Yoo  arc  positive  sernidefinite  solutions  to 
the  two  Riccati  equations,  and  the  spectral  radius  ol'  the  product  of  the  Riccati 
solutions,  p(X«Y«.),  is  less  than  y2,  then  the  controller  which  satisfies  Y  i-"’ 

given  by: 
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(2.64) 


wheie: 


Xoo  —  AooXc 


vdooLooy 


u  -  FooX.., 


=  A  +  +  H2F00  +  Z00L00C2 

Fc«  =  -B2'‘Xo 
Loo  =  -YooC2'^' 

Z,„  =  (l-Y-^YooXoo)-l 


(2.65) 


The  number  of  controller  states  is  equal  to  the  order  of  the  design  model.  While 
closed  loop  stability  is  guaranteed,  controller  stability  is  not.  Many  physically 
motivated  problems  do  yield  a  stable  eo’nponsator. 

2.8  REDUCED  ORDER  OBSERVER  BASED  DESIGN 

A  minimal  order  dcsigit  algorithm  developed  in  [2.14]  considers  the  linear 
time-invariant  system  of  Fig.  2.U!. 


Fig.  2.10  Minimal  Order  Design  Model 


where  (A,B)  is  stabilizable,  (A,C)  and  (A,Fli)  are  detectable,  and  B  and  H2  have  full 
column  rank. 


The  minimal-order  Hoo  design  algorithm  generates  a  controller  that  stabilizes 
the  closed-loop  system  and  bounds  the  H„-norm  of  the  transfer  function  from  the 
disturbance  (w^)  to  the  controlled  outputs  by  y.  The  controller  has  the 

Luenberger  observer-based  structure  shown  in  Fig,  2.11  and  has  dimension  equal  to 
the  number  of  plant  states  minus  the  number  of  plant  measurements. 


F'ig.  2.11  Minimal  Order  Controller  Structm'c 


The  controller  state,  Xq,  is  an  estimate  of  a  linear  transformation  of  the  plant  state, 
Tx,  where  T  is  the  transformation  matrix.  The  controller  parainclcrs  Kj-,  F,  T,  M,  and 
N,  must  satisfy  the  Luenberger  constraints 

TA  -  FT  =  (2.66) 

G  =  'I’li  (2.67) 

NT  +  MC  =  K,„  (2.68) 


for  any  K<,  such  that  A-BKj,  is  Kurwitz. 

According  to  [2.14],  the  robust  controller  parameter  Kj,  is  given  by 

Kc  -  (2.69) 


v^'ith  Pj  >  0  satisfying 


PjA  -I-  a'^’Pj  h'J'Hi  (l+a)Y'^PjG,G}Vi 


P,b[hJhJ'‘b'^Pi  <  0  (2.70) 


.-{4 


subject  to  the  existence  of  a  P2  ^  0  that  satisfies 


P./  +  F’^i^2  +  n''h2H2N  +  0- 


The  design  parameters  a,  Kj-,  and  F  represent  the  freedom  available  to  the  designer 
and  are  chosen  such  that  a  is  a  positive  real  scalar,  F  is  any  stable  matrix,  and  Kf  is 
completely  arbitrai'y.  The  design  tuning  parameter  a  is  used  to  tighten  the  norm 
bound  between  the  disturbance  and  the  controlled  output  [2.15].  An  initial  bound  yis 
selected  that  allows  the  existence  of  a  positive  semi-definite  solution  to  eq.  (2.70).  The 
controller  parameter  K(.is  now  computed  using  Pj,  and  then  the  Luenberger 
constraints  in  cqs.  (2.66)-(2.(i8)  arc  solved.  Sylvester  equation  (2.66)  i.s  solved  for  T, 
and  then  G  is  computed  using  eq.  (2.67).  Parameters  N  and  M  are  computed  using  a 
variation  of  eq.  (2.68); 


IN  M1  =:  K 


T 

C 


-1 


(2.72) 


Note  that  the  choice  of  K,-  and  F  provide  the  design  freedom  to  obtain  a  T  that  insures 


the  existence  of 


T 

C 


.  11'  there  exists  a  solution  to  eq.  (2.71),  a  series  of  y-reduction 


iterations  is  performed  on  eq.  (2.70)  to  find  the  smallest  bound,  y,,.;,,  ,  for  which 
positive  semi-definite  solutions,  Pj  and  P^,  exist.  If  no  solution  exists  to  oq.  (2,71),  yis 
increased,  and  the  procedure  i.s  repeated. 


2.9  STRUCTURED  SINGULAR  VALUE  SYNTHESIS 


The  structured  singular  value  (p)  framework  provides  a  unifying  measure  which 
can  be  used  to  simultaneously  address  stability  and  performance  robustness 
specifications.  If  p  is  less  than  unity  for  a  properly  scaled  system,  then  the 
spocification.s  arc  met.  It  is  desirable  to  be  able  to  address  these  multiple  objectives 
directly  within  a  design  method,  p-.synthosi.s  provides  for  the  direct  incorporation  of 
robust  stability  and  performance  goals  into  a  design  by  combining  H.x,  design  with 
structured  singular  value  analysis  [‘2.16|.  M’he  p-synthesis  problem  is  described  by 
the  attempt  to  find  u  controller  that  minimizes  an  upper  hound  on  the  structured 
singular  value, 
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(2.73) 


min  inf  sup  C5(DM(K)l>l). 

K  De  D  (0 

M(K)  is  the  weighted  closed  loop  transfer  function  shown  in  Fig.  2.8.  One  approach  to 
this  problem  is  the  DK-iteration,  it  calls  for  alternately  minimizing  sup  a(nM(K)D  h 
for  either  K  or  1)  while  holding  the  other  constant.  First  the  controller  synthesis 
problem  is  solved  using  H«,  design  on  the  nominal  design  model,  P^Cs).  p-analysis  is 

then  performed  on  the  closed  loop  transfer  function  M(K),  producing  values  of  the  D 
scaling  matrices  at  each  frequency.  The  resulting  frequency  response  data  is  fit  with 
an  invertable,  stable,  minimum  phase  transfer  function  which  becomes  part  of  the 
nominal  sjmthesis  structure.  With  D  fixed,  the  controller  synthesis  problem  is  again 
solved  by  performing  an  Hoo  design  on  the  augmented  system.  The  DK-iterations  are 

continued  until  a  satisfactory  controller  is  found  or  a  mininium  is  reached.  Fig.  2.12 
shows  a  flow  diagram  for  the  DK-itoration  . 


Fig.  2.12  DK-iteration  Process 


The  resulting  controller  order  is  the  order  of  the  desi,..!  plant  and  weighting  matrices, 
as  well  as  the  order  of  the  D-scale  transfer  function  fits.  With  each  iteration,  the  D- 
scale  frequency  response  data  from  tlie  previous  iteration  is  eoinbined  with  the 


current  values,  and  then  the  transfer  function  fit  is  performed  on  the  combined  data. 
This  approach  avoids  a  built-in  increase  in  controller  order  that  would  result  if  at  each 
iteration  new  D-scale  fit  transfer  functions  were  au^'nented  into  the  synthesis  model 
from  *^he  previous  step.  It  is  important  to  note  that  the  DK-iteration  is  not 
guaranteed  to  converge  to  a  global  minimum,  but  practical  experience  has  shown 
that  the  method  works  well  for  a  broad  class  of  problems  [2.17], 

2.10  BALANCED  REALIZATIONS  AND  TRUNCATION 

If  a  system  G(s)  =  C(sl  -  AV^B  +  D  is  minimal  and  stable,  a  coordinate 
transformation  matrix  T  exists  such  that  the  transformed  system  has  controllability 
and  observability  grammians  that  are  equal  and  diagonal  [2.18], 

The  controllability  grammian,  F,  and  observability  grammian,  Q,  are  the 
hermitian  solutions  to  the  following  fiyapunov  equations. 

AF -I- PA'l' -1- BB'f  =  0  (2.74) 

A'l'Q  +  QA C'l'C  =  0  .  (2.75) 

The  Cholesky  factorization  of  Q  is 

Qi  =  li'l'H  .  (2.76) 

RPR'f  is  a  positive  definite  matrix, 

llltri’  =  RlMfl'  where  U'l’U  =  I , 

and  I  is  a  diagomd  matrix  of  the  1  lankel  singular  values  in  descending  order 

L  =  diag(ai,  0-2,  . .  a„),  ni  >  02  .  .  >  a„  >  0  (2.77) 

The  transformation  which  balances  the  system  is 

T  =  i-i/^u'rR . 
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(2.78) 


The  controllability  and  observability  grammians  of  the  transformed  system  are  equal 
and  diagonal  [2,19]. 


P  =  TPTT' =  I-l/^LrrRPRTui  1/2  =  x  (2.79) 

Q  =  (XT')  -lQT-1  =  rl/2iri'(R'r)  -hRT’RV’  -lUi:-l/2  =  I  .  (2.80) 

By  changing  the  state  space  coordinates  to  %  =  Tx,  then  the  transformed  system 
becomes 

z  =  A  z  +  Bu  y  =  Cz  +  Du  ,  (2.81) 

where 

A  =  TAT-l,  B  =  TB,  and  C=CT-1  ,  (2.82) 


Because  the  Hankel  singular  values  are  arranged  in  descending  order,  the  last 
elements  of  the  coordinate  vector  z  correspond  to  the  least  controllable/ohservable 
part  of  the  system.  The  system  in  eq.  (2.81)  can  be  partitioned  into  a  k  dimensional 
vector  Z2  and  an  n-k  dimensional  vector  Z]. 
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(2.83) 


(2.84) 


The  controllability  and  obsei’vability  grammians  of  the  partitioned  system  are 


P  =  Q  = 


Xu  0 
0  X22 


(2.85) 


A  reduced  order  model  may  now  he  formed  by  removing  the  partition  7.2  from  the 
state  vector  z.  This  approach  is  known  as  balanced  truncation.  The  reduced  order 
system  has  controllability  and  observability  grammian.s  of  Xu>  that  is  the  new 


system 
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7,1=A|iZ]+  Bju  y  =  Cjz^  +  Du  ,  (2.86) 


satisfies  the  Lyapunov  equations 

Ai  jIii  +  IjiAji'*'  +  =  0 

(2.87) 

An’l'l„  +  IiiAii+CiTCi=0  . 

(2.88) 

The  reduced  order  system  Or  =  Ci{sl  -  Anh^Bl  +  B  preserves  the  stability  of  the 
original  system.  The  error  between  the  original  and  reduced  order  system  has  an  Loo 

en'or  bound  given  by  [2. 19] 

!ICi(j6))  -  (ir(j(o)  ll|^^  ■^^^n-k+1'^  ^n-k+2'^  .  .  -  +  Cfn)  (2.89) 

The  direct  feedthrough  term  of  the  reduced  order  system  does  not  affect  the 

A. 

grammians,  so  the  D  term  can  be  adjusted  to  improve  the  approximation.  To  force  an 
exact  match  at  high  frequencies,  set  Or^w)  =  (l(j(o)  with  (o  =  0°.  Then 

f)  =  l)  .  (2.90) 

To  force  an  exact  match  in  steady  state  gain,  set  (lr(jw)  --  G(ja))  with  o)  =  0.  Then 

D=  C(-A)1B+  D  -  C-,(-  A]  j)-iBi  .  (2.91) 

To  approximate  a  match  at  a  single  frequency  or  set  of  frequencies,  (Op  0)2,  .  .  .  O),^, 
the  direct  feedthrough  term  is  rhoson  to  .solve 

m 

IIG(j0),)-(v(j(0i)l|2  ,  (2.92) 

D  i=i 

which  has  the  solution 

1>=  ReiGljco.l  -  A)  IBs-  I)  -  flpjo),!  -  AjD-^T^)  ■ 

i-l 


(2.93) 


Other  variations  of  balanced  trucation  include  optimal  Hankel  norm  approximations 
[2.19]  and  frequency  weighted  balanced  truncation  [2.20], 
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3  CONTROL  DESIGN  METHODOLOGY 


This  section  details  a  design  methodology  for  manual  flight  control  systems.  The 
approach  is  based  on  an  inner/outer  loop  control  structure.  Control  redistribution  is 
governed  by  a  control  selector.  The  control  selector  is  based  on  generalized  inverses 
which  normalize  the  control  effectiveness  with  respect  to  generalized  inputs.  An 
inner  loop  is  formulated  that  equalizes  the  plant  dynamics  across  the  flight  envelope. 
This  inner  loop  is  a  static  output  feedback  or  low  order  compensator  that  is  a  function 
of  flight  parameters  such  as  altitude,  Mach  number,  and  angle  of  attack.  The  outer 
loop  controller  is  designed  using  advanced  techniques  to  achieve  performance  and 
robustness  goals.  This  controller  should  not  require  gain  scheduling  because  of  the 
equalization  performed  by  the  inner  loop. 

3.1  CONTROL  SELECTOR 


The  control  selector,  sometimes  referred  to  as  pseudo-controls,  has  two  functions. 
The  first  is  to  normalize  control  effectiveness  by  transforming  generalized  rotational 
rate  commands  into  actuator  position  commands.  The  second  is  to  take  advantage  of 
available  control  redundancy  by  allowing  for  control  redistribution  without  changing 
the  linear  dosed  loop  pcrlbrmance.  The  basic  idea  of  the  control  selector  is  to  redefine 
the  control  contribution  to  the  state  equation  13.1,3.2], 

H5  =  B’6*.  (3.1) 


where  B  is  the  actual  control  effectiveness  matrix,  5  is  the  vector  of  control 
deflections,  11“  i.s  the  generalized  control  effectiveness  matrix,  and  6“  is  the  generalized 
control  vector.  The  iictual  control  can  now  be  defined  in  terms  of  the  generalized 
control. 


5=:T8'. 

'fhe  transformation,  'f,  is  the  control  selector,  it  i.s  defined  simply  by 

T  =  NlBNl^^B'- 


(3.2) 
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(3.3) 


The  operation  (  )^  is  a  pseudo-inverse  and  N  is  a  matrix  that  may  be  used  to  combine 
controls  or  emphasize/de-emphasize  a  control  channel  in  the  case  of  redundant 
effectors.  Because  the  B  matrix  in  eq.  (3.3)  is  a  function  of  flight  condition  and 
aircraft  state,  the  control  selector  is  a  function  of  parameters  such  as  Mach  number, 
altitude,  and  angle  of  attack. 

3.2  INNER  EQUALIZATION  LOOP 

An  inner  equalization  loop  is  used  to  account  for  the  changes  in  plant  dynamics 
with  flight  condition.  The  goal  is  to  make  the  input/output  behavior  of  the  closed  loop 
system  uniform  for  all  operating  conditions  by  using  a  nonlinear  static  feedback 
matrix  .  This  inner  loop  feedback  can  be  derived  by  scheduling  linear  compensators  or 
through  direct  nonlinear  methods  such  as  dynamic  inversion. 

This  equalization  can  be  quantified  by  the  concept  of  relative  error.  Relative  error 
is  defined  as 


Am  ;=  (P  -  Po)Po^^  (3.4) 

where  Pq  is  the  plant  used  for  outer  loop  design  and  P  is  the  equalized  plant.  This 
idea  of  relative  error  gives  us  a  mathematical  framework  for  analyzing  the 
effectiveness  of  an  inner  equalization  loop.  A  Robustnr.ss  Theorem  by  Safonov  and 
Chiang  can  provide  a  weak  sufficient  condition  for  stability  based  on  relative  error 
[3.3]: 


If  a(Ajn)  <  1  for  0)  <  (0^,  then  the  dosed  loop  -system  will  be  stable  provided  that 
the  control  bandwidth,  o)ij,  is  loss  than  0^. 


Thus,  if  a  system  with  an  outer  loop  controller  designed  fur  Pq  i.s  closedi  loop  stable. 


then  closed  loop  stability  should  be  preserved  if 
relative  error  is  -suflicicntly  small. 


is  replaced  by  P  provided  that  tVio 


Explanation:  Consider  the  system  with  equalized  plant,  P,  outer  loop  controller,  Kf^i, 
and  relative  error,  A^-  The  bandwidth  of  the  control  system  is  defined  as  the 
frequency  range  where  the  loop  transfer  function  gain  is  big  13.3).  That  is. 
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(3.5) 


a(PKoi)  »  1  V  to<  tojj 
A  sufBcient  condition  for  stability  is 

^A„i)  a(PKoi(I  +  PK„])-1)  <  1  (3.6) 

It  follows  from  eq.  (3.5)  that  for  (o  < 

n(PK„,(I  + PKoj)!)-  1  (3.7) 

So  for  frequencie.s  whoro  tho  loop  transfer  function  f^ain  is  hig,  a  sufTicicnt  condition 
for  stability  i.s 

d(A,„)  <  1  (3.6) 

This  explanation  is  not  a  proof  in  a  mathematical  sense  because  of  the  loose  definition 
of  what  is  considered  big.  Strong  claims  .should  not  be  made  about  robustness  at 
frequencies  where  the  loop  gain  i.s  not  very  large.  In  addition,  the  theorem  says 
nothing  about  what  happen.s  at  frequencies  above  the  control  bandwidth  where  the 
relative  error  may  be  greater  than  unity.  At  the.se  frequencies  there  is  no  guarantee 
that  eq.  (3.6)  will  be  satisfied,  It  is  also  important  to  remember  that  eq,  (3.6)  is  only  a 
sufficient  condition  for  stability  which  can  be  very  conservative.  Considering  all  of 
these  factors,  the  primary  utility  of  the  Robust nv-ss  Theorem  for  this  work  is  that  it 
provides  a  framework  for  quantifying  whether  or  not  an  inner  loop  has  been 
successful  in  equalizing  a  sot  of  plants.  More  rigorous  post-design  analysis  is 
petdormed  to  provide  well  defined  bounds  for  stability  and  performance  robustness. 


Fig.  3.1  shows  the  implementation  of’ the  equalization  loop  using  generalized 
controls. 


P 


CS  is  the  control  selector  and  is  the  equalizing  controller.  With  the  equalization 
loop  in  place,  P  should  approximate  a  linear  .system  across  the  flight  envelope. 

3.3  OUTEK  ROBUST  PERFORMANCE  LOOP 


The  ouior  robust  performance  loop  uses  a  fixed  output  feedback  compensator  to 
achieve  performance  and  robustness  to  parametric  uncertainties,  The  inner 
equalization  loop  along  with  the  control  selector  lineahze  the  aircraft  response  to 
generalized  commands.  using  these  generalized  eommand.s  as  control  inputs,  a 
single  output  feedback  compensator  can  satisfy  design  goals  for  a  broad  range  of 
operatii'.g  conditions.  The  intractable  problem  of  gain  scheduling  a  high  order 
compensator  is  eliminated.  Any  linear  de.sign  method  can  now  he  used  to  achieve 
robustness  and  performance  goals.  Fig.  3.2  shows  the  outer  loop  implementation. 
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4  VISTA  F-16  LATERAL/DIRECTIONAL  DESIGN 


The  following  sections  describe  the  design  and  analysis  of  a  full  envelope,  manual 
flight  control  system  for  the  lateral/directional  axes  of  the  Variable  Stability  In-b’light 
Simulator  Test  Aircraft  (VISTA)  F-16  test  aircraft.  This  example  shows  tlie  utility  of 
the  inner/outcr  loop  control  structure  as  well  as  how  various  state-of-the-art 
technologies  can  be  brought  together  in  a  comprehensive  design.  First,  the  aircraft 
and  specific  design  requirements  are  described.  Next,  (he  design  approach  and  results 
are  presented,  including  linear  and  nonlinear  results.  Finally,  the  linear  aircraft 
models  and  resulting  controllers  iire  given  in  an  appendix  so  that  these  results  may  be 
verified,  duplicated,  or  compared  to  other  approaches. 

4.1  MODEL  DESCRIPTION 

VISTA  is  an  advanced  development  program  in  tlie  l'’Iigiit  Dynamics  Directorate  of 
Wright  Laboratory.  VIS'l'A  is  a  modified  F-16  with  the  capability  to  simulate 
advanced  aircraft  configurations  and  test  advanced  flight  control  concepts.  Fig.  4.1 
shows  a  three-view  of  the  F-16  with  AIM-9  heat  seeking  air-to-air  missiles  on  the 
wingtip.s.  The  aircraft  has  leading  and  trailing  edge  flaps,  an  all  moveable  horizontal 
tail,  and  a  single  rudder  1'ho  pilot  controls  include  a  force-feel  side  stick,  rudder 
pedals,  and  throttle. 

A  high  fidelity,  six  degree  of  freedom,  nonlinear  simulation  mode!  lias  been 
developed  for  the  VIH'l'A  F-16  vehicle.  The  Vlri'l’A  nonlinear  model  is  written  as  a 
series  of  FORTRAN  subroutines  and  is  implerrented  and  validated  in  a  generic 
nonlinear  simulation  environment.  Specific  application  modules  include  .subroutines 
describing  the  equations  of  motion,  actuators,  and  sensors.  The  actuation  models  for 
all  of  the  control  surface  deflections  include  fourth  order  linear  dynamics,  hinge 
moment  effects,  and  position  and  rate  limits.  Additional  subroutines  doscrihe  the 
propulsion  system,  weight  and  moments  of  inerfda  variations,  and  the  atmospliere. 
The  aerodynamic  data  exist  for  a  wide  range  of  Mach  numhors,  altitudes,  and  angles 
of  attack  and  sideslip.  The  high  fidelity  model  allows  for  the  simulation  of  long 
duration  and  large  amplitude  maneuvers  with  extreme  accuracy,  'fhe  nonlinear 
model  is  used  to  generate  linear  models  for  control  law  design  and  (.o  generate 
nonlinear  time  histories  to  evaluate  control  designs  A  diagram  of  llu'  VlS'l’A  I*'- 16 
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4.2  Nonlinear  Model  Block  nia^rrain 
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A  trim  solution  can  be  derived  for  the  aircraft  at  a  spL'cifu:  Mach  number  and 
altitude  at  which  all  linear  and  angular  accelerations  are  zero.  'I'he  solution  to  oq. 
(2.4)  is  found  through  a  numerical  routine  implemented  in  tht;  generic  simulation 
environment.  Once  a  trim  condition  is  established  for  the  nonlinear  aircraft  model,  a 
linear  model  can  be  generated  to  capture  the  perturbational  dynamics  around  the 
equilibrium  point,  binear  models  are  generateci  I'or  the  o|)en  loop  aircrall  using  two- 
sided  perturbations  to  evaluate  the  necessary  partial  derivatives.  M'ho  linearization 
procedure  yields  a  tenth  order  model  of  op<m  loop  aircraft  dynamics.  The  mode'  state.s 
are  transformed  into  conventional  aircrall  state.s:  Buler  angles,  body  axi.s  rotational 
rates,  forward  speed,  angles  of  attack  and  sideslip,  ai>d  altitude.  For  a  wings  level 
trimmed  condition,  the  linear  model  consists  of  decoupled  11  (Ih  o*'der  longitudinal  and 
1 I I  •)<  <>»*•  s  I //li  •!  '» 1  ii.? 
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by  removing  the  trajectory  states:  the  neutrally  stable  heading  angle  and  the  roll 
angle.  The  lateral/directional  linear  model  becomes: 
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^  ■■  1 
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=  Y|5j^  +  sincx  p  -  com.  r  -k  Ysiri'^nT  +  YsdfSdf  +  Y^rSk 
p  =  L|j(^  +  LpP  +  L.r  +  Lfiiri'^irr  -+  LfiDF'^UF  + 
r  =  N|j[3  +  Npp  +  N,.r  +  Nfii/ivSi)'!’  +  Nfii)i.’8[)i.'+  N§r6r  ; 

wiioro  (5  is  anf^lc  ol’ sideslip,  a  is  aiiRlo  of  attack,  p  is  body  axis  roll  rate,  r  is  body  axis 
yaw  rate,  bin'  is  din'erontial  liori/.ontal  tail  dellection,  8i>i.’  is  difi’orential  flap  deflection, 
and  5k  is  mdder  deflection. 

The  measured  outputs  are  roll  and  yaw  rates,  aiiRlc  of' attack,  and  sideslip  angle.  A 
stability  axis  ixjll  rate,  p,  out  put  (■•an  h(‘  constructed  from  these  measurement, s 


!>  cosex  +  r  simx . 


(4.2) 


It  is  assumed  for  this  dcssign  tint  angh?  of  attack  and  sideslip  angle  can  either  bo 
measured  directly  or  reconstructed  from  inertial  data.  Only  one  weight  and  store 
configuration  and  the  trim  throttle  setting  are  used  in  the  design  and  analysis,  Stores 
con.sist  of  two  AlM-9  missile.s  on  the  wing  tips. 


4.2  FLYINCJ  QUAI^rriKS  REQUIREMENTS 

A  high  order  transfer  function  which  includes  controller  dynamics  and  actuator 
dynamics  is  approximated  across  a  frequency  range  of  interest  by  a  low  order  fit.  The 
parauu'ters  of  tlu!  low  order  e(|uivalent  system  are  used  to  determine  whether  the 
aircraft,  respon.se  will  be  acci'f)tal)le  t.o  a  pilot.  The  following  requirements,  taken  from 
section  2.2,  must  be  met  for  hevol  1  flying  cjualities:  Dutch  roll  frequency  ,  0)j)  >  1 
rad''s,  Dutch  roll  damping,  >  O.'l,  roll  mode  time  constant,  'I'ii  <  1.0  sec,  spiral 
mode  stable  or  time  to  double  >  12  s(’c.  and  equivalent  time  delay,  ij,  <  0.10  sec.  For 
high  angle  of  attack  flight  the  criteria  for  Dutch  roll  frequency  is  relaxed  to  O)])  ^  0.7 
rad/s  and  roll  mode  time  constant  is  rtdaxed  to  Tg  <  l.C)  sec. 


4.3  CONTROL  SELECTOR  DESIGN 


The  generalized  inputH  lor  this  design  arc  body  axis  roll  acceleration  cornniand  and 
body  axis  yaw  acceleration  command.  The  true  control  inputs  are  asymmetric 
horizontal  tail,  asymmetric  flaps,  and  rudder.  'I’he  cotrespon.ding  control  cflectivencss 
matrices  are; 

Yhi/i'  Y^dp  Y*k  ■ 

W  ■■  I'.SDT  l-ftDF  hsK  (4.,^) 

_  N.sdt  Nsuk  Nsk  . 


(4.4) 


It  is  important  to  note  that  the  pseudo-inverse  in  the  control  selector,  (.'q.  (3,3),  does 
not  explicitly  accovint  for  limitation.s  m  control  deflection.  Undesirable  results  are 
therefore  possible  if  the  control  selector  transformations  ari'  applied  naively.  One 
example  is  canceling  roll  inputs  from  (he  asymmetric  flaps  and  a.symmotric 
horizontal  tail  that  together  achieve  a  very  small  increment  in  yaw.  Huch  cases  can 
catise  unreasonable  control  deflections.  If  N  in  eq.  (3.3)  is  chosen  to  be  id('ntity  in  the 
transformation  equations,  this  is  I'xacfly  what  happens  for  the  VISTA  I*'- lb.  This 
behavior  can  easily  be  prevented  by  combining  the  asymnudric  flaps  and 
asymmetric  horizontal  tail  into  a  single  control  efl'ector  which  we  will  call  aileron  f6^\). 

oa 

---  N  (4.5) 


The  matrix  N  fixes  the  proportion  between  asyinmetric  hoi'i/ontal  tail  ami 
asymmetric  flap  commands.  Because  the  primary  purpose  of  tho  horizontal  tail  is 
pitch  control,  a  ratio  of  1/4  is  used 


O.^if)  ().() 
1.0  ().<) 
DO  1.0 


(4.6) 


This  c-hoiL-i'  is  based  upon  enpniie(*riiipr  judgement  and  is  non-unique.  More  olefiant 
solutions  tiiat  niaxiiui/.e  contfol  elleel iveiiess  iiri'  possible.  A  lew  possibilities  are 
described  in  [4.1]. 

Tile  ciuitrol  selector  equation  is 
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d  ~  'fd'  where  iS 

and  d'  - 

Pc  1 
1 

-  bit  _ 
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and  the  ri'sultinp  iransf'onnation  is 

'r=N(BN)'^lb  (4.8) 

'I'he  control  derivativi  s  m  eq  '  (..qi  mx«  stored  in  a  tabular  database  as  a  lunction  of 
Mach,  altitude,  and  aiiftle  ot  attack  'Phe  control  selector  is  found  at  any  flight 
condition  by  performing  a  (able  lookup  of  these  parameters  and  using  eq.  (4.8)  to 
calculate  '!' 

4.4  INNKR  LOOP  DESIGN 

1  he  goal  ot  the  inner  equali/ation  loop,  as  described  in  section  8.2,  is  to  linearize  the 
mput/output  behavior  of  the  .system  across  the  flight  envelope.  The  dominant 
nonlineanties  which  must  he  handled  are  the  changes  in  the  aerodynamic  stability 
and  control  derivatives  with  flight  condition.  In  order  to  achieve  our  goal  of  plant 
eciuali/.ation,  a  controller  iiuir.t  he  found  which  itself  is  a  function  of  flight  condition. 

4.4.1  INNER  LOOP  FORMUlJVTION 

Dynamic  inversion  is  used  to  develop  the  inner  loop  equalization  control  law.  For 
inner  loop  synthesis  we  uichule  only  the  iioiilinearities  associated  with  changes  in  the 
aerodynamic  parameti'i  s  as  a  f  unctiou  of  Mach  (M),  altitude  (h),  and  angle  of  at  tack 
b/.i.  With  (lie  control  selector  implementi'd  and  neglecting  actuator  dynamics,  the 
aircraft  lateral  equations  ofniotion  are. 
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Because  the  control  distribution  matrix  has  a  column  rank  of  2,  only  two  output? 
may  be  chosen  for  the  dynamic  inversion  formulation.  'I'hc  body  axis  roll  and  yaw 
rate  states  are  chosen  as  the  controlled  outputs  as  they  represent  the  dominant  fast 
dynamics  of  the  open  loop  system.  Those  variables  must  only  be  differentiated  once 
for  the  control  to  appear  in  the  output  equations. 
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The  inverse  dynamics  control  law  can  now  be  written  as 
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(4.10) 


(4.11) 


where  V  is  a  matrix  that  represents  the  desired  linear  dynamics.  Note  that  the  f) 
portion  of  the  dynamics  in  eq.  (4.9)  is  not  ob.served  by  cq.  (4. 10).  The  inner  loop 
equalization  will  therefore  not  he  perfect  and  the  internal  dynamics  of  dynamic 
inversion  will  show  up  as  inner  loop  relative  error. 


The  choice  of  the  desired  dynamics,  v,  i.s  a  critical  step  m  the  formulation  of  the 
inner  loop.  If  unreasonable  parameters  are  chosen,  the  control  system  will  bo 
sensitive  to  unmodeled  dynamics  and  have  a  tendency  towards  actuator  rate  and 
position  limit  saturation.  For  the  VISTA  F-16  design,  a  linear  quadratic  regulator 
design  is  performed  at  a  flight  condition  which  is  considered  c.e.ntroL.  Central  simply 
means  that  the  chosen  flight  condition  is  at  some  dynamic  pressure  between  the 
minimum  and  maximum  value.s  for  the  design  envelope.  This  choice  is  based  upon 
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engineering  judpnent  and  is  made  such  that  the  relative  error  between  the  open  loop 
central  model  and  all  other  models  is  reasonably  small.  The  regulated  dynamics  at 
this  condition  become  the  desired  dynamics  for  the  dynamic  inversion  calculations: 


Aiom  “  ^\vntral  '  ^ 


(4.12) 
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(4.13) 


The  inner  equalization  loop  can  be  represented  as  a  linear  state  feedback, 
compensator  of  the  form: 


Pc 
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~  - 

_  r  _ 

(4.14) 


where 


^noni21'^''p^^  '^iioin23’^ 
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(4.15) 


The  aerodynamic  parameters  in  eq.  (4.15)  are  stored  in  a  tabular  database.  The 
inner  loop  feedback  is  derived  at  any  flight  condition  by  performing  a  table  lookup  of 
these  parameters  and  calculating  the  gain  matrix,  The  equalized  system  is 

shown  in  Fig.  4.3. 


Fig.  4. .‘I  Equalized  SyKtom 


Relative  ei*ror  must  bo  calculated  for  the  system  which  is  to  be  regulated  by  the  outer 
loop.  The  transformation  matrices  Tj  and  T2  are  included  so  that  tlie  ouiei  loop 
provides  regulation  in  dit'  stability  axis  rather  than  in  the  body  axis.  The  maLrices 
are: 


T,  - 


siacx  COSO, 
-cosot  sir.cK 


i  1  0  0 

_  I 

^  L  0  eos(x  siiKX 


(4.16) 


By  regulating  in  the  stability  axis,  tlie  control  has  a  built  in  dependence  on  angle  of 
attack  that  will  alleviate  some  of 'he  nei;d  for  sclKiduling  with  (x.  Tj  is  a  unitary 
matrix  and  therefore  does  nut  allect  singular  vah.es. 


The  two  transfer  funxhions  of  copcern  for  the  relative'  error  analysi.s  are  ^hPo'l'j 
and  T2  PT].  is  (he  .system  inodel  that  is  usi;d  for  the  outer  loop  des'gn.  it  is  the 
regulated  central  system  with  third  order  actuator  dyiiamics  included  in  the  iunen 
loop,  P  is  the  equalized  plant  which  is  formed  at  each  analysis  point.  'I'he  internal 


dynamics  of  dynamic  inversion  and  the  effects  of  actuator  dynamics  arc  the  tv/o 
primary  sources  of  relative  en'or  between  these  transfer  functions. 


Two  different  sets  of  desired  dynamics  are  used  in  this  application.  One  is 
formulated  from  a  central  model  at  a  low  angle  of  attack  (low-cx)  flight  condition  and 
one  is  formulated  from  a  central  model  at  a  high  angle  of  attack  (high-a)  flight 
condition.  These  desired  dynamics  reflect  the  different  robustness  and  performance 
requirements  at  high  and  low  angles  of  attack. 

4.4,2  LOW  ANGLE  OF  AITACK  INNER  LOOP  RESULTS 

The  central  point  for  the  low  angle  of  attack  de.sign  is  the  trim  flight  condition  at 
Mach  0.6,  altitude  20, 000  feet.  The  dynamic  pressun'  at  this  point  is  245.1  psf. 
is  the  solution  to  the  linear  quadratic  regulator  problem  with  identity  state  and 
control  weightings.  Details  on  this  design  are  presented  in  Appendix  4. 

Fig,  4.4  shows  the  relative  error  results  for  four  test  points,  'fho  test  points  shown 
in  Table  4,1  represent  a  wide  range  of  operating  conditions. 

Table  4, 1  Test  Points  for  Low-cx  Linear  Analysis 


Test  i’oint 

Mach 

Altitude  (ft) 

(X  (deg) 

q  (psO 

Central 

0.6 

20,000 

4.  .‘3 

245  1 

1 

0.4 

25,000 

12.6 

87.0 

2 

0.5 

20,000 

6.2 

170.2 

;i 

0,6 

15,000 

5,6 

501,1 

A 

J 

0.7 

10,000 

2.0 

49S.2 

I  rcqucncy  (rad/scc) 

Fig.  4.4  Relativo  Error  Rosull.^  Low  a 

While  relative  error  is  small  (less  than  db)  lor  the  frequency  range  of  interest 
(<  10  rad/s),  it  is  not  less  than  unity  for  all  conditions  tested.  'Phis  is  primarily  the 
result  of  interactions  between  the  equalization  loop  and  actuator  dynamics.  These 
interactions  are  not  included  in  the  dynamic  inversion  calculations.  Relative  error 
less  than  unity  is  only  a  suHicient  condition  for  closed  loop  stability.  Because  this  is  a 
conservative  test,  errors  that  are  small  but  greater  than  unity  can  bo  tolerated  with 
satisfactory  stability  and  perfornianco  results,  'fho  errors  shown  in  Fig.  4,4  are 
considered  small  for  frequencies  of  interest. 

4.4.3  HIGH  ANGLE  OF  ATTACK  INNER  LOOP  RESULTS 

The  central  point  for  the  high  angle  of  attack  design  is  the  trim  flight  condition  at 
Mach  0.3,  altitude  25,000  feet.  The  dynamic  pressine  at  this  point  is  49.5  psf.  is 
the  solution  to  the  linear  quadratic  regulator  problem  with  identity  .staUi  weighting 
and  a  control  weighting  of  Details  on  this  design  are  presented  in 

Appendix  4. 

Fig.  4.5  shows  the  relative  error  results  for  four  test  points.  'Phe  test  points  shown 
in  Table  4.2  represent  a  wide  range  of  high-cx  operating  conditi{)ns. 
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Table  4.2  Test  Points  for  High-a  Linear  Analysis 


Frequency  (rad/s) 

Fig.  4.5  I-telative  Frror  K.esuU.s  High-u 


The  relative  error  is  leas  than  db  for  frequencies  less  than  10  rad/s.  As  in  the  low-a 
results,  this  is  primarily  the  result  of  interactions  between  the  equalization  loop  and 
actuator  dynamics.  I’he  errors  shown  in  Fig.  4.5  are  considered  small  for  frequencies 
of  interest. 
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4.5  OUTER  LOOP  DESIGN 


This  section  describes  the  formulation  and  results  of  the  outer  loop  control  design 
problem.  Two  separate  outer  loop  designs  are  peiToimed  to  reflect  the  din'erent  Hying 
qualities  requirements  for  the  low  and  high  angle  of  attack  flight  regimes.  The 
formulation  of  the  two  designs  ai'o  the  same,  only  the  values  of  the  design  parameters 
differ.  Section  4.0  describes  how  the  low-cx  and  high-(x  designs  are  blended  to  create  a 
continuous  compensator. 

4.5.1  OUTER  LOOP  SYNTHESIS 

The  outer  loop  compensator  is  designed  to  achieve  two  primary  goals,  Hying 
qualities  and  robustness.  These  objectives  can  be  achieved  directly  by  using  a  special 
formulation  of  p-synthesis.  Fig.  4.6  ahow.s  the  design  model  for  p-synthesis. 


Fig.4.6  Design  Model  for  p'.synthe.sis 


The  dynamic  compensator  that  satisfies  these  goals  is  found  using  a  DK-iteration 
based  around  the  regulated  central  plant,  'fhe  design  model  can  be  redrawn  to 
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illuvslrale  the  n-Hynthesis  p»’oblcin.  Kig.  4.7  shows  an  alternato  reprcsentalioji  of  the 
design  model  which  parallels  the  general  form  shown  in  Fig.  2.8. 


Fig.  4.7  M-A  Form  of  Design  Model 


15y  using  the  above  form,  the  VKSTA  F-l(i  {.l-Sytilbesis  problem  can  be  described  as 
the  attempt  to  solve 

min  inf  sup  d(DM(Kj^i.,^y„ )!)■').  (4.17) 

F|,i..syM  F'- ^  (I) 


4.5.1. 1  IDEAL  MODEL  GENERATION 


As  described  in  section  2.2,  Hying  qualities  are  the  primary  measures  of 
performance  for  a  manual  flight  control  system.  For  the  VLSTA  F-16  design,  a  pilot 

lateral  stick  input  translates  into  a  stability  axis  roll  rate  command,  Pcojn,  una  a  pilot 
pedal  input  commands  a  sideslip  angle,  |icoiii.  A  low  order  equivalent  system  fit  of  the 
complementary  sensitivity  function  will  thend'ore  drive  the  Hying  qualities  results  for 


this  design.  An  ideal  model  of  tiic  desired  ajrcraft  response  to  pilot  inputs  can  be 
formulated  from  the  ideal  low  order  equivalent  system  transfer  function  parameters. 
By  forcing  the  complementary  sensitivity  function  to  take  the  frequency  response 
shape  of  this  ideal  model ,  flying  qualities  can  be  included  in  the  design  process  |4.2]. 
The  ideal  model  for  this  design  has  the  form: 
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where  co^  represents  the  desired  Dutch  roll  frequency,  represents  the  desired 
Dutch  roll  damping,  and  Tg  represents  the  desired  roll  mode  time  constant.  Two 
different  sets  of  these  parameters  are  used  to  reflect  the  dc.r;ired  dynamics  at  low  and 
high  angle  of  attack  conditions. 


4.5. 1.2  PERFORMANCE  WEIGHTING 

Hoo  optimization  in  p-synthe.sis  provides  a  direct  way  of  minimizing  the  error 

between  the  ideal  model  and  the  complementary  sensitivity  function.  A  weight, 
Wp,  forces  this  error  to  be  small  at  frequencies  Ic.ss  tliun  10  rad/s.  The  performance 
weight  is  .synthesized  such  that  the  error  between  the  ideal  model  and  tiie 
complementaiy  sensitivity  function  is  bounded  by  Wp-^  if  the  de.sign  achieves  a  p  of 
less  than  unity.  High  weightings  at  low  frequencies  I'orce  the  steady  state  tracking 
error  to  be  small.  The  error  hound  must  he  tight  at  frequencies  between  1  and  10 
rad/s,  the  region  that  dominates  the  transient  response  of  the  closed  loop  sy.stem  to 
pilot  commands.  If  the  performance  error  is  n()t  reduced  adequately  in  this  frequency 
region,  higher  order  dynamics  v/ill  show  up  in  the  transient  response,  destroying  flying 
qualitie,s. 

4.5. 1.3  ACTUATOR  WEIGHTING 


Actuator  dynamics  can  play  a  significant  role  in  closed  loop  performance  and 
robustness.  If  actuator  dynamics  are  significantly  faster  than  the  desired  control 
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bandwidth,  they  may  be  left  out  of  the  compensator  synthesis  problem.  This  is  often 
not  the  case,  as  it  is  not  with  the  VIS'I'A  F-16.  In  order  to  include  the  domii.ant 
actuator  effects,  a  third  order  approximation  of  actuator  dynamics  is  included  in  the 
outer  loop  design  model.  Practical  considerations  such  as  rate  saturation,  unmodeled 
dynamics,  and  time  delays  are  included  by  weighting  actuator  position,  rate, 
acceleration,  and  jerk  (acceleration  rate)  through  the  design  weight  Generalized 

actuator  models  and  static  weights  are  used  in  the  synthesis  model  to  minimize  the 
order  of  the  synthesis  model.  The  generalized  actuator  model  is; 


,  (U).7)(S5.())^  . 

~  (s+19.7)(s^+2(().71X{-;r,.0)s+(6ri.0)‘'^^ 


(4.20) 


where 


(4.21) 


4.5.1,4  PAIIAME TEH  UNCERTAINTY  WEIGHTING 

Kquali'/ation  errors  and  uncertainty  in  aerodynamic  stability  derivatives  drive  the 
requirements  I'or  robustness  to  parametric  uncortuinlies.  Kobustaoss  to  parameter 
variations  is  directly  incorporated  into  the  p-synthesis  design  model  through  the 
weights  G-Adcs'  As  described  in  section  4.4,  the  inner  loop  does  not  do  a 

perfect  job  of  equalizing  the  plant.  In  this  dc.sigu,  tlio  iiiternul  dynamics  of  dynamic 
inversion  consist  of  the  A  matrix  parameters  Yj^,  sinix,  and  cosix.  By  treating  these 
parameter.^  as  uncerLaintics  in  the  p-synthesis  design,  the  overall  control  sy.stem  will 
be  more  robust  to  the  relative  error.s  due  to  equalization.  The  weighting  matrices  in 
the  -synthe.sis  model  are  scaled  such  that  the  maximum  level  of  uncertainty  in  the  A 
matrix  is  ropresented  as: 


(4.22) 


A  trade-off  exists  between  nominal  and  robust  performance  in  this  p-synthesis 
design.  Flying  qualities  requirements  dictate  a  (  ;ht  bound  on  the  allowable  error 
between  the  ideal  model  and  the  coniplementary  .sensitivitv  function.  Robust 


performance  and  robust  stability  to  parameter  uncertainties  must  be  simultaneously 
achieved  in  a  design  that  yields  p  less  than  unity.  Robust  performance,  in  this  case, 
means  that  the  tight  bound  on  performance  error  must  be  met  even  in  the  presence 
of  these  parameter  uncertainties.  The  tight  performance  bound  was  formulated  to 
force  nominal  performance,  but  this  error  bound  is  too  stringent  of  a  requirement  in 
the  sense  of  robust  performance.  It  is  impossible  to  reduce  p  to  less  than  unity  when 
a  large  amount  of  parameter  uncertainty  is  introduced  to  the  problem.  The  next 
logical  step  seems  to  be  to  relax  the  performance  weight.  The  problem  is  that  this 
causes  the  flying  qualities  specifications  to  be  violated  at  the  nominal  condition.  No 
design  will  be  accepted  if  it  cannot  meet  specs  at  the  nominal  condition.  The 
compromise  solution  is  to  dramatically  reduce  the  ernount  of  parameter  uncertainty 
introduced  in  the  design  model  and  perform  a  comprehensive  post-design  analysis  to 
ensure  that  adequate  stability  and  performance  robustness  ba.s  been  achieved. 
Section  4.7  presents  the  results  of  this  aiialysis. 

4.5.2  LOW  ANGLE  OF  ATTACK  OUTER  LOOP  RESULTS 


The  ideal  model  parameters  for  the  low-a  design  arc  0)j)  =  3.0  rad/s,  =  O.Tl,  and 
Tg  =  0.33  seconds.  The  performance  weight  is: 


(s  +  30) 
(s  +  0.03) 


1*2x2 


(4.23) 


The  generalized  actuators  are  weighted  equally.  'I’he  generalized  actuator  position 
weight  is  l.GxlO’’^  the  rate  weight  is  0.0x10  '^  ,  the  acceleration  weight  is  6.0x10'^, 
and  the  acceleration  rate  weight  is  6.()xl()-7.  The  parameter  uncertainty  weights 
used  in  this  design  arc: 


^Adea 


'0.1  0.1  ■ 
0  0 

0  0 


and  CAd.;.s  = 


0.05  0  0  ■ 

0  0,20  U 


(4.24) 


Fig.  4.8  shows  the  final  structured  singular  value  upper  and  lower  bound.s  that 
result  from  the  DK-iteration. 


iTcqUL'iicy  (rad/s) 

4.8  Low-a  |.l  Bounds  for  f)ulx'r  Loop  Design 

The  resulting  controller  is  stable  and  22nd  order.  I'he  controller  transfer  function  is 
shown  in  polc/zero  form  in  Appendix  4.  Using  balanced  truncation,  the  outer  loop 
controller  is  reduced  to  15th  order.  This  controller  order  reduction  is  an  iterative 
process.  First,  a  balanced  realization  is  found  for  the  full  order  compensator.  The 
balanced  grammians  are  then  examined  to  determine  if  there  is  any  natural  drop-oil’ 
in  the  size  of  the  Hankel  singular  values.  An  initial  tmneation  is  made  and  analysis  is 
performed  to  check  if  the  structured  singular  values  of  the  closed  loop  system  have 
been  significantly  perturbed,  A  new  order  is  chosen  based  on  this  analysis  and 
iterations  are  continued  until  the  minimal  controller  order  is  found  that  does  not 
deteriorate  the  p  bounds.  It  should  be  noted  that  a  great  deal  of  effort  was  not  put 
into  controller  order  reduction,  and  other  approaches  which  could  yield  a  lower  order 
compensator  most  likely  exist.  The  most  promising  direction  is  probably  in  variations 
of  frequency  weighted  balanced  truncations.  Fig.  4.9  shows  the  |i  upper  and  lower 
bounds  for  the  system  with  full  and  reduced  order  controllers. 


(if) 
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Fig.  4.9  Low-a  |i  Bound.s  for  Full  and  Reduced  Order  Controllers 


The  reduced  order  controller  transfer  function  is  shown  in  Appendix  4.  'I'hc  reduced 
order  result  matches  the  full  order  result  very  closely,  therefore  little  or  no 
degradation  ol' robustness  or  performance  should  be  expected. 

Fig.  4.10  shows  the  tran-sfer  function  between  \)  and  Idea)  model 

(dashed  line)  and  for  the  closed  loop  sy.stem  (solid  line)  at  the  design  condition.  The 
magnitude  responses  match  very  clo.sely  for  frequencies  less  than  10  rad/s. 


ITcqucncy  (rad/scc) 

Fig  4.10  how-fx  lV(icoii,  for  Ideal  Model  and  Outer  Loop  Design 
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Fig.  4.1 1  shows  the  transfer  function  between  n  and  Ucorn  ideal  model  (dashed 

line)  and  for  the  '  sed  loop  system  < solid  line)  at  the  design  condition.  The 
compensated  system  matches  the  desired  first  order  response  very  closely  for 
frequencies  less  than  10  rad/s. 


♦  ♦ 

Fig.  4.11  Low-(i:  (.i/Moom  Ideal  Model  and  OuU  r  Loop  Design 


Figures  h.V2  and  4.  II'  show  linear  responses  fur  unit  step  commands  in  sideslip  and 
otanility  axis  roll  ra"  L'did  lines  represent  the  closed  loop  .system  response  including 
full  order  actuator  models  ■.vhiie  the  dashed  line  represent  the  response  of  the  ideal 
ruodel.  Results  are  shown  for  the  do.'  ign  condition  as  well  as  for  test  conditions  1-4 
h  orn  'Fable  4.1 . 
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lime  (sec)  lime  (sec) 

fig.  4.12  Low-a  Pcom  H^i.spon.se  at  Design  and  Test  Points 
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fig.  4.13  Low  it  Ke.spoiise  at  Design  and  'Pest  Points 

Tht  time  response  plots  show  that  both  tracking  and  respon.se  decoupling  arc 
maintained  across  the  low-oi  design  (light  envelope,  demonstrating  the  success  oCthe 
equalization  in  preserviwjj  nominal  pcrfhnnanco. 


(is 


4.5.3  HIGH  ANGLE  OF  ATl’ACK  OUTER  LOOP  RESULTS 


The  ideal  model  parameters  for  the  hi{,'h-(x  design  arc  =  1.0  rad/s,  =  1.0,  and 
Tr  =  1.2  seconds.  The  performance  weight  is: 


_  0.75(s  +  30)  . 

"  (s  +  0.03)'  • 


(4.25) 


The  generalized  actuators  are  weighted  equally.  The  generalized  actuator  position 
weight  is  1.5x10  the  rate  weight  is  O.O.’vlO  ,  the  acceleration  weight  is  G.OxlO  -'’, 
and  the  acceleration  rate  weight  is  0.0x10  The  parameter  uncertainty  weights 
u.sed  in  ♦’  »s  design  are; 


0.1 

0.1  ^ 

0.05 

t) 

0  ■ 

0 

0 

and 

0 

0.26 

0  _ 

0 

0 

(4.26) 


Fig.  4.14  shows  the  final  structured  singular  value  upper  and  lower  hounds  that 
result  from  the  l)K-iteration. 


Fig.  4. 14  High-a  ^  Hounds  for  Outer  hoop  Design 


'fhe  resulting  controller  is  stable  and  22nd  order.  Using  balanced  truncation,  the 
outer  loop  controller  is  reduced  to  13th  order.  '^I'he  full  and  reduced  order  controller 


transfer  functions  aix*  shown  in  Ajjpendix  4.  Pip.  4.  Ii5  shows  the  fl  upper  and  lower 
bounds  for  the  system  with  full  ajid  reduced  (.rder  controllers. 


Pig.  4. 15  Kigh-a  (I  Bounds  for  h'ull  and  lieduced  Order  Controllers 


The  reduced  order  result  matches  the  full  order  result  very  closely,  therefore  little  or 
no  degradation  of  robustness  or  performance  should  be  expected. 

Fig.  4.16  shows  the  transfer  function  between  p  and  for  the  ideal  model 
(dashed  line)  and  for  the  closed  loop  system  (solid  line)  at  the  design  condition.  The 
magnitude  responses  match  very  closely  for  frequencies  less  than  10  rad/s. 


I•rcqucllcy  (rail/scc) 


Fig.  4. 16  High-rt  [Vllcom  Ideal  Model  and  Outer  Loop  Design 

Fig.  4.17  shows  the  transfer  function  between  p  and  for  the  ideal  model  (dashed 
line)  and  for  the  closed  loop  system  (solid  line)  at  the  design  condition.  The 
compensated  system  matches  the  desired  first  order  response  very  closely  for 
frequencies  less  then  10  rad/s. 


ticqucDv'y  nad/.sc(. ) 


Fig.  4.17  Hi{^h-(X  |.i/|.icom  ’deal  Model  and  Outer  L-oop  Design 
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Fig.  4.18  and  4.19  show  linear  responses  lor  unit  step  commands  in  sideslip  and 
stability  axis  roll  rate.  Solid  lines  represent  the  closed  loop  system  response  including 
full  order  actuator  models  while  the  dashed  line  represent  the  response  of  the  ideal 
model.  Results  are  shown  for  the  design  condition  as  well  as  for  test  conditions  1-4 
from  Table  4.2. 


Fig.  4.18  High-a  pcom  liesponso  at  Design  and  Test  Points 
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Fig.  4. 19  High-(x  Response  at 


I)(‘sign  and  'I’est  Points 
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The  time  response  plots  show  that  both  tracking  and  response  decoupling  are 
maintained  across  the  high-a  design  flight  envelope,  demonstrating  the  success  of  the 
equalization  in  preserving  the  nominal  performance  of  the  central  design  point. 

4.6  CONTROLLER  IMPLEMENTATION 

The  low-a  and  high-ot  controllers  must  be  combined  to  form  a  complete  full- 
envelope  control  system.  A  blending  parameter,  C,  can  be  introduced  to  create  a 
transition  region  in  angle  of  attack. 


C  -  1 

_  17.5  -  (X 

17'.5  -  12.5 

C  =  0 


if  (x<  12.5 
if  12.5  <a<  17.5 
if  (X  >  17.5  . 


(4.27) 


The  combined  low/high-(x  desired  dynamics  and  outer  loop  controllers  can  now  be 
represented  as  a  function  of  the  blending  parameter. 


~  ^l(ivv-a^ '  '^high-a^  ^  ^ 


1) 


(4.28) 

(4.29) 


The  implementation  of  the  outer  loop  controller  is  shown  in  fig.  4.20.  P  represents 
the  equalized  system  including  the  control  selector  and  equalization  feedbacks. 


The  outer  loop  compensator  accepts  inputs  of  sideslip  and  stability  axis  roll  rate  error 
and  generates  a  sideslip  acceleration  and  stability  axis  roll  acceleration  command 
signal. 

Fig.  4.21  shows  the  overall  structure  of  the  VISTA  F-16  inner/outer  loop  controller. 


CONTROLLER 

Fig.  4.21  Inner/Outer  Loop  Controller 


The  combined  controller  receives  inputs  of  sideslip,  body  axis  roll  rate,  body  axis  yaw 
rate,  sideslip  command,  and  stability  axis  roll  rate  command.  The  controller 
generates  actuator  commands  for  asymmetric  flap,  asymmetric  horizontal  tail,  and 
rudder  dellections.  As  described  earlier,  the  equalization  gains,  and  the  control 
selector,  CS,  are  functions  of  Mach,  altitude,  and  angle  ol’  attack.  The  outer  loop 
controller,  is  a  dynamic  compensator  that  is  a  linear  function  oi'angle  of  attack. 

4.7  ROBUSTNESS  ANALYSIS 

A  robustness  analysis  model,  shown  in  Fig.  4.22,  is  built  to  rellect  the  robustness 
requirements  for  structured  uncertainty,  unstructured  uncertainty,  and  performance. 
The  resulting  M-A  form  block  diagram  has  inputs  and  outputs  that  reflect  a  fictitious 
performance  block  (wj  -  zj),  a  structured  uncertainty  block  (  -  zo),  .vn  unstructuied 

uncertainty  block  at  the  input  to  the  actuators  (w;^  -  Zj),  aiid  an  um  h  ucturod 


uncertainty  block  at  the  output  of  the  sensors  (W4  -  z^).  Because  of  the  lack  of  highly 
accurate  models  for  the  different  uncertainties,  analysis  is  performed  considering 
each  type  separately.  This  approach  allows  problem  areas  and  trends  in  robustness 
to  be  isolated  and  identified. 


Fig.  A.2'2.  Robu.stness  Analysis  Model 

4.7.1  LOW  ANCLE  OF  ATTACK  ROBUSTNESS  RESULTS 

4.7. 1.1  LOW  ANGLE  OF  ATTACK  STRUCTURED  UNCERTAIN  IT 

Structured  uncertainty  requirements  for  the  ViS'l'A  r  - 16  lateral/directional  design 
are  driven  by  perturbations  in  aerodynamic  parameters.  Seven  stability  derivatives 
and  seven  control  dcrivative.s  are  identified  for  robustness  analysis.  We  can  write  the 
perturbed  state  equations  as: 


X  =  (A  +  AA)x  +  (B -f- AB)u  , 


(4.30) 


where 
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AA  = 


0  0  AYsr 

AB  =  ALsht  ALsdf  ALsr 
..  ANsdT  ANfinir  ANfiR 


(4.31) 


AYp  0  0 

AL)j  ALp  AL, 

_  AN[j  ANp  AN,-  _ 


J 


The  structure  of  the  uncertainty  can  be  captured  by  rewriting  the  perturbed  state 
equations  as; 


X  =  Ax  +  Bu  +  B^A^pZ2 


(4.32) 


(4.33) 


so  whencKA^  B)  <  1,  the  maximum  parameter  uncertainty  in  the  .system  matrices  is 
represented  as: 


AA  =  nA(:A  AB  =  BPa 


(4.34) 


The  level  of  parameter  uncertainty  is  captured  by  the  matrices: 
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(4.35) 
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(4.36) 
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ALsdt  0  0 

ANsdT  0  0 

0  ALaop  0 
0  ANsdp  0 
0  0  AYfiR 

0  0  ALsr 

0  0  ANgR  - 


(4.37) 


The  level  of  uncertainty  for  a  parameter  at  each  operating  condition  is  captured  as  a 
percentage  of  its  nominal  value.  The  uncertainty  levels  used  for  stability  robustness 
analysis  are  shown  in  Table  4.3.  The  values  are  based  upon  the  accuracy  of  the 
VISTA  F-16  aerodynamic  database  and  the  uncertainty  database  describedi  in  [4.3]. 


Table  4.3  Structured  Uncertainty  Levels 


stability  derivatives 

control  derivatives 

AYp  =  0.15  Y(i 

AYsr  =  0.15  Ysr 

AL(3  =  0. 10  L(} 

ALsdt  '=  O’ 15  Lsdt 

ALp  =  0.30  Lp 

ALsdf  =  0.10  Lsdf 

AL,-  =  0,20  L, 

ALsr  =  0.40  Lsr 

AN(5  =  0.30  N|i 

ANgDT  =  0.15  NftDT 

ANp  =  0.50  Np 

AN^df  =  0.20  Nsdf 

AN,.  =  0.15  N,. 

ANsr  =  0.15  Nsr 

Stability  robustness  to  structured  uncertainties  is  tested  for  the  cmtral  plant  and 
the  four  test  points  described  in  Table  4.1.  Fig.  4.23  shows  the  resulting  upper  and 
lower  bounds  on  the  structured  singular  value. 
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i-Vcqucncy  (rad/s) 

Fig.  4.23  Low-(x  |i  Bounds  for  Robust  Stability  to  Structured  Uncertainty 

For  all  conditions  tested,  jd  is  less  than  one  for  all  frequencies.  The  desired  level  of 
stability  robustness  to  structured  uncertainty  has  been  achieved. 

4.7.1.2  LOW  AT^GLE  OF  ATT  ACK  UNSTRUCTURED  UNCERTAIN  IT 

Unstructured  uncertainty  at  the  actuator  input  represents  both  unmodeled 
dynamics  and  saturation  effects.  Because  structured  .singular  value  stability 
analysis  methods  have  only  recently  been  wddely  accepted,  little  data  is  available  on 
what  appropriate  levels  of  uncertainty  should  be,  'fhe  reijuirements  for  robustness 
for  unstructured  uncertainty  for  this  design  are  derived  from  the  classical  gain  margin 
specification  of  6  db  in  each  loop.  This  roughly  corresponds  to  a  ^)()'?^  nuiltiplica''ve 
uncertainty  stability  tolerance  in  individual  loops.  This  requirement  is  extended  to  the 
case  of  simultaneous  multiplicative  perturbations  by  relaxing  the  level  of  uncertainty 
to  30%.  Fig.  4.24  shows  .structured  singular  value  analysis  results  for  the  design 
condition  and  test  points  1-4  u.sing 


Wact  - 


0.3  0  0  ' 

0  0.3  0 

0  0  0.3 


(4;i8) 


7H 


Fig,  4.24  Low-a  jj  Bounds  for  Unstructured  Uncertainty  at  Plant  Input 


The  desired  level  of  robustness  has  been  achieved  to  unstructured  uncertainty  at  the 
plant  input  for  all  conditions  tested.  Analysis  also  reveals  a  trend,  the  lower  dynamic 
pressure  test  cases  exhibit  lower  levels  of  robustness  for  this  type  of  uncertainty. 
The  steady  degradation  of  robustness  is  a  result  of  higher  feedback  gains  that  arc 
required  for  plant  equalization  at  low  dynamic  pressure  conditions. 


Levels  of  unstructured  uncertainty  at  the  sensor  output  are  driven  by  the  quality 
of  infonnation  available  in  each  sensor  channel.  The  body  axis  rotational  rates  can 
be  measured  by  gyros  which  are  preci.se  and  reliable.  A  .sideslip  signal,  on  the  other 
hand,  normally  must  be  estimated  or  recon.structed  usi  g  complementary  filtering. 


TT 


Uncertainty  in  the  sideslip  measurement  is  the  driving  requirement  for  robustness 
analysis  at  the  plant  output.  Structured  singular  value  analysis  is  performed  for  a 
10%  multiplicative  uncertainty  in  body  axis  rotational  rate  signals  and  a  40% 
multiplicative  uncertainty  in  sideslip  measurement. 


W 


sens 


0.4  0  0 

0  0.1  0 
0  0  0.1 


The  results  for  the  design  condition  and  test  points  1-4  are  sh($wn  in  Fig.  4.25. 


(4.39) 
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Fig.  4.2f>  Low-a  ji  Bounds  for  Unslructured  IJncortuinty  at  Plant  Output 

The  desired  level  of  robustness  has  been  achieved  to  unstructured  uncertainty  at  the 
sensor  outputs  for  all  conditions  tested.  Fig.  4.25  shows  that  the  same  reduction  of 
robustness  margins  at  low  uynamic  pressures  is  present  in  the  output  uncertainty 
analysis  as  it  is  with  input  uncertainty  analysis. 

4.7. 1.3  LOW  ANGLE  OF  AITACK  ROBUST  PERFORMANCE 


Performance  robustness  is  defuied  for  this  problem  as  achieving  a  small  relative 
error  between  the  ideal  nwdel  and  the  actual  response  in  the  presence  of 
uncertainties  in  aerodynamic  parameters.  The  levcd  of  structured  uncertainty 
included  in  this  analysis  i.s  reduced  from  that  used  for  robust  stability  analysis.  The 
stability  robu.stness  analysis  uncertainties  were  chosen  to  he  worst  case  te.sts  ul 
closed  loop  stability.  It  is  not  rea.sonanie  to  expect  performance  rotnistuess  (.o  the 
same  worst  case  set.  Performance  robustness  analysis  is  performed  u.sui^; 
structured  uncertainty  levels  of  25‘/f ,  the  si/.e  of  those  shown  in  Table  4  2  'I'hc 


performance  weight,  Wp,.rf,  nsed  for  this  analysis  is  not  the  same  as  the  weight  used 
lor  control  synthesi.s.  'I’hat  design  weight  was  chosen  witli  otlier  consuiecatidus  in 
mind,  namely  good  nominal  performance.  The  inverse  of  the  pei  forniar.ce  weight, 
shown  in  Fig.  4.26,  bounds  the  allowable  error  between  the  analysis  medei  and  the 
idea!  nuHlef  in  the  presence  of.structureil  uncertainties 


KO 


STcquL'iicy  (rad/soc) 

Fig.  4.26  Performance  Error  Bound  for  Robustness  Analysis 

The  performance  weight  dictates  that  the  steady  state  error  to  commands  must  be 
less  than  10%.  Fig.  4.27  shows  the  results  of  performance  robustness  analysis  for 
the  design  point  and  test  points  1-4. 


I'rctiucncy  (raU.^'i) 


Fig.  4.27  Ia)w-(X  p  Bounds  for  Performance  Robustness 
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The  structured  singular  values  are  less  than  unity  for  all  frequencies,  so  performance 
robustness  has  been  achieved  for  all  conditions  tested. 

4.7.2  HIGH  .\NGLE  OF  ATTACK  ROBUSTNESS  T  "JLTS 


4.7.2. 1  HIGH  ANGLE  OF  ATTACK  STRUCTURED  UNCERTAINTY 

Stability  robustness  to  structured  uncertainties  is  tested  for  the  high-a  central 
plant  and  the  four  test  points  described  in  Table  4.2.  The  levels  of  uncertainty  in 
Table  4.3  are  used  for  the  analysis.  Fig.  4.28  shows  the  resulting  upper  and  lower 
bounds  on  the  structured  singular  value. 


1-icqucncy 

(rad/.scc) 

P‘':g.  4.28  High-tt  (.1  Bounds  for  Robust  Stability  to  Structured  Uncertainty 


For  all  conditions  tested  p  is  les.s  than  one  for  all  frequencies.  'I'he  desired  level  of 
stability  robustness  to  structured  uncertainty  has  been  achieved. 


4.7.2.2  HIGH  ANGLE  OF  ATTACK  UNSTRUC  nJRED  UNCERTAINTT 

Robustness  to  a  30%  per*  irbation  at  the  actuator  input  is  tested.  Fig.  4.29  shows 
structured  singular  value  analysis  results  for  tlie  central  plant  and  high-(x  test  points 
1-4. 
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Frequency  (rml/scc) 

Fig.  4.29  High-a  [.1  Bounds  I’or  Unstructured  Uncertainty  at  Plant  Input 
The  desired  level  of  robustness  to  uncertainty  at  the  plant  input  has  been  acliieved. 

Fig.  4.30  shows  the  results  of  robustness  analysis  for  a  10%  uncertainty  in  roll  and 
yaw  rate  measurements  and  a  40%  unceiiainty  in  sideslip  measurement. 


Fig.  4.30  High-(x  }i  Bounds  for  Unstructured  Uncertainty  at  Plant  Output 


For  all  conditions  tested,  the  desired  level  of  robustness  to  output  uncertainties  has 
been  easily  achieved  for  the  central  plant  and  the  four  test  cases. 

4.7.2.3  HIGH  ANGLE  OF  ATIACK  ROBUST  PERFORMANCE 

Robust  performance  analysis  is  performed  for  the  Vn^h-a.  central  and  test  points 
using  the  same  weights  and  levels  of  structured  uncertainty  described  in  section 
4. 7. 1.3.  Fig.  4.31  shows  the  results  of  this  analysis. 


ficqucncy  (radAscc) 

Fig.  4.31  High-a  |i  Bounds  for  Performance  Robustness 

The  flight  control  system  maintains  good  steady  state  trucking  and  satisfactory 
flying  qualities  in  the  face  of  uncertainties  in  aerodynamic  and  control  derivatives. 

4.7.3  ROBUSTNESS  ANALYSIS  OF  LOW/HIGH-(x  BLENDING 

The  controller  blending  introduced  in  section  4.6  creates  a  potential  stability 
problem.  While  it  has  been  proven  that  the  low'-cx  and  high-a  flight  conditions  are 
stabilized  by  their  respective  controllers,  nothing  has  been  demonstrated  about  the 
stability  properties  of  the  blending.  Eqs.  (4.28)  and  (4.29)  can  be  rewritten  such  that 
the  stability  of  the  control  blending  becomes  an  uncertainty  problem  which  can  be 
solved  with  structured  singular  value  methods.  liedefine  the  blending  paramek'f  a.s 
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C  =  0.5  (I  +  A) ,  where  A  =  5l,  and  ‘o(A)  <  1  (4.40) 

The  desired  dynamics  and  outer  loop  controller  can  now  be  represented  as  a  function 
of  the  uncertain  diagonal  repeated  block  A. 

V=  0.5(V}^igh.(^  +  +  0.5A(V[^jgj^.(y  -  (4.41) 

=  0-5(K^ihigh.<x'*'  ^iiow-«^  ^i-low-a^  i^A2) 


Fig.  4.32  shows  the  robustness  analysis  model  for  control  blending. 


Fig.  4.32  IiOW/'High-{X  Blending  Robustness  Model 

If  the  .system  can  tolerate  a  with  an  infinity  norm  less  than  unity,  the  control 
blending  is  guaranteed  not  to  de.stabilize  the  closed  loop  system  for  any  value  of  C 
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between  0  and  1.  Fig.  4.33  shows  the  results  of  structured  singular  value  analysis  for 
the  ten  flight  conditions  described  in  Tables  4.1  and  4.2. 


The  control  blending  is  guaranteed  not  destabilize  the  system  for  any  point  tested. 
Notice  that  peaking  occurs  in  the  stuctured  singular  value  plot  for  conditions  where 
angle  of  attack  is  less  than  12.6  degrees  and  greater  than  20  degrees.  In  the  region 
including  angle  of  attack  from  12.6  to  20  degrees,  the  structured  singular  values  are 
very  well  behaved.  The  blending  region  described  in  section  4.6  should  provide  a 
smooth  transition  between  the  low-  and  high-a  controllers. 
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4.8  NONLINEAR  ANAI.YSIS 


In  order  to  f)e  jiidf^ed  acceptable,  a  manual  flight  control  system  must  not  only 
meet  the  linear  design  criteria,  it  must  exhibit  satisfactory  properties  in  a  dynaniic 
nonlinear  simulation  environment.  A  sot  of  nonlinear  evaluation  manuevers  is 
created  to  test  the  ability  of  the  control  system  to  handle  large  amplitude  and  rapidly 
changing  inputs.  These  manuevers  arc  performed  at  a  range  of  flight  conditions 
including  both  low  and  high  angles  of  attack.  Tracking  performance  and  response 
decoupling  are  the  primary  measures  of  ir.erit  during  these  manneveis. 

4.8.1  CONVENTIONAI.  MANEUVERS 

(1)  Loaded  Roll 

'fhe  aircraft  is  ba.nked  tit)  deg  to  tlie  right  and  loadc'd  to  'A-6  g's.  While 
maintaining  the  g  loading,  a  roll  is  executed  at  6t)-120  deg/s  to  bring  the 
air’Taft  to  60  deg  loll  bank,  'fhe  aircraft  is  then  unloaded  and  the  vehicle  is 
rolled  back  to  wings  level.  Fig.  4.84  shows  a  loaded  roll  maneuver.  Fig.  4.88 
shows  the  time  histories  for  a  loaded  roll  performed  from  a  starting  condition  of 
Mach  0.7  and  altitude  20,000  feet,  an  off  design  point.  The  time  re.sponse  plots 
show  that  tracking  performance  is  maintained  for  large  amplitude  commands. 
Sideslip  is  regulated  lo  less  that;  1  degree  during  the  higliJy  coupled  maneuver. 

(2)  Higli  Kate  Roll 

A  large  amplitude  stability  axis  roll  rate  input  is  applii.'d  to  roll  the  aircrall 
through  860  deg.  fig.  4.85  shows  a  high  ixite  roll  performed  a!  Mach  0.6  and  all.itude 
20,000  feet.  'I'he  time  response  plots  in  Fig.  4.80  sliow  that  for  a  250  deg/s  stability 
axis  roll  rate  command,  g;)od  tracking  and  response  decoupling  are  maintained.  The 
maximum  .sideslip  is  loss  than  2.5  deg. 

(8t  Steady  Fidesl ip 

A  1  deg/s  .sideslip  command  is  applied  until  a  sideslip  of  10  deg  is  achieved,  'i'he 
command  is  then  reversed  to  bring  the  aircraft  to  10  deg  sideslip  in  the  opposite 
direction.  Fig.  4.86  .shows  a  steady  sideslip  manuevor  performed  at  Mach  0.6  and 


altitude  20000  feet.  The  time  responses  in  Kig.  4.40  .sh(iw  that  the  sideslip  command 
is  accurately  tracked  while  a  wings  level  attitude  is  maintained.  Notice  that  at  10  deg 
sideslip,  the  rudder  is  close  to  position  saturation. 

(4)  Turn  Reversals 

Large  amplitude  roll  commands  are  apphed  to  rapidly  bring  the  aircral't  to  60 
deg  left  roll  angle,  then  to  60  deg  right  roll  angle,  back  to  60  deg  left,  and  then  to  wings 
level.  Fig.  4.37  shows  a  turn  reversal  manuever  performed  at  Mach  0.6  and  altitude 
20,000  feet.  The  time  responses  in  Fig.  4.41  show  that  the  rapidly  changing  large 
amplitude  stability  axis  roll  rate  commands  test  the  ability  of  the  control  .system  to 
provide  tracking  and  turn  coordination  while  avoiding  actuator  .saturations,  'I’lie  peak 
sideslip  during  this  manuever  is  less  than  3.0  deg. 
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4.8.2  HIGH  ANGLE  OF  AITACK  MANEUVERS 


(1)  High-cx  High  Hate  Roll 

A  large  amplitude  stability  axis  roil  rate  input  is  applied  to  roll  the  aircraft 
through  360  deg.  Fig.  4.42  shows  a  high  rate  roll  performed  at  Mach  0.3,  altitude; 
25,000  feet,  and  23  deg  angle  of  attack.  The  arrow  in  this  figure  indicates  the 
aircraft’s  velocity  vector.  The  time  respon.se  plots  in  Fig.  4.44  show  that  for  a  45 
deg/s  stability  axis  roll  rate  command,  good  tracking  and  response  decoupling  are 
maintained.  The  maximum  sideslip  is  less  tiian  4.1  deg.  Notice  tliat  the  rudder  rate 
and  position  saturates  at  the  beginning  of  this  maneuver  and  asymmetric  Haps  also 
rate  saturate.  These  limits  are  encountered  because  of  the  physical  limitations  of  the 
vehicle  at  elevated  angles  of  attack.  It  is  important  to  note  that  the  control  system 
does  not  drive  Oio  system  unstalih'  when  the.se  saturation:;  occur. 

(2)  High-(x  Steady  Sideslip 

A  1  deg/s  sideslip  command  is  applied  until  a  sideslip  of  10  deg  is  achieved.  The 
command  is  then  reversed  to  bring  the  aircraft  back  to  U  dog  sideslip.  Fig.  4.43  shows 
the  aircraft  during  a  steady  sideslip  maneuver  performed  at  Mach  0.3,  altitude 
20,000  feet,  and  23  deg  angle  of  attack.  The  time  responses  in  Fig.  4.45  show  that 
the  sideslip  command  is  accurately  tracked  while  a  wings  level  attitude  is  maintained. 
Notice  that  at  10  deg  sideslip,  the  maximum  control  delloction  takes  place  in 
asymmetric  flaps  rather  than  rudder.  3'his  is  because  at  high  angle.s  (ji'  attack,  a 
sideslip  or  stability  axis  yaw  requires  a  body  axis  roll. 


4.43  Hi<^h-aiSteady  Sideslip  Maneuver 
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4.9  CONCLUSIONS  AND  LESSONS  i/lIAilNED 


This  example  deinenstraLes  an  application  of  dynamic  inversion  and  structured 
singular  value  synthesis  to  a  manual  ilight  control  problem  for  the  VhSTA  I<'-](').  'I'he 
design  is  perfonned  using  the  inner/outer  loop  control  metiiodulogy  described  in  section 
3.  Detailed  linear  and  nonlinear  analysis  shows  that  the  resulting  control  laws 
provide  the  desired  performance  and  robustness  across  a  wide  flight  envelope. 

A  number  of  lessons  learned  can  he  derived  from  the  experience  of  designing 
lateral  control  laws  for  the  VlSl'A  F-13.  Dynamic  inversion  is  a  useful  approach 
which  eliminatOvS  the  need  fur  gain  scheduling,  hut  it  may  yield  poor  results  if  applied 
naively.  The  selection  of  the  desired  dynamics  has  a  significant  efi'ect  on  the  overall 
robustness  of  a  design.  The  relationship  betwiien  the  desired  dynamics  and 
robustness  is  a  topic  worthy  ol‘ further  ivsearch.  Another  lesson  learned  is  that  it  is 
difficult  to  include  a  high  level  of  robustness  in  the  p-syntliesis  design  model  and  still 
achieve  a  high  performance  control  law.  'I'liis  is  duo  primarily  to  two  factors:  physical 
limitations  due  to  trade-offs  between  performance  and  robustness,  and  the 
conservatism  that  comes  with  representing  real  perturbations  as  complex  in  the 
design  model.  Fortunately,  control  laws  designed  using  a  reduced  uncertainty  set 
provide  satisfactory  levels  of  robust  ness. 
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APPENDIX  4 


All  of  the  linear  models  included  in  the  appendix  are  pivcn  in  the  form; 


Yii  sirux  -cos(x 


N|i  Np  N, 


A  =  !  L(5  hp  L;-  ,  H  =  hfiirr  hsDF  hsK 


Yxdt  Ysdf  Ysr 


NfiDT  NsdF  Nfilj 


(A4.1) 


LOW  ANGLE  OF  ATl'ACK  DESIGN  AND  ANALYSIS  IVLATRICES 

The  state  space  representation  of  the  low-a  central  model  is: 


■  -1.G8H5O-01  7.r,9'l!)(>-02  .9.!)523o-01 
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■  2.3384(-02  :i.7()19e-Oa  2.r.281o-02 

hec/Rm/  =  •^•'24C)4(‘+01  2.95070+01  ().()951e+0() 

L  -2.3973C+00  -5.27G4e-01  -2.GG22e+0() 


(A4.3) 


The  LQR  gain  matrix  described  in  section  d.4.2  is; 


r  -8.:{248o-01  l,722Io-01  1.02.Glo  01  “I 
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The  low-a  desired  dynamics  for  the  dynamic  inversion  calculations  are 
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1()W-(X  test  point  1/ 1  is; 
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low-a  test  point  #2  is: 
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low-a  test  point  #2  is; 
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low-a  test  point  #4  is; 
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The  state  space  representation  of  the  hii;h-a’.  mi/.ra/  model  is; 
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The  high-a  desired  dynamics  for  the  dynainic  inversion  calculations  are  ■ 
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high-a  test  point  #1  is; 
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high-a  test  point  #2  is: 
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high-a  test  point  #3  is: 
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-1.7573e-f01  -l.C894e+01  3.7761e■^00 

.  9.9705e-01  -1.8()26e-01  -4.3128o-01. 

. -1.1659O-1-00  2.5008e-01  -2.2423e-i90 

(A4.1G) 


high-a  test  point  #4  is; 


-3.6933C-02 
-1.26790-1-01 
2,084 1(^-^00 


5.0019O-U1  -8.6457e-0] 
-1.7i45o-01  i.i820o-i-00 
-1.3896e-01  -2.9045e-01 


9.05610-03 

-7.978Ge'i-00 

-1.7620O-01 


-3. 11490-03  1.5815e-02  ' 
-5.40C'7e-i-00  2,2032e-i-00 
3.1292e-01  -9.5343e-01 


(A4.17) 
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OUTER  LOOP  CONTROIXER  MATRICES 


The  full  and  reduced  order  outer  loop  controllers  have  the  tbilowing  form: 


Kg 


numl.1  numI2 
num21  num22 

don 


where 


num  =  K(s  +  zi)(s  4  Z2)...(s  4-  4-  z„) 


and 


den  =  is  4-  p])(s  4-  p2)...(s  4-  p„.i)(s  (  p„) 


(A4.18) 


(A4.19) 


(A4.20) 


Table  A4.1  Full  Order  Low-cx  Outer  Loop  Controller 


numll 

numl2 

num21 

num22 

den 

K  (P22) 

2.5721e+03 

-8.3453e+01 

4.6524e+01 

4.0843e+03 

-6.9474e+02 

Z2I  (P2l) 

-6.9474e+02 

-6.9443e+02 

-6.9477e+02 

-6.9444e+02 

-6.94G9e+n2 

Z20  (P2o) 

-6.9469e+02 

-6.9468e+02 

-6.9469e+02 

-6.94690+02 

Zl9  (P19) 

-6.9469e+02 

ZI8  <^Pl8^ 

-9.065  le+01 

mujmii 

-8.8507e+01 

-4.6643e+01+ 

4.5350e+01i 

-9.0650e+01 

Zi7  (P17) 

-4.7124e+01+ 

4.5254e+01i 

-4.9275e+01+ 

4.5883o+01i 

-4.3328e+01 

+4.5444c4-01i 

-4.6643e+01- 

4.5350e+01i 

-8.9030e+01 

Z16  (Pie) 

-4.9275e+01- 

4.5883e+01i 

-4.3328e+01- 

4.5444e+01i 

215  (P15) 

-2.0982e+01+ 

2.5265e+01i 

6.5360e-01+ 

2.6938e+01i 

-1.0606e+00 

+4.4836e+01i 

-1.9293e+01 

-2.097  le+01- 
2.5276e+01i 

214  (P14) 

-2.0982e+01- 

2.5265e+01i 

6.5960e-01- 

2.6938e+01i 

-1.0606e+00- 

4.4836e+01i 

-5.5413e+00+ 

1.0974e-f-01i 

-6.4138e+00+ 

1.0886e+01i 

213  (P13) 

-1.6045e+01 

-1.5847e+01+ 

1.7131e+00i 

-3.2032O+01 

-5.5413e+00- 

1.0974e+01i 

-5.4138e+00- 

1.0886e+01i 

Z12  (P12) 

-1.6385e-}-01 

-1.5847e-*01- 

1.7131e+00i 

-1.8881e+01 

+1.2690e+01i 

-1.2267e+01+ 

2.0950e+00i 

-1.6021e+01 

211  (Pll) 

-5.301  le+0()+ 
4.2606e+00i 

1.6723e+00+ 

4.9997e+00i 

-1.8881es-01- 

1.2690e+01i 

-1.2267e+01- 

2.0950e+00i 

-1.1339e+0l 

ZlO  (PIO) 

-5.301  ie+00- 
4.2606e+00i 

1.6723e+00- 

4.9997e+00i 

-1.6102e+01 

-9.7888e+00 

29  (pg) 

-1.40l3e40b+ 

3.4410e+00i 

-3.8339e+00+ 

1.4163e+00i 

-5.3073e+00 

+4.2694o+00i 

-5.5771e+00- 

1.2845e+00i 

-4.7406e+00 

28  (PS^ 

-1.4013e+00- 
3.441  Oe-tOOi 

-2.6205e+00 

-1.8338e+00 

27  (P?) 

-1.8380e+0b 

-1.8380e+00 

-1.83a0e+00 

-1.8380e+00 

zg  (pg) 

-1.8379e+00 

-1.8292e+00 

-1.8380e+00 

-1.8365ei00 

25  (P5) 

-7.3952e-04 

1.1408e-02 

-7.0672O-04 

-2.3573C-03 

-8.2088e-04- 

1.8614e-04i 

Z4  (P4) 

Q  /xooOi,  no  1 
-0  C/C”V/ 

1.5362e-G7i 

-7.6371e-03 

-7.-'.055e-03 

-7.1433e-03 

-3.52GGe-02 

23  (P3) 

-3.5464e-02 

-3.5228e-02 

-7.1055e-03 

-3. 52290-02 

22  (P2) 

-7.1054e-03 

-7.1055e-03 

-3.5229e-02 

-3.5236e-02 

-7.1055e-03 

21  (pl) 

-7.1055e-03 

-3.5229e  02 

-7.1061e-03 

-3.5229e-02 

-7.1055e-03 

'^^i'ahle  A^'  .2  Order  Hi/.'h-oc  ( )ui.<'r  I -ooii  ( ’nnirollcT 


nunill 

ruuulC 

iuan21 

num22 

(ien 

K  (P2fi) 

1.426604-02 

-9.67356+00 

-3.57416+00 

3.47506+02 

■4.76606+02 

Zl^,5  (P25) 

-4.7660O-1-02 

-4.76560+02 

-4.76606+02 

-4.76556+02 

-1.0012e+02 

7-24  (P24) 

-4.7659e+G2 

-4.76596+02 

-4.7659o+0T^ 

-4.76596+02 

-9.99996+01 

723  (P23) 

-4.76596+02 

-4.76596+02 

-4.7659e+02 

^  -4.76596+02 

-4.76596+^02 

7-22  (P22^ 

-1.00126+02 

-1.00016+02 

-  1.0002O+02 

^■-4.7.3466+01^ 
4.523769  011 

-4./659e+02  + 

1 .66406  !  01 

-4.76906+01  + 
4.5164e+01i 

-5.38246  1-01 1 
4.67026+011 

-5.2048o+()l  + 
4.66236+011 

-4.7346e+()l- 
4.52376  t-O  11 

-4.76596+02- 

1.66406-101 

720  ^P20) 

-4.7699e+01- 

+.5iG4o+01i 

-5.3824e+01- 

4.6702e-,4)li 

-5.20486+01- 
4.66236  1-011 

-9.99996  1-01 

-1.10826+01  + 
1.51 28690 11 

719  (PIO) 

-1.10906+01+ 

1.5150o+0li 

4.18446i()0» 

2.76726+011 

3.82006  lOl 

-2.372.6.H01 

-1.10826+01 

1.51286+011 

7.18  (Pl8^ 

-i.l096c+01- 

1.5150c+01i 

4.18440+00- 

2.7G72e+0h 

-2.38586+0  1 

-1.41656  1-01 

-2.37246+01 

Z17  (piv) 

-2.37456+01 

-2.37316+01 

-1.1245.M-01  + 
1.25166-1011 

-3,679369  00+ 
7.72396+OOi 

-1.54146+01 

716  (PlG^ 

-1.5330ef01 

-9.85206+001 

2.5888o+00i 

■1.12456+01- 

1.25l6oi{)li 

-3.6793r;+00- 

7.7239c->-00i 

-3.68296+00+ 

7  70756+001 

715  ^Pis) 

-1.11966+01 

-9,85206+00- 

2.58886+001 

1.5306c  1  01 

-9.00166+(K)i 

2.93296+001 

-3.68296+00- 

7.70756+001 

714  (P14) 

-3.33016+00+ 

2.9883e+00i 

-3.24480+00+ 

6.5110c+00i 

-2.99276+00+ 

9.2275e-)0()l 

-9.U016e+00- 

2.93296+001 

-8.16826+00+ 

2.4218e+00i 

zi3  (Pia) 

-3. 3301 0, >-00- 
2.9883c+00i 

-3,24486+00- 
6.51  lOo+OOi 

-2.99276+00- 

9.227569-001 

•8.1682e+00- 

2.4218e+00i 

Z12  (Pl2) 

-3.56346+00+ 

2.9423c+00i 

-5.93276+00 

-3.5404c+()()+ 
2.97776 1  OOi 

-4.604869  00 

-2.00846+00 

711  (Pll) 

-3,56346+00- 

2.94236+001 

-2.04696+00 

-3.54046+00- 
2.9777. 'i-OOi 

-1.9,3676  1-00 

-3.05816+00  + 
5  4()67.'  0.3i 

7l0  (PlO^ 

-3,()39(l6+()U 

.  _  . . 

-2.98596+00 

-3.030.86100 

.3.(),5H1<' I'OO-  ' 
5  ■100.7, •  f),3i 

Zfj  (P9) 

-3.05706 1  00  t- 
1.250(;{'  Oi.i 

-3.()520ri  00 

■.3.0;  -’06  I  (K.i 

.3. 052,1). ■ )  00 

3  Oi.liOciOU 

78  (P8^ 

-3.05201H  00- 
1  25066-05i 

■;i.05),'H'  ■  00 

3. 01. 206  id! ! 

3.0;',2!-r(00 

)  ..3,34  ,(  02 

Z?  (py) 

■2.03636-Oi 

-1.3502.' 1^2 

1.33  18.'  02 

1.33766  02 

-2.5  1 57+  (j3 

76  (Po) 

1.33-176  1)2  1  2.01556  01 

■  1.95251'  (1  1 

■11  1966  04 

7.60,5.36-05 

75  (Ps) 

-2,52286-03 

-4  70‘^26  03 

-2.5)396  (13 

-4  693  1.' 0.3 

-2  03416-01' 

Z4  (P4) 

-2,53576-02 

-4.67256-03 

-4. 07456-0.3 

■4.67256-03 

-2.53566-02 

73  (P3) 

-2.53596-02 

-2.67  126-03 

-1.67256  0,3 

2.53576-02 

-2.53576-02 

72  ^P2) 

-4,6739e-03 

-2.53576  02 

-2,53576  02 

-2.53566-02 

-4 .67256-03 

Zl  (pl) 

-4.67256-03 

•2. 53;, (,6-02 

2  ;i3r)(;6-02 

-2.52 166-03 

-4.67256-03 

Table  A4.3  Reduced  Order  Low-a  Outer  Loop  Controller 


numll 

numl2 

num21 

num22  den 

K 

-9.1780e-03 

-1.3913e-04 

-5.2125e-04 

-7.9018e-06 

Zl5  (Pl5) 

2.8279e+05 

-5.97486405 

9.18816404 

5.16926408 

-7.07826+02 

Zu  (pi4) 

-7.1262e402 

-6.6096e402 

-7.07826402 

-6.9537e+02 

Zl3  (P13) 

-8.7801e+01 

-8.11996401 

-1.22556402 

-4.65906401+ 

4.648064011 

-8.7837e+01 

Z12  (P12) 

-4.5073e+01+ 

4.5590e+01i 

-4.87896401+ 

4.6977e401i 

-3.70436401 

+5.480264011 

-8.11396+01 

Zll  (pii) 

-4.5073e+01' 

4.5590e+01i 

-4.87896401- 

4.G977e401i 

-3.70436401- 

5.480264011 

-8.11396401 

-2.13176+01+ 

2.43416+011 

ZlO  (Pio) 

-2.1323e+01+ 

2.4620e+01i 

6.7676e-0l4 

2.6868e401i 

-8.9503e-01 

+4.460564011 

-5.49336400+ 

1.106964011 

-2.13176+01- 

2.46416+011 

Zg  (pg) 

-2.1323e+01- 

2.4620e+01i 

6.7676e-01- 

2.6868e401i 

-8.9503e-01- 

4.4605e401i 

-5.49336400- 

1.10696+011 

-5.37106+00+ 

1.0982e+01i 

Z8  (Ps) 

-5.0555e+00+ 

3.6437e+00i 

1.748064004 

5.0408e400i 

-1.76326401+ 

1.708464011 

-1.2468C+01 

-5.37106+00- 

1.09826+011 

2?  (P?) 

-5.05556400- 

3.6437e400i 

1.74806400- 

5.0408e400i 

-1.76326401- 

1.708464011 

-4.62696+00+ 

1.7637e+00i 

-5.20056+00+ 

9.0751.6-011 

26  (Pe) 

-1.40846400+ 

3.443964001 

-3.39986400+ 

1.237264001 

-5.21336400 

+3.709964001 

-4.6269e+00- 

1.76376+001 

-5.2005e+00- 

9.07516-011 

Z5  (Ps) 

-1.40846400- 

3.4439e400i 

-3.39986400- 
1 .237264001 

-5.21336400- 

3.709964001 

-2.5910e+00 

-3.32486-02 

Z4  (P4) 

-3.32i2e-02 

1.1408e-02 

-3.3213e-02 

-8.9760e-03 

-8.20886-04+ 

1.86146-041 

23  (P3) 

-8.9757e-03 

-8,96396-03 

-8.97526-03 

-2.35736-03 

-8.20886-04- 

1.8614e-041 

22  (P2) 

-7.3952e-04 

-3,3448e-02 

-7.12756-03 

-7.17106-03 

■7.12736-03 

21  (pl) 

-7.1264e-03 

-7.0672e-04 

-3.32186-02 

-8,97446-03 

Table  A4.4  Reduced  Order  High-a  Outer  Loop  (^ontroIhM- 


numll 

numl2 

nuni21  1 

num22 

den 

K 

-2.0973e-04 

4.4100e-03 

mkmmm 

r  -5.7H49e-02 

Zl2  (P12) 

9.01756+05 

3.78846+03 

2.49206+03 

8.10136+03 

-9.35156+02 

zil  (Pll) 

■^lESSiSi 

-6.9228e+02 

-6.7712e-^02 

-9.339  le+02 

-8.2879e+02 

zio  (Pio) 

-5.59916+01 

1 

-5.45756+01 

-5.0619e+01 

-2.09116+01  + 
1.1223e+01i 

-7.98526+00+ 

].7716c+01i 

Zg  (pg) 

-8.01426+00+ 

1.7728e+01i 

5.40516+00+ 

2.6664e+0H 

3.69316+01 

-2.09116+01- 
1.1 2236+0  li 

!  -7.98526+00- 
1.77166+Oli 

28  (Ps) 

-8.01426+00- 

1.7728e+01i 

5.40516+00- 

2.6664e+01i 

-5.31056+00+ 

1.4401e+01i 

-3.98236+00+ 

8.1873e+00i 

-3.97416+00+ 

8.1974e+00i 

Z7  (P7) 

-3.14266+00+ 

7.0307e+00i 

-5.3105e+00- 
1.4401 6+0  li 

-3.9823e+00- 

8.1873e+00i 

■  3.9741 6+00- 
8.1974e+()0i 

26  (Pg) 

-2.38536+00+ 

2.9159e+00i 

-3.14266+00- 

7.0307e+00i 

-4.17286+00+ 

6.0155e+00i 

-2.41936+00+ 
6.61666-0  li 

-2.l405e+00 

25  (Ps) 

-2.38536-f-OO- 

2,9159e+00i 

-1.6954e+00 

I 

-4.17286+00- 
6.0l55e  fOOi 

-2.4193e+00- 

6.6I666-OI1 

-2.6792e02 

Z4  (P4) 

-2.6804e-02 

-2.71936-02 

-2.5690C-03 

-2.G817e-02 

-5.8955e-03 

23  (P3) 

-2.57946-03 

-2.7396e-03 

-5.9017e-03 

-2.57956-03 

-2.034.16-04 

22  (P2) 

-5.8947e-03 

-5.92836-03 

-5.9316e-03 

-7.6053e-05 

2l  (Pl) 

-2.08636-04 

2.01556-04 

!  -1.95256-04 

-1.11956-04 

-2.57156-03 

5  SUPERMANJi:UVERABLE  F-18  DESIGN 


This  section  presents  the  design  of  a  full  envelope  flight  control  system  for  a 
supermaneuverabie  F-18  aircraft.  The  purpose  of  the  design  is  to  take  advantage  of 
advanced  linear  control  design  techniques  while  addressing  practical 
implementation  issues.  A  description  of  the  aircraft  model  is  presented,  followed  by 
detailed  descriptions  of  the  flight  control  system  design  and  analysis.  The  specific 
numbers  used  in  the  design  are  indicated  in  Appendix  5. 

5.1  MODEL  DESCRIPTION 


The  supermaneuverabie  aircraft  model  described  in  this  section  is  based  upon 
a  modified  version  of  the  F-18  aircraft.  The  supermaneuverabie  aircraft  is  a  twin 
engine  fighter  aircraft  with  a  moderately  swept  wing,  twin  canted  vertical  tails,  and 
a  large  leading  edge  root  extension.  Leading  and  trailing  edge  flaps  are  used  to 
control  the  wing  camber  for  maximum  maneuvering  perfomiance  throughout  the 
flight  envelope.  The  horizontal  tail  is  located  below  the  wing  to  provide  increased 
longitudinal  stability  at  high  angles  of  attack.  The  twin  vertical  tails  provide 
positive  directional  stability  beyond  maximum  trimmed  angles  of  attack.  An  F-18 
aircraft  is  shown  in  Fig.  5.1. 


The  aircraft  model  is  augmented  with  two-dimensional  thrust  vectoring 
nozzles  that  provide  pitch  and  yaw  moments  when  deflected  symmetrically  and  a 
roll  moment  when  deflected  asymmetrically.  The  aerodynamic  control  inputs  to  the 
aircraft  dynamics  are  the  elevators,  the  ailerons,  the  rudders,  and  the  leading  and 
trailing  edge  flaps.  The  aerodynamic  surfaces  are  useful  at  normal  flight 
conditions,  v/here  there  is  adequate  aerodynamic  control  surface  effectiveness.  The 
thrust  vectoring  inputs  are  useful  at  high  angle  of  attack,  low  dynamic  pressure 
operating  conditions,  where  the  traditional  aerodynamic  control  effectiveness  is 
inadequate.  The  pilot  inputs  include  a  control  stick  and  rudder  pedals. 
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P’ig.  5.1  F- 18  Aircraft 


5.1.1  N0NLINK4R  MODEL 

The  aircraft  model  is  a  nonlinear  FORTRAN  model.  The  model  consists  of 
separate  modules  describing  the  atmosphere,  nonlinear  equations  of  motion, 
aerodynamics,  engines  including  thrust  vectoring  noz^les,  variable  geometry  inlets, 
sensors,  and  actuators  which  include  rate  and  position  limits.  The  high-fidelity 
model  was  developed  in  L51],  which  gives  more  detail  than  presented  here.  This 
section  is  presented  for  completeness;  however,  it  is  only  a  summary  of  the  model 
description  given  in  [5.1].  The  nonlinear  model  structure  of  the  supermaneuverable 
aircraft  is  the  same  as  the  F-16  model  structure  shown  in  Fig.  4,2. 

The  central  component  of  the  nonlinear  model  contains  the  six  clegree-of- 
freedom  aircraft  equations  of  motion  given  by  eqs.  (2.1)  and  (2.2).  A  complete 
description  of  the  equations  of  motion  and  atmospheric  components  is  given  in 
section  2.1.  The  primary  inputs  to  the  equations  of  motion  are  the  forces  and 
moments  generated  by  the  aerodynamic,  propulsion,  and  inlet  systems.  'I'he 
primary  outputs  are  the  state  variables  of  the  aircraft. 
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The  aerodynamic  actuator  component  of  the  nonlinear  model  transforms 
actuator  commands  into  aerodynamic  surface  positions.  There  are  five  pairs  of 
aerodynamic  surfaces;  three  pairs  for  active  control  and  two  pair  scheduled  for 
optimum  performance.  The  ailerons,  rudders,  and  elevators  are  used  for  stability 
augmentation  and  flight  path  manipulation.  The  leading  and  trailing  edge  flaps  are 
scheduled  versus  angle  of  attack,  Mach  number,  and  dynamic  pressure  to  provide 
optimal  lift-to-drag  ratio,  improve  high  angle  of  attack  stability,  and  alleviate  air 
loads  on  the  wing  during  high-g  maneuvers.  All  aerodynamic  surfaces  are  powered 
by  hydro-mechanical  servo-actuators.  The  leading  edge  flaps  use  rotary  hydraulic 
drive  motors  to  power  mechanical  actuators  while  all  other  surfaces  use  linear 
piston  actuators.  The  dynamic  models  of  the  aerodynamic  surface  actuators  and 
corresponding  position  and  rate  limits  are  given  in  Appendix  5.  The  resulting 
aerodynamic  surface  positions  from  the  actuator  models  are  input  to  the 
aerodynamic  model. 

Baseline  aerodynamic  data  are  obtained  from  wind  tunnel  and  flight  test 
data  covering  Mach  numbers  ranging  from  0.2  to  2,  altitudes  from  sea  level  to 
60,000  feet,  and  angles  of  attack  and  sideslip  ranging  from  -12  to  90  degrees  and  -20 
to  20  degrees  respectively.  Thrust  vectoring-induced  aerodynamic  data  are  obtained 
from  nonaxisymmetric  nozzle  wind  tunnel  tests  and  predictive  jet  interference 
analytical  methods.  Linear  extrapolation  is  used  to  define  aircraft  properties 
outside  of  the  ranges  of  the  wind  tunnel/flight  test  model.  The  aerodynamic  data  are 
contained  in  tabular  format  and  linear  interpolation  is  used  for  traditional  force  and 
moment  aerodynamic  coefficient  build-up.  Thrust  vectoring-induced  aerodynamic 
effects  contributions  are  added  to  static  and  dynamic  baseline  aerodynamic 
coefficients  to  obtain  total  aerodynamic  coefficients.  The  total  coefficients  are 
translated  from  the  stability  axes  to  the  body  axes  and  modified  to  account  for 
center  of  gravity  offset. 

A  modified  engine  model  is  used  to  provide  increased  thrust  to  the  modified 
F-18  airframe  model,  and  the  airflow  is  scaled  to  match  standard  F404-400  engine 
airflow.  The  model  represents  the  controlled  closed-loop  response  of  the  nozzle  and 
engine.  Engine  dynamics  are  modeled  as  a  first  order  lag  with  a  time  constant  of 
0.8  second  and  rate  limits  of  ±22  deg/sec  to  model  engine  control 
acceleration/deceleration  schedules. 
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Symmetric  and  asymmetric  motion  of  axisymmetric  gimbaled  nozzles  pro^^ide 
roll,  pitch,  and  yaw  moments  by  vectoring  the  engine  exhaust.  Over/underturn 
effects  of  the  vectored  flow  are  modeled,  and  the  gimbaled  nozzle  actuators  are 
assumed  to  have  second  order  dynamics  with  position  and  rate  limits.  The  nozzle 
actuator  models  are  given  in  Appendix  5.  The  inputs  to  the  propulsion  model  are 
thrust  commands,  flight  conditions,  and  thrust  vectoring  commands,  and  the 
outputs  are  propulsion  forces  and  moments.  Thrust  and  thrust  vectoring  forces  and 
moments  are  transformed  into  the  body  axes  and  center  of  gravity  offset  is  taken 
into  account. 


The  F-18  model  is  modified  with  a  high  performance  variable  geometry  inlet 
model.  The  inlet  model  consists  of  a  fixed  ramp  and  an  adjustable  second  ramp  to 
provide  external  supersonic  flow  compression  at  olf-design  conditions.  Auxiliary  by¬ 
pass  and  inlet  doors  in  the  subsonic  diffuser  provide  additional  airflow  control  at 
low  speed,  high  airflow  conditions.  The  inputs  to  the  inlet  model  arc  the  flight 
parameters  and  the  engine  airflow  command.  The  main  outputs  are  the  inlet  forces 
and  moments  due  to  pressure  recovery  and  ram  and  spillage  drag.  All  resultant 
inlet  forces  and  moments  are  resolved  to  body  axes  and  modified  for  center  of 
gravity  offset. 


The  sensor  component  includes  a  rate  sensor  unit,  an  accelerometer  sensor 
unit,  an  jur  data  unit,  and  an  inertial  sensor  set.  The  rate  gyro  sensors  provide  roll, 
pitch,  and  yaw  rate  measurements  in  the  body  axes,  i'hey  are  assumed  to  be 
perfectly  accurate  and  have  infinitely  fast  dynamics  since  the  actual  .sensor 
dynamics  are  well  beyond  the  control  system  bandwidtli.  The  acceloromoters 
provide  acceleration  measurements  at  the  sensor  locations,  and  they  are  also 
assumed  to  be  perfect  with  infinitely  fast  dynamics,  i’he  air  data  unit  measures 
dynamic  and  ambient  pressure  and  local  angles  of  attack  and  sideslip  using  pitot- 
static  probes  and  aerodynamic  vanes  respectively.  The  inertial  sensor  unit  provides 
angular  orientation  and  inertial  po.sition  and  velocity  measurements.  Air  data  are 
augmented  with  inertial  data  to  generate  actual  angle  of  attack  and  sideslip 
feedback  signals.  Kalman  filtering  is  used  to  account  for  local  wind  ei'fects  on  the 
air  data.  Again,  the  air  data  and  inertial  sensor  units  are  assumed  to  provide 
perfect  measurements  with  infinitely  fast  dynamics. 


Dynamics  of  the  primary  structural  modes  are  modeled  to  generate 
incremental  accelerations  and  rotational  rates  at  the  sensor  locations.  The  four 
modeled  structural  modes  are  the  first  wing  symmetric  and  asymmetric  modes  and 
the  first  vertical  and  lateral  fuselage  modes.  The  inputs  to  the  structural  model  are 
normal,  lateral,  and  roll  accelerations.  Each  modeled  structural  mode  is 
represented  by  NASTRAN  generated  second  order  dynamic  models. 

5.1.2  LINEAR  MODEL 

The  linear  control  law  design  process  begins  by  generating  a  set  of  linear 
models  at  different  trimmed  flight  conditions.  As  discussed  in  section  2.1,  the 
nonlinear  model  is  trimmed  at  a  fixed  flight  condition  and  linearized  by  perturbing 
the  model.  The  resulting  linear  equations  representing  the  longitudinal  and 
lateral/directional  dynamics  are  given  by  eqs.  (2.11)  and  (2.14),  respectively.  Linear 
models  generated  at  different  flight  conditions  allow  the  use  of  well  established 
linear  design  methods  and  represent  the  aircraft  dynamics  at  a  range  of  Mach 
numbers  and  altitudes,  and  therefore  dynamic  pressures.  The  points  chosen  for 
linear  model  generation  must  represent  a  wide  range  of  dynamic  pressure,  since 
dynamic  pressure  is  used  in  the  scheduling  of  the  inner  loop  feedback  gains. 

This  chapter  addiesscs  robust  stabilization  for  the  full  conventional  flight 
envelope  of  the  supermaneuverable  aircraft  described  preriously.  The  conventional 
flight  envelope  is  defined  as  the  subsonic  Mach  regime  between  sea  level  and  40,000 
ft.  Dynamic  pressure  ranges  from  50  psf  to  1000  psf.  These  limitations  correspond 
to  angles  of  attack  up  to  approximately  25  deg.  Fig.  5.2  represents  the  flight 
envelope  considered  for  flight  control  design. 

The  flight  conditions  chosen  for  linear  longitudinal  axis  design  are  chosen 
such  that  all  linear  models  represent  the  full  range  of  dynamic  pressures  for  the  full 
conventional  envelope.  Likewise,  the  lateral/directional  linear  models  are  generated 
to  cover  the  same  dynamic  pressure  range  but  at  different  Mach  number  and 
altitude  combinations  than  those  of  the  longitudinal  models.  Different  design 
conditions  are  chosen  for  the  longitudinal  design  and  the  lateral/directional  design 
to  demonstrate  the  lack  of  meaning  of  "design  condition"  for  full  envelope  control 
law  design.  One  common  set  of  design  conditions  could  just  as  easily  have  been 
chosen  for  both  the  longitudinal  and  lateral/directional  controller  designs. 
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Fig,  5,2  Flight  Envelope 

Appendix  5  lists  the  flight  conditions  depicted  in  Fig.  5.2  and  information  about 
them. 

5.2  CONTROL  SELECTOR  DESIGN 

As  presented  in  section  3.1,  the  control  selecto;  is  a  mapping,  given  by  eq. 
(3.2),  from  the  generalized  controls  to  the  actual  control  effectors.  The  generalized 
and  actual  controls  for  the  supermaneuverable  vehicle  are  given  by 
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where  p  is  the  roll  acceleration,  q  is  the  pitch  acceleration,  f  is  the  yaw  acceleration, 
5e  is  the  S3Tnmetric  elevator  position,  8dt  is  the  asymmetric  elevator  position,  8a  is 
the  aileron  position,  5r  is  the  rudder  position,  5pTV  is  the  symmetric  pitch  thrust 
vectoring  nozzle  position,  Srtv  is  the  asymmetric  pitch  (roll)  thrust  vectoring  nozzle 
position,  and  Syrv  is  the  yaw  thrust  vectoring  nozzle  position. 


Using  eq.  (5. 1),  the  linear  dynamics  of  eqs.  (2, 11)  and  (2.14)  are  written,  as: 
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The  dynamics  of  the  body-axis  rates  are  a  subset  of  eq.  (5.2)  and  are  given  by: 
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Consider  the  following  partitioning  of  the  control  effector  vector: 
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Now  eq.  (5.3)  becomes 
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rates,  eq.  (3.1)  is  written  as; 
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5  is  given  by  eq.  (5.1),  and  the  generalized  control  effectiveness  matrix,  B*,  is 
identity  since  the  generalized  controls  correspond  to  the  rotational  degrees-of- 
freedom  accelerations. 


A  daisy-chain  method  is  used  to  generate  thrust  vector  commands.  Thrust 
vectoring  is  used  only  when  the  aerodynamic  surfaces  are  not  able  to  generate  the 
necessary  forces  and  moments  required  for  commanded  maneuvers.  Therefore,  the 
computation  of  aerodynamic  control  commands  is  independent  of  thrust  vectoring 
control  commands.  The  control  selector  is  then  defined  using  eqs.  (3.2)  and  (3.3) 
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where  and  are  used  to  prioritize  the  use  of  redundant  control  effectors. 
Since  the  aile.rons  c<'nt.ribute  more  to  the  roll  acceleration  and  the  first  priority  of 
the  horizontal  taii  should  be  pitch  control,  the  differential  horizontal  tail  command 
is  reduced  by  weighting  the  command  to  be  a  quarter  of  the  other  aerodynamic 
commands.  There  is  no  redundancy  for  the  thrust  vectoring  control  effectors,  and 
thus,  the  prioritiz  uion  matrices  become 
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.001. 

Computation  of  the  control  .selector  equation  (5.9)  depends  on  flight  condition. 
Therefore,  elements  of  and  are  found  using  linear  interpolation  of 

stored  table  values. 


Nonlinear  elements,  such  as  position  and  rate  limits,  are  required  to 
implement  the  daisy-chain  thus  making  the  control  selector  nonlinear.  Fig.  5.3 
shows  the  structure  of  the  nonlinear  control  selector. 
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Fig.  5.3  Nonlinear  Control  Selector 


A  limited  aerodynamic  surface  command  is  generated  from  a  rotational 

acceleration  command  (6  )  via  the  aerodynamic  control  selector  (Tae,.^)  and 

aerodynamic  surface  rate  and  position  limits.  An  acihevable  rotational  acceleration 
* 

vector  (djjvail)  is  computed  from  the  limited  aerodynamic  surface  command  using  the 
aerodynamic  control  distribution  difference  of  the  commanded  and 

achievable  rotational  acceleration  vectors  (es*)  is  transformed  to  a  thrust  vector 
command  (6^ygg)  using  the  thrust  vector  control  selector  (T^ygc). 

The  bottom  line  of  transformation  to  generalized  controls  is  seo'^  by 
combining  the  state  dynamics  and  the  control  selector.  The  resulting  longiti 
and  lateral/directional  linear  state  dynamical  equations  arc: 

(; 


I 


(5.12) 


based  on  the  above  control  selector  definition.  Therefore,  eqs.  (5.11;  and  (5.12) 
represent  the  form  of  the  models  used  for  linear  control  design. 


5.3  LONGITUDINAL  AXIS  CONTROLLER 


The  longitudinal  axis  manual  flight  control  system  generates  elevator  and 
pitch  thrust  vector  commands  from  angle  of  attack  commands.  These  control 
effector  commands  stabilize  the  short  period  aircraft  dynamics  across  a  wide 
conventional  flight  envelope.  The  open-loop  linear  .“short  period  output-state 
equation  is  defined  by; 


a 

-  q  - 

—  Along 

a 

-  q  - 

+  Biong 

r  “1 

a 

-  q  - 

=  f^long 

a  ' 

-  q  - 

with 

5e 


Glong  —  (Along, Blong»C]ong>0)> 


ong)^long»'“'long> 


(5.13) 


where  Along  and  Biong  are  defined  in  (2.11),  and  Ciong  is  obviously  identity.  Taking 
into  account  the  function  of  the  control  selector,  eq.  (5. 13)  becomes 
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where  B]g„g  is  defined  in  eq.  (5.11). 

A  minimal-order  design  algorithm  is  used  for  the  desicn  of  an  inner  loop 
equalization  controller.  Structured  singular  value  synthesis  is  used  to  design  outer 
loop  robust  performance  controllers.  Different  control  laws  are  found  for  high  and 
low  dynamic  pressure  conditions,  and  controller  commands  are  blended  for  a  small 
region  of  dynamic  pressure.  The  inner  and  outer  loop  control  designs  are  presented 
followed  by  robustness  and  flying  qualities  analyses  for  the  entire  longitudinal  axis 
flight  controller, 
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5.3.1  INNER  LOOP  DESIGN 


The  structure  of  the  aircraft  with  a  closed  inner  loop  controller  is  shown  in 
Fig.  5.4. 


P 


The  pitch  acceleration  command,  Qc,  is  made  up  of  an  inner  loop  pitch  acceleration 
command,  and  an  outer  loop  pitch  acceleration  command,  P  is  a 

linear  model  of  the  inner  equalization  controller,  Ke,|,  and  a  combination  of  the 

open-loop  aircraft  dynamics  and  the  control  selector,  Gjong-  Dependency  of  an 
element  on  flight  condition  is  represented  by  (M,h). 

The  purpose  of  the  inner  loop  equalization  controller  is,  as  the  name  suggests, 
to  equalize  the  transfer  function  between  angle  of  attack  response  and  outer  loop 
pitch  acceleration  command  across  the  flight  envelope.  In  other  words,  the  goal  of 
the  inner  loop  is  to  I’educo  the  variation  of  the  aircraft  dyi^amic.s  between  operating 
conditions,  thereby  reducing  the  modeling  uncertainty  between  flight  condition 
dependent  aircraft  models.  Successfully  reducing  the  modeling  uncertainty 
provides  a  greater  possibility  that  one  rohnst  performance  controller  will  provide 
stability  and  performance  across  all  operating  conditions.  Therefore  an  effective 
inner  loop  controller  must  be  designed  (o  reduce  the  relative  error  between  P  and  Py 
as  defined  by  eq.  (3.4). 

Since  P  is  a  function  of  flight  condition,  the  equalization  module  must  be 
scheduled  against  slowly-varying  flight  condition-dependent  parameters  L5.2j,  such 
as  dynamic  pi'essure,  altitude,  or  Mach  number.  Cain  scheduling  of  the  inner  loop 
controller  is  necessary  because  the  aircraft  dynamics  vary  greatly  across  a  wide 


envelope,  making  it  impossible  to  design  a  single  robust  controller.  In  this  design, 
the  inner  loop  equalization  controller  gains  are  scheduled  across  the  flight  envelope 
as  functions  of  dynamic  pressure. 


The  minimal-order  Hoc  design  algorithm  described  in  section  2.8  is  used  to 
design  a  longitudinal  axis  inner  loop  equalization  controller  for  the  design  model 
shown  in  Fig.  5.5. 


Fig.  5.5  Inner  Loop  Design  Model 


A  frequency  dependent  weight,  Wa,  is  used  to  obtain  the  desired  response  of  the 
angle  of  attack  to  pitch  acceleration  commands.  If  the  design  weight  has  the 
realization  Wa  =  (Aw,Bw,Cw,Dw),  then  the  parameters  of  the  design  model  from 
Fig.  2.10  become 


A  = 

Along 

,  B  ^ 

* 

^long 

- 1 

< 

0 

_ J 

.  0  . 

C  =  [  Ciong  0  ], 


Gi- 


^lonK^W 


B 


w 


,  0], 


■  inner 

u  =  Qc  ,  X  =  y  = 


a 

.qj 


.  outer 

,  Wi  =  q^.  ,  zi  =  tt, 


^ri.l5) 


and  H2  is  the  control  signal  weighting  used  as  a  design  parameter.  Since  the 
control  system  will  eventually  track  angle  of  attack  commands,  the  performance 
output  distribution  matrix  is 

C*  =  [  1  0  ].  (5.16) 


For  the  flight  conditions  of  Fig.  5.2,  the  magnitudes  of  the  open  loop  pitch 
acceleration  to  angle  of  attack  transfer  functions  are  shown  in  Fig.  5.6.  There  is  a 
wide  variation  in  aircraft  dynamics  across  the  flight  envelope.  The  high  dynamic 


pressure  flight  condition  models  have  lower  magnitudes  at  low  frequency  than  the 
low  dynamic  pressure  flight  condition  models.  Since  the  purpose  of  the  inner  loop 
controller  1 ;  to  equalize  these  dynamics,  the  design  weight,  W^,  is  chosen  as  the 
inverse  of  some  transfer  function  that  lies  among  the  span  of  open  loop  dynamics.  It 
is  suggested  that  the  inverse  design  weight  dynamics  be  chosen  close  to  the 
dynamics  of  the  high  dynamic  pressure  flight  conditions  so  that  positive  feedback  is 
not  used  to  diminish  the  fast  dynamics  of  the  high  dynamic  pressure  conditions. 


Fig.  5,6  Pitch  Acceleration  to  Angle  of  Attack  Open  Loop  Dynamics 

The  weights  used  for  the  inner  equalization  design  are 

Aw  =100  Bw  =  l  Cw  = '45918  1)^  =  51020 

Hg  =  .035  .  (5.17) 


With  the  design  weight  augmented  to  the  open-loop  aircraft  dynamics,  the 
number  of  design  states  become.s  three,  and  the  numb(o-  of  inputs  and  outputs 
remains  at  one  and  two  respectively.  Therefore,  the  order  of  the  equalization 


controller  is  one  (#  of  states  -  #  of  measurements).  The  resulting  controller 
parameters,  with  design  parameter  a  =  1  and  design  weights  given  by  eq.  (5.17),  are 

F  =  -40  K^=:[l  1]  T  =  [0  .0247  Oj 

N  =  N(q)  M  =  [Mi(q)  M2(q)].  (5,18) 

Only  three  control  parameters  require  scheduling  with  dynamic  pressure  (q).  The 
elements  ofN  and  M  in  eq.  (5.18)  are  fit  using  least-squares  polynomial  curve-fits, 
linear  in  dynamic  pressure; 

N(q)  =  -.312  q-(-461 

M^fq)  = -.058  q  +  50.5  = -.006  q -t- 8.11  (5.19) 

The  controller  given  by  eqs.  (5.18)  and  (5.19)  provides  closely  matched  frequency 
responses  of  the  closed  inner-loop  transfer  function  (P)  from  to  a  across  the 
flight  envelope  as  shown  in  Big.  5.7. 

It  is  assumed  that  equalizing  the  pitch  acceleration  to  angle  of  attack 
dynamics  will  also  equalize  the  pitch  acceleration  to  pitch  rate  dynamics  since  they 
are  not  independent.  The  open  loop  pitch  acceleration  to  pitch  rate  transfer 
functions  are  shown  in  Fig.  5.8.  Closing  the  inner  loop  with  the  controller  given  by 
eqs,  (5.18)  and  (5.19)  approximately  equalizes  the  pitch  acceleration  to  pitch  rate 
dynamics  as  shown  in  Fig,  5.9. 

To  analyze  the  variation  of  aircraft  dynamics  at  different  flight  conditions, 
the  relative  error  defined  by  eq.  (3.4)  is  computed  for  a  given  linear  model  (P)  aiad  a 
nominal  linear  model  (Pq).  In  order  to  evaluate  the  relative  error  between  the 
equalized  models  of  Fig.  5.7,  the  closed  inner-loop  model  (P)  at  Mach=0.95  and 
altitude=20  kft  is  arbitrarily  chosen  as  the  nominal  model  (Pq).  Fig.  5.10  shows  the 
maximum  singular  values  of  the  model  errors  relative  to  Pq  across  the  flight 
envelope  at  design  conditions  given  in  Fig.  5.2.  Since  tlie  closed  inner-loop  errors 
relative  to  the  Mach=0.95  and  altitude=20kft  flight  condition  are  less  than  unity  for 
all  frequencies  and  across  the  flight  envelope,  Pq  is  used  as  the  open  outer-loop 
design  plant  for  the  robust  performance  controller. 
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Fig.  5.7  Pitch  Acceleration  to  Angle  of  Attack  Closed  Inner-Loop  Dynamics 


Fig.  5.8  Pitch  Acceleration  to  Pitch  Rate  Open  Loop  Dynamics 


5.3.2  OUTER  LOOP  DESIGN 

The  structure  of  the  system  with  the  outer  loop  controller  is  shown  in  Fig. 

6.11. 


Fig.  5.11  Outer  Loop  Control 


Po  is  the  closed  inner  loop  nominal  model  (or  the  open  outer  loop  model),  and  the 
robust  outer  loop  performance  controller  is  Kj,.  The  purpose  of  the  outer  loop 

performance  controller  is  to  generate  a  pitch  acceleration  command  from  an 

angle  of  attack  command  (ac),  that  produce.s  a  desired  robust  aircraft  angle  of  attack 
response  (a). 


This  problem  is  formulated  as  an  implicit  model-following  problem  [5.3], 
where  the  ideal  model  to  be  followed  is  chosen  to  be  a  second-order  transfer  function 


*^ideal  _  _ _  ...  ^ideal 

S  +  2  ®  ■)'  ^deal 


(5.20) 


with  and  chosen  based  on  fl3dng  qualitie.s  requirements  from  section  2.2. 
Using  this  ideal  model-following  approacli,  the  outer  loop  porformance  controller 
design  model  is  developed  in  Fig.  5.12. 


Fig.  5.12  Outer  Loop  Design  Model 


Po  is  the  design  plant  described  previously,  and  Wp  and  are  performance  and 
input  uncertainty  design  weights,  respectively. 

Structured  singular  value  (|i)  synthesis  is  used  to  design  outer  loop 
performance  controllr  rs  for  the  system  in  Fig.  5.12.  As  described  in  section  2.9, 
p-synthesis  is  a  combination  of  general  Hao  feedback  design  and  structured  singular 
value  analysis.  Therefore,  the  design  model  of  Fig.  5.12  is  organized  into  the 
structure  of  the  standard  Ho«  design  model  of  Fig.  2.9  and  the  standard  structured 
singular  value  M-A  format  of  Fig.  2.7.  The  resulting  p-synthesis  model  is  shown  in 
Fig.  5.13. 


Fig.  5.13  Structured  Singular  Value  Syrithesis  Model 


Pw  is  the  weighted  outer  loop  design  plant,  or  the  weighted  Pq. 


The  robust  performance  controller  (Kjj)  for  Pw  is  designed  using  p-synthesis 
with  the  performance  and  input  uncertainty  weights  set  as  follows: 


Wp 


25s  +  50 
s  +  5 


10s  +  1000 
-  s  +  10000  • 


(5.21) 


The  process  o.  choosing  these  design  weights  is  somewhat  empirical.  The  tracking 
of  the  ideal  angle  of  attack  response  is  impf'rtant  at  frequencies  most  sensitive  to 
the  pilot  (low  frequencies),  whereas  robustness  to  input  uncertainties  is  important 
at  different  frequencies  (high  frequencies).  Therefore,  the  performance  weight  is 
chosen  as  a  low-pass  filter,  and  the  input  uncertainty  weight  is  chosen  as  a  high- 
pass  filter.  The  selection  of  the  poles,  zeros,  and  gains  of  the  filters  is  iterative. 
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Initial  filter  values  are  chosen,  a  controller  is  designed,  and  time  responses  are 
examined.  The  filter  values  are  tuned  until  satisfactory  tracking  occurs  with 
minimal  control  effort. 

Unlike  responses  in  the  lateral/directional  axes,  the  desired  pitch  response  is 
not  unifonn  across  the  flight  envelope.  The  pilot  would  like  to  feel  a  faster  response 
at  higher  speeds.  To  account  for  the  different  performance  requirements,  two 
robust  performance  controllers  are  designed  using  two  different  ideal  models  based 
on  eq.  (5.20):  1)  a  low  dynamic  pressure  controller  using  a  slow  ideal  model 

(Cideai=0-S»  (0ideai=3  rad/sec),  and  2)  a  high  dynamic  presscu'e  controller  using  a  fast 
ideal  model  (Cideal=Q'8.  Widcal=5  rad/sec).  The  resulting  controllers  are  13th  order 
and  of  the  following  form: 

K(s  +  Zl)(s  4-  Z2)-(s  +  Z12)(s  +  zia) 

(s  +  pi)(s  +  P2)-(S  +  Pi2)(s  +  P13)  ' 

where  tlie  numerator  and  denominator  coefficients  are  given  in  Appendix  5. 
Implementation  of  the  two  outer-loop  perfomiance  controllers  is  discussed  later  in 
this  section. 

As  seen  by  eq.  (5.22),  p-synthesis  typically  generates  controllers  of  an  order 
much  higher  than  the  original  plant  because  of  the  design  weights  and  the 
frequency  dependent  scalings  discussed  in  section  2.9.  The  additional  state 
variables  in  the  controllers  are  a  result  of  the  design  method.  Therefore,  good 
performance  is  possible  with  a  reduced-oi’der  controller.  The  structured  singular 
values  of  the  closed  loop  system  of  Kig.  5.13  are  computed  with  the  13tli-order 
controller  implemented.  Fig.  5.14  and  Fig.  5.15  show  that  the  maximum  singular 
values  for  the  low  and  high  dynamic  pressure  controllers  are  0.95  and  0.85, 
respectively.  The  balanced  truncation  method  of  ^^ection  2.10  is  used  to  reduce  the 
controller  order  without  degrading  the  performance  and  robustness  of  the 
controller.  A  maximum  number  of  transformed  controller  states,  z,  of  eq.  (2.81),  are 
truncated  such  that  the  maximum  structured  singular  value  does  not  degrade  above 
unity  for  the  closed-loop  system  of  Fig.  5.13.  The  controller  is  reduced  from  13th 
order  to  4th  order  without  significant  performance  or  robustness  deterioration. 
Therefore,  using  the  notation  of  section  2.10,  the  dimension,  k,  of  the  truncated 
transformed  controller  state  vector,  Z2,  i^  9,  and  the  dimension  of  the  reduced-order 
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state  vector,  21,  is  4.  Fig.  5.14  shows  the  structured  singular  values  of  the  closed- 
loop  system  with  full-  and  reduced-order  low  dynamic  pressure  controllers. 


10-2  10-'  ,  ip»  10'  102 

0)  (rad/sec) 


Fig.  5.14  Closed- Loop  System  with  Low  q  Controller 


It  is  seen  that  reduction  of  the  full  13th-order  controller  to  a  4th-order  controller 
does  not  destroy  the  original  robustness  and  performance.  Pig.  5.15  shows  similar 
results  for  the  high  dynamic  pressure  controller,  The  final  4th-order  low  and  high 
dynamic  pressure  controllers  are 


lo  1.30x102(3  +  6.02xl02)(s  +  3.84x100  ±j6.25xl00) 

N  is  +  5.91x101  ±  j5.20xl0l)(s  +  3.63xl0-2)(,s  +  6.64xl00) 

N  1.20x102(3  2.42x103)(s  -h  5.4Qxl00  ±  j7.26xl00) 

~  (s  -f-  7.18x101  ±  j4.05xl0"l)(s  +  7.23xl0-2)(s  -f-  2.05xl00) ' 
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B’ig.  5.15  Closed-Loop  System  with  High  q  Controller 

It  is  known  that  at  high  dynamic  pressure  conditions,  l^'should  be  used,  and 

at  low  dynamic  pressure  conditions,  K|"should  be  used.  However,  the  definition  of 

high  and  low  dynamic  pressure  is  unclear.  The  multiple  dynamic  compensators  of 
eq.  (5.23)  are  implemented  through  a  blending  of  controller  commands  over  a  range 
of  dynamic  pressure.  B'ig.  5.16  schematically  shows  the  combination  of  controller 
commands  generated  from  the  high  and  low  dynamic  pressure  controllers  of  eq. 
(5.23). 
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B^ig.  5. 16  Outer  Loop  Controller  Blending 
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Implementation  of  the  controllers  is  achieved  through  a  simple  linear  blending 
parameter  (5b)  that  gener«-^s  a  combination  of  controller  commands  for  a  region  of 
dynamic  pressure.  Fig.  i^.l7  shows  the  blending  parameter  as  a  function  of  dynamic 
pressure. 


Fig.  5.17  Blending  Faramcter  Function 

The  dynamic  pressure  blending  region  between  170  psf  and  200  psf  and  the  piece- 
wise  linear  function  were  chosen  to  maximize  the  use  of  the  "fast"  high  dynamic 
pressure  controller.  Once  the  high  dynamic  pressure  controller  caused  an  elevator 
rate  saturation  foi  a  5  deg.  angle  of  attack  step  input  (-200  psf),  the  low  dynamic 
pressure  controller  command  was  blended  with  the  high  dynamic  pressure 
controller  command.  The  piece-wise  linear  function  was  chosen  for  simplicity,  and  a 
blending  region  of  about  30  psf  seemed  to  minimize  the  transient  blending  effects  in 
nonlinear  simulations. 

6.3.3  ROBUSTNESS  ANALYSIS 

Structured  singular  value  analysis  techniques  are  used  to  analyze  stability 
robustness  of  the  designed  control  system  to  uncertainties  corresponding  to 
unmodeled  actuator  and  sensor  dynamics,  parameters  in  the  plant  model,  and 
blending  of  high  and  low  dynamic  pressure  outer  performance  controllers. 

The  reduced  order  actuator  models,  given  in  Appendix  5,  are  used  in  the 
nonlinear  model  and  the  linear  analysis  model,  and  the  difference  between  the  high- 
order  and  low-order  models  represents  the  unmodeled  actuator  dynamics  [5.1].  The 
actual  error  dynamics  between  these  models  are  fit  to  real  rational  transfer 
functions  that  become  the  weighting  functions  for  the  actuator  uncertainty 
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(5.24) 


.63s^  +  3.03s  +  .078 

Wp  =  2 

s  +  68.4s  +  1900 


Wg  is  the  elevator  uncertainty  weight.  The  absc.ace  of  a  pitch  thrust  vector 
uncertainty  weight  is  due  to  the  fact  that  thrust  vectoring  is  used  only  during 
aerodynamic  effector  saturation.  Since  saturation  is  a  nonlinear  phenomenon, 
linear  robustness  analysis  is  inappropriate. 


The  angle  of  attack  and  pitch  rate  sensor  dynamics  are  captured  entirely  as 
unstructured  uncertainty.  The  sensor  dynamics  are  estimated  from  flight  test  data 
of  high  performance  aircraft  [5.1]  and  fit  to  real  rational  transfer  functions 
representing  measurement  uncertainty  weighting  functions 


_  21.9s^+  112QS  + 91100 
“  '  s\  574s  +  1140000 

.745s^  +  152s^'  +  95.9s  +  1.38 
^  3  +  626s  +  173000s  +  235000 


(5.26) 


Wq  is  the  angle  of  attack  sensor  uncertainty  weight,  and  is  the  pitch  rate  sensor 
uncertainty  weight. 


An  extensive  aerodynamic  uncertainty  database,  developed  from  wind  tunnel 
and  flight  test  data  [5.1],  is  used  to  generate  structured  uncertainty  models  of 
stability  and  control  derivatives.  These  models  are  then  translated  into  state-space 
element  uncertainty  models.  The  state-space  element  uncertainties  depend  on 
flight  condition.  However  for  simplicity,  the  uncertainty  model  for  each  state-space 
element, 


AA- 


±.02 1  Za  I 
±.04  I  Ma  I 
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±2lMql 


AB  = 


±.22 1  ZsE  I 
±.04  I  Msk  I  .  ’ 


(5.26) 
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is  held  constant  at  the  worst  flight  condition  case  thus  making  the  parameter 
uncertainty  models  conservative  in  a  sense.  The  additive  uncertainty  structure 
enters  the  state  equation  as 

X  =  (A  +  AA)x  +  (B  +  AB)u  .  (5.27) 

Recall  from  section  5.3.2  that  for  better  flying  qualities,  high  and  low 
dynamic  pressure  robust  performance  controllers  are  designed  and  implemented  by 
blending  both  controller  commands  with  the  blending  parameter  8b.  The  actual 
blending  parameter,  5b,  varies  between  0  and  1.  However,  it  is  normalized  such 

that  the  blending  parameter  for  analysis, 6b,  varies  betv/een  -1  and  1.  Therefore, 

analyzing  robustness  to  8b  is  equivalent  to  analyzing  robustness  for  all  outer-loop 
controller  command  combinations. 

Fig.  5.18  shows  the  robustness  analysis  model  including  the  structure  of  the 
uncertainty  as  well  as  the  normalization  weighting  elements.  To  show  the 
uncertainty  structure,  note  that  the  diagram  is  broken  down  to  the  scalar  loop  level 
except  for  the  plant  parameter  uncertainty  loops  v/hose  structure  is  given  by  oq. 
(5.27). 

Fig.  5.19  shows  the  robust  stability  plots  of  the  structure  shown  in  Fig.  5.18 
for  all  robustness  analysis  conditions  given  in  Appendix  5.  The  robust  stability 
plots  are  the  mixed  real/complex  structured  singular  values.  It  is  interesting  to  note 
the  absence  of  a  peak  n.ear  the  short  period  frequency.  This  suggests  insensitivity  of 
the  short  period  dynamics  to  the  given  uncertainty  structure.  The  peak  at 
approximately  40  rad/sec  corresponds  to  sensitivity  to  angle  of  attack 
measurements.  If  a  redesign  of  the  control  system  were  within  the  scope  of  this 
task,  inclusion  of  the  angle  of  attack  measurement  weight  in  the  design  model 
might  increase  robustness  to  angle  of  attack  measurements. 


6.3.4  FLYING  QUALITIES  ANALYSIS 


Several  flying  qualities  measures  of  section  2.2  are  used  to  analyze  the 
handling  characteristics  of  the  aircraft  and  the  designed  control  system.  The 
measures  for  pitch  response  are  given  in  terms  of  low-order  equivalent  systems 
(LOES)  that  represent  pitch  rate  (q)  and  normal  acceleration  at  the  instantaneous 
center  of  rotation  ni  responses  to  pilot  stick  deflection  inputs  (Sp).  Only  considering 
the  short-term  pitch  response  from  eq.  (2.10),  the  LOES  become 


i  t-  Te " 

Sp  s  +  2^gp(0^pS  +  (O^p 


+  2CspWspS  +  ^ 


sp 


(5.28) 


where  ^sp  is  equivalent  short-period  damping,  cOgp  is  equivalent  short-period 
frequency,  T02  corresponds  to  the  pitch  rate  zero  ,  and  xq  and  Tn  are  equivalent  pitch 
time  delays.  Only  the  parameters  of  the  stick  to  pitch  rate  transfer  function  are 
considered  for  analysis  since  once  T02  is  defined  and  it  is  assumed  the  time  delays 
are  equal,  the  transfer  functions  only  differ  by  the  steady-state  gain.  The  complete 
scheduled  flight  control  system,  full-order  actuator  models,  and  second-order 
aircraft  short-period  models  are  used  to  generate  high-order  closed-loop  aircraft 
linear  models.  LOES  are  generated  from  these  high-order  models  using  an 
equivalent  system  transfer  function  matching  program. 


The  Control  Anticipation  Parameter  (CAP)  defined  by  eq,  (2.18)  is  computed 
using  the  appropriate  LOES  parameters,  and  Fig.  5.20  shows  that  all  flight 
conditions  meet  Category  A,  Level  1  flying  qualities  requirements  for  CAP  and  aigp 
except  at  the  Mach  0.4,  6  kft  flight  condition  which  is  slightly  below  Level  1.  Figs. 
5.21  and  5.22  show  that  Level  1  requirements  are  met  for  cosp,  T02,  Csp)  and  tq. 

The  mismatch  defined  by  eq.  (2.15),  ranges  from  0.009  to  4.376  which 
suggests  a  close  LOES  approximation.  Section  2.2  describes  the  mismatch  as  well 
as  maximum  unnoticeable  added  dynamic  bounds.  Figs,  5.23  and  5.24  show  that 
the  stick  deflection  to  pitch  rate  added  high  order  dynamics  are  well  within  the 
maximum  allowable  unnoticeable  added  dynamics  represented  by  the  dashed  lines. 


Fig.  5.20  Control  Anticipation  Paramoter 


Fig.  5.21  Damping  and  Frequency 
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Fig  5.22  Equivalent  Time  Delay 
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Fig.  5.23  Umioticeable  Added  Dynamic  Gain 
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Fig.  5.24  Unnoticeablc  Added  Dynjunic  Phase 

The  mismatch  values  and  unnoticeable  dynamics,  Figs.  5.23  and  5.24,  show  that  the 
LOES  estimates  the  full-order  short-period  system  well  across  the  flight  envelope 
for  frequencies  of  interest  Since  the  phugoid,  or  long-term,  response  of  the  aircraft 
is  stable,  a  detailed  analysis  of  long-term  response  flvfTig  qualities  is  not  necessary. 
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The  nonlinear  time  responses  in  section  5.5  show  that  control  system  does  not 
destabilize  the  phugoid  mode. 


6.4  LATERAI/DIRECTIONAL  AXES  CONTROLLER 

The  lateral/directional  axes  manual  flight  control  system  generates 
differential  elevator,  aileron,  rudder,  and  roll  and  yaw  thrust  vector  commands  from 
sideslip  and  stability-axis  roll  rate  commands.  Those  control  effector  commands 
stabilize  the  roll  and  Dutch  roll  aircraft  dynamics  across  a  wide  conventional  flight 
envelope.  A  robust  oigenstructure  assignment  design  algorithm  is  used  for  the 
design  of  the  inner  loop  equalization  controller.  Structured  singular  value  synthesis 
is  used  to  design  outer  loop  robust  performance  controllers.  The  inner  and  outer 
loop  control  designs  are  presented  followed  by  robustness  and  Hying  qualities 
analyses  for  the  entire  lateral/directional  axes  flight  controller. 

6.4.1  INNER  LOOP  DESIGN 


The  desired  oigenstructure  for  the  inner  loop  is  chosen  to  consist  of  a  real  pole 
and  a  complex  pair  of  poles  that  correspond  to  the  desired  equivalent  roll  and  Dutch 
roll  poles  from  the  desired  flying  qualities.  The  desired  eigenvalues  are 

Xl<^=-1.67  and  = -2.1  ±  2.14  j  ,  (5.29) 

The  desired  eigenvectors  corresponding  to  the  desired  eigenvalues  are 
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where  x  represents  the  unspecified  clemeuts  in  the  eigenvectors.  Since  the  control 
effectiveness  matrix  is  replaced  by  the  generalized  control  effectiveness  matrix,  the 
plant  is  already  in  the  form  of  cq.  (2.54).  The  gain  matrix  from  the  three  outputs  to 
the  two  generalized  inputs  is  found  for  each  flight  condition,  i'he  gain  from  the  roll 
I’ate  to  the  generalized  yaw  acceleration  command  is  constrained  to  he  zero,  since 
the  gain  would  otherwise  be  in.sigiii.ficant.  This  clement  can  be  constrained  to  be 
zero  dii’ectly  in  the  equation.s  used  to  find  the  feedback  gain  matrix.  Tho^  feedback 
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matrix  eq.  (2.58)  is  rewritten  in  terms  of  a  Kronecker  product  and  a  row  stacking 
operator,  so  the  elements  of  the  feedback  gain  matrix  are  written  in  terms  of  a 
vector.  Each  row  is  solved  independently  for  each  of  the  feedback  matrix  elements, 
and  any  element  can  be  constrained  to  be  zero  by  removing  the  corresponding  row  of 
the  matrix  equation,  and  by  removing  that  element  from  the  vector  of  elements. 


Each  of  the  elements  in  the  inner  loop  gain  matrix  is  plotted  with  dynamic 
pressure,  and  a  polynomial  is  found  to  approximate  the  gain:  at  each  of  the  points. 
The  yaw  rate  to  roll  acceleration  and  yaw  acceleration  gains  are  also  functions  of 
altitude,  so  they  are  scheduled  as  a  function  of  static  pressure  ratio  as  well.  For  a 
particular  dynamic  piessure,  the  gain  is  found  at  discrete  altitudes  and  interpolated 
to  find  the  specific  gain  at  the  corresponding  altitude.  The  plots  of  the  inner  loop 
gains  appear  in  Fig.  5.25.  The  three  curves  in  Figs.  5.25c  and  5.25e  represent  the 
gains  at  three  different  altitudes. 
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Fig.  5.25  inner  Loop  Gains 


The  polynomial  gain  functions  of  Fig.  5.25  are  of  the  following  form; 


gain  -  C2  q^  +  ci  q  +  cq  , 


(5.31) 


where  the  coefficients  Cj,  i=:0,l,2,  are  given  in  Appendix  5. 
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Fig.  5.25  (cont.)  Inner  Loop  Giiins 


5,4.2  OUraR  LOOP  DESIGN 


-Since  the  inner  loop  I'cedliack  gains  and  the  control  selector  aie  both 
scheduled  v/ith  flight  condition,  the  ijsncr  loop  dynamics  are  invariant  with  respect 
to  flight  condition.  For  the  outer  loop  feedback  design,  a  nominal  flight  condition  is 
chosen,  so  nominal  plant  dynamics  and  the  corresponding  inner  loop  feedback  gains 
can  be  used  in  the  synthesis  model.  The  flight  condition  at  Mach  0.5  and  20,000 
feet  altitude  is  chosen  as  the  nominal  design  point.  The  dynamic  prcssurii  at  this 
condition  is  170  psf,  and  the  trim  angle  of  attack  is  5.1  degrees.  The  control 
effectiveness  matrix  is  replaced  by  the  generalized  contn;!  eflectiveness  given  in  eq. 
(5.5). 
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The  synthesis  model  used  for  the  outer  loop  control  law  design  is  shown  in 
Fig.  5.26. 


Fig.  5.26  Outer  Loop  Design  Model 


Included  in  the  model  are  the  nominal  plant  model  including  the  control  selector, 
the  irmer  loop  feedback  gains,  Keq,  and  the  output  transformation,  Tl,  given  by 


-10  0  - 
.  0  cos  a  sin  a  J  ’ 


(5.32) 


Performance  is  achieved  by  including  an  ideal  model  of  the  roll  rate  and  sideslip 
angle  responses  in  the  synthesis  model.  By  using  target  flying  qualities 
parameters,  the  ideal  model  is; 


_ _ _ _ 9 _ 

Pc  ~  (s2  +  4.2s  +  9)  ’ 


Hidsi-’ 


(0.6s  +  1) 
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Uc 


(5.33) 


To  achieve  performance,  an  error  signal  is  generated  between  the  ideal  model 
and  the  actual  model  responses.  This  error  is  weighted  with  frequency  dependent 
weights,  shown  in  Fig.  5.27,  and  given  by 


Wp  =  diag(WpP  Wp^^), 

w  P  -  0-Q5(s  +  500)  w  -  0-QQ5(s  +  500) 
~  (s  +  0.1)  -  (s  +  0.1) 


Fig.  5.27 


Performance  Weights 


(5.34^ 


The  input  matrix  is  replaced  by  the  generalized  input  matrix  in  the  nominal 
design  model  to  represent  the  control  selector  matrix.  Hecausc  of  this,  the  real 

o  »•/%  #*»*»♦->  1 »  r*» »  ^  »  *-v  » I- <-»  f*  #»»» 
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yaw  acceleration.  Since  the  real  actuator  models  cannot  be  included  in  this 
arrangement,  the  actuator  models  used  in  the  design  are  fictitious,  and  are  given  by 


6_ _ 20 

~  (s  +  20) 
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To  prevent  unrealistic  control  commands,  the  actuator  rates  are  weighted  in  the 
problem  formulation  by 


Wact  =  0.005  h , 


(5.36) 


Robustness  to  unstructured  uncertainty  is  achieved  by  including  a 
multiplicative  uncertainty  at  the  plant  input.  The  uncertainty  boimds  are  shown  in 
Fig,  5.28,  and  are  given  by 


Wa  =  diag(W{,.Wv), 


10(s  +  5) 

(s  +  50,000)  • 


(6.37) 


Fig.  5.28  Uncertainty  Bounds 


The  synthesis  model  is  rearranged  to  form  the  interconnection  structure  in 
Fig.  5.29.  Here,  zi  is  the  weighted  error  between  the  ideal  model  and  the  output  of 
the  plant,  z2  is  the  weighted  actuator  rate  vector,  and  zg  is  the  output  from  the 
unstructured  uncertainty  scaling  Wa-  The  perturbation  block  includes  a  2  by  4 
performance  block  between  r  and  zj  and  z2  and  a  2  by  2  uncertainty  block  between 
d  and  z3. 
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Fig.  5.29  Lateral  Interconnection  Structure 
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jj-synthesis  is  applied  to  the  interconnection  structure  by  first  designing  an  H« 
control  law.  The  initial  bound  on  the  Hoo  norm  is  1.325.  |.i-analysis  is  performed  on 
the  closed  loop  system.,  and  third  order  D-scaling  matrices  are  fit  to  the  frequency 
dependent  D  scales  used  to  find  the  upper  bound  of  )i.  These  are  used  to  scale  the 
original  plant.  An  Hoo  control  law  was  found  for  the  scaled  plant,  with  an  H«  bound 
of  0.552.  (i-analysis  was  performed  again,  and  third  order  D-scales  were  used  a 
second  time.  The  final  control  law  was  found  with  an  Hw  nonn  bound  of  0.545. 


The  resulting  |i-synthesis  controller  has  24  state  variables  and  is  stable.  The 
controller  is  balanced  and  reduced  to  10  state  variables  by  truncating  the  14  state 
variables  with  the  smallest  Hankel  singular  values.  The  resulting  lUth  order 
controller  includes  2  large  negative  polos,  one  of  which  is  -  ^OOO.O,  the  other  of  which 
is  -472.8.  They  are  residualizod  to  form  the  final  8th  order  controller.  The  full- 
order  and  reduced-order  controllers  arc  given  in  Appendix  5. 

6.4.3  ROBUSTNESS  ANALYSIS 


Although  the  control  law.s  are  designed  with  a  specific  accounting  for 
unstructured  uncertainty  at  the  plant  input  to  account  for  the  uncertainty  in  the 
actuators,  there  are  many  sources  of  uncertainty  in  the  dynamic  aircraft  model  that 
were  not  accounted  for  directly  in  the;  design.  3’o  try  t  o  account  directly  for  ail  of  the 
sources  of  uncertainty  would  be  computationally  intensive  and  re.sult  in  a  very  large 
order  controller.  Robustnes.s  must  he  verified  against  all  sources  of  uncertainty 
after  the  design  is  completed  to  a.sse.ss  the  need  for  rede.sign  with  more  rigorou.s 
treatment  of  different  types  of  uncertainty. 

While  fictitious  control  inputs,  actuator  models,  and  uncertainty  descriptions 
are  used  in  the  design  process  to  be  able  to  incorporate  the  control  selector  into  the 
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nominal  plant  model,  the  actual  control  input  matrix,  dynamic  actuator  models,  and 
uncertainty  models  of  the  actuators  are  used  in  the  analysis  of  robustness.  The 
weights  used  for  the  unstructured  uncertainty  at  the  plant  input  are  plotted  in  Fig. 
5.30,  and  are: 

Wb  =  Wdt  =  l~3Ty  Wa  =  ,  (5.38) 


Fig.  6.30  Weights  for  Actuator  Uncertainty 


Note  that  only  the  aerodynamic  actuators  are  used  in  the  analysis  for  the  same 
reasons  discussed  in  Section  5.3.3. 

In  addition  to  unstructured  uncertainty  at  the  input  to  reflect  the  modeling 
errors  and  neglected  dynamics  in  the  actuators,  there  is  also  corresponding 
unstructured  uncertainty  at  the  output  of  the  plant  to  reflect  the  modeling  errors 
and  neglected  dynamics  in  the  sensors.  The  measurements  assumed  for  this  design 
are  the  body  axis  roll  and  yaw  rate  signals,  and  the  sideslip  angle  signal.  Although 
it  is  now  practical  to  generate  a  sideslip  angle  signal  for  feedback  control,  the  signal 
is  often  corrupted  with  noise  and  often  has  significant  errors.  The  robustness  of  the 
closed  loop  system  to  errors  at  the  plant  output  is  also  tested.  The  weights  for  the 
unstructured  uncertainty  at  the  plant  output  are  plotted  in  Fig.  5.31,  and  arc: 


(6.39) 


Wp  =  Wr  =  ^,-200). 


“  (s  +  1000)- 


_ _a. _ L.,i.a  1.1.11 
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CO  (rad/sec) 


Fig.  5.31  Weights  for  Sensor  Uncertainty 
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Finally,  the  stability  and  control  derivatives  in  the  linear  state  space  model 
are  uncertain  to  various  degree.s.  Errors  in  modeling  and  wind  tunnel  testing,  as 
well  as  variations  in  moments  of  inertia  and  mass  with  fuel  burn  all  contribute  to 
the  uncertainty  in  the  stability  and  control  derivatives.  The  level  of  uncertainty 
assumed  for  each  stability  or  control  derivative  is  shown  in  Table  5.1,  and  is  baaed 
on  confidence  in  each  of  the  derivatives.  Each  uncertainty  i.s  given  as  a  percentage 
of  the  nominal  value. 

The  interconnection  structure  for  the  robustness  analysis  consists  of  the 
eighth  order  outer  loop  regulator,  the  inner  loop  feedback  gain  matrix,  the  control 
selector  matrix,  the  actuator  models,  including  a  fourth  order  linear  model  to 
represent  the  elevator  and  second  order  linear  models  to  represent  the  ailerons  and 
the  rudder,  and  the  third  order  lateral/directional  model.  Tlic;  third  order  rather 
than  the  fourth  order  iateral/directional  model  is  used  because  of  the  marginally 
stable  spiral  mode.  If  the  structured  singular  value  were  evaluated  with  the 


Table  5.1  Laterai/Directional  Parameter  Uncertainty  Bounds 


Stability  Derivatives 

Control  Derivatives 

Derivative 

Derivative 

HIHBHlHi 

15 

Y6r 

15 

10 

b5DT 

15 

Lo 

30 

LSa 

10 

Lr 

20 

b6R 

40 

Np 

30 

Nsdt 

15 

No 

50 

N6a 

20 

Nr 

15 

N5r 

15 

margi.nally  stable  spiral  mode  included  in  the  dynamics,  it  would  ho  very  large 

since  very  small  perturbations  would  make  this  mode  unstable.  Tins  would  corrupt 
the  analysis  since  an  unstable  spiral  mode  can  be  acceptable  with  respect  to  Level  1 
flying  qualities. 


The  structured  singular  value  is  assessed  at  four  different  flight  conditions  te 
detex'mine  robust  stability.  The  flight  conditions  used  for  robustness  analysis 
appear  in  Appendix  5.  The  perturbation  matrix  included  a  complex  3  by  3  matrix 
for  the  unstructured  uncertainty  at  the  input,  a  complex  3  by  3  matrix  for  the 
unstructured  uncertainty  at  the  plant  output,  and  14  scalar  uucei  tainties  for  the 
uncertain  stability  and  control  derivatives.  First,  the  scalar  perturbations  to  the 
stability  and  control  derivatives  are  assumed  to  be  complex.  Since  the  parameter 
values  are  always  real  in  the  actual  plant,  this  analysis  is  conservative.  Next,  the 
scalar  perturbations  are  constrained  to  be  real.  This  analysis  i.s  much  loss 
conservative. 

The  structured  singular  values  are  plotted  in  Fig.  5.32  for  each  of  the  four 
analysis  conditions.  Where  the  parameter  perturbations  are  assumed  to  be 
complex,  the  plots  are  shown  as  solid  lines,  and  for  parameter  perturbations 
assumed  to  be  real,  the  plots  are  shown  as  dashed  lines. 


a) 


b.  Mach  0.5,  20  kfi 


CO 


c.  Mach  0.4,  30  left 


d.  Much  0.9,  30  kft 


Fig,  5.32  liobu.st  Stability  Ploty 


At  each  flight  condition,  the  structured  singular  value  sliows  a  large  peak  near  3 
rad/sec  when  the  parameter  uncertainty  is  assumed  to  be  complex.  This  peak 
corresponds  to  some  perturbation  that  will  make  the  closed  loop  Dutch  roll  poles 
near  that  frequency  unstable.  The  peaks  range  from  1.3  to  2.0,  meaning  that  the 
true  plant  can  only  tolerate  between  1/1.3  and  1/2.0  of  the  allowed  perturbation. 

V/hen  the  structured  singular  value  is  computed  with  the  parameter 
perturbations  constrained  to  be  real,  no  peaks  appeared  izi  the  structured  singular 
value  near  3  rad/sec,  meaning  that  the  parameter  perturbation  that  the  plant  is 
sensitive  to  is  complex,  a  situation  which  can  never  physically  occur.  The  maximum 
values  of  the  structured  singular  values  now  lie  between  0.65  and  1.3.  For  two  of 
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the  cases  shown,  the  value  of  n  does  exceed  unity,  meaning  that  for  some 
combination  of  parameter  perturbations,  as  well  as  perturbations  in  the  actuators 
and  sensors,  the  plant  will  go  unstable.  Although  robustness  for  all  of  the 
simultaneous  uncertainty  is  not  guaranteed,  the  analysis  shows  that  the  closed  loop 
system  is  not  extremely  sensitive  to  any  of  the  uncertainty,  as  many  types  of 
uncertainty  were  included  in  the  analysis,  and  the  value  of  p  only  exceeded  one 
slightly. 

5.4.4  FLYING  QUALITIES  ANALYSIS 

The  flying  qualities  at  several  flight  conditions  are  tested  by  fitting  the  actual 
high  order  response  to  a  low  order  equivalent  system  (LOES)  and  evaluating  the 
parameters  of  the  LOES.  The  flight  conditions  tested  are  listed  in  Appendix  5.  The 
flight  conditions  are  chosen  to  be  a  subset  of  those  used  for  design,  and  are  chosen  U) 
represent  the  entire  range  of  dynamic  pressure.  The  response  of  interest  is  the 
fourth  order  equivalent  system  from  the  stability  axis  roll  rate  command  to  the 
stability  axis  roll  rate  given  by  eq.  (2.19).  The  high  order  response  consists  of  the 
fourth  order  linear  plant  model,  the  actuator  dynamics,  the  control  selector  matnx, 
the  inner  loop  gain  matrix,  and  the  eighth  order  outer  loop  controller,  Only  the 
aerodynamic  control  effectors  are  used  in  the  analysis. 

To  get  a  good  equivalent  system  fit,  the  high  order  roll  rate  response  is  first 
fit  to  a  first  order  LOES  to  get  the  equivalent  roll  mode  time  constant.  Next,  the 
high  order  sideslip  angle  response  is  fit  to  a  second  order  L(3ES  to  jjet  the 
equivalent  Dutch  roll  frequency  and  damping.  The  fourth  order  roll  rate  LOES  is 
found  by  fixing  these  parameters  and  using  the  remaining  free  parameters  in  eq. 
(2.19)  to  find  the  fourth  order  equivalent  system. 

Fig.  5.33  shows  the  equivalent  Dutch  roll  damping  and  frequency  for  each  of 
the  flight  conditions  tested.  For  each  point,  the  parameters  are  well  within  the 
Level  1  boundaries.  In  addition,  all  of  the  points  tested  are  very  close  together, 
although  the  points  tested  represent  a  wide  range  of  dynamic  pressure.  This  is 
because  of  the  gain  scheduled  inner  loop.  Figs.  5.34  and  5.35  show  the  equivalent 
roll  mode  time  constant  and  eqviivalent  time  delay,  respectively,  each  plotted  as  a 
function  of  dynamic  pressure. 
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boundaries,  and  each  is  fairly  constant  with  dynamic  pressure  because  of  the  gain 
scheduled  inner  loop. 


5.5  NONLINEAR  ANALYSIS 


Using  the  nonlinear  simulation  with  the  control  laws  implemented  as 
FORTRAN  subroutines,  several  maneuvers  are  simulated.  First  the  outer  loop 
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controller  blending,  described  in  section  5.3.2,  is  analyzed  with  angle  of  attack  unit 
step  commands  given  to  the  longitudinal  contiol  system.  Next,  three  highly  coupled 
maneuvers  are  performed  and  analyzed  to  evaluate  the  control  system  during 
realistic  operating  conditions. 

Angie  of  attack  step  inputs  are  used  to  demonstrate  outer  loop  controller 
blending  performance  of  the  longitudinal  control  system  during  nonlinear 
simulations  at  high,  intermediate,  and  low  dynamic  pressure  conditions.  The  three 
different  conditions  are  chosen  such  that  three  different  combinations  of  outer  loop 
controller  commands  are  used  based  on  the  structure  shown  i.n  Fig.  5.16  and  Fig. 
5.17.  Fig.  5.36  shows  the  nonlinear  responses  of  the  aircraft  model  to  identical 
angle  of  attack  unit  step  inputs  at  the  tliree  different  conditions.  The  high  and  low 
dynamic  pressure  conditions  use  only  the  high  and  low  dynamic  pressure 
controllers,  respectively.  However,  the  middle  dynamic  pressure  condition  uses  a 
combination  of  both  controllers  resulting  in  an  angle  of  attack  response  between  the 
high  and  low  dynamic  pressure  conditions.  Note  that  thrust  vectoring  is  not 
required  for  this  benign  maneuver. 

Next,  a  3g  loaded  roll  is  performed  at  Mach  0.3  and  20,000  feet.  This 
maneuver  is  a  coupled  maneuver  that  excites  dynamic  modes  in  the 
lateral/directional  and  longitudinal  axes.  The  maneuver  consists  of  the  aircraft 
entering  a  3g  turn,  with  a  roll  reversal  to  a  3g  turn  in  the  opposite  direction  without 
unloading  the  aircraft.  The  maneuver  is  shown  in  Fig.  5.37.  The  load  factor  is 
maintained  near  3g's  throughout  the  maneuver,  while  the  sideslip  angle  never 
exceeds  1.5  degrees.  This  maneuver  does  not  saturate  the  aerodjmamic  surfaces, 
therefore  thrust  vectoring  is  not  commanded. 

A  high  rate  roll  is  performed  at  the  flight  condition  of  Mach  0.5  and  20,000 
feet.  A  250  deg/sec  roll  rate  is  commanded  for  a  360  deg.  roll.  The  state  responses 
and  control  deflections  are  plotted  in  Fig.  5.38.  The  stability  axis  roll  rate  response 
is  satisfactory  with  good  turn  coordination  in  that  sideslip  remains  less  than  2  deg. 
Thrust  vector  commands  are  quickly  generated  since  the  ailerons  and  rudders 
immediately  rate  saturate  to  the  point  of  position  saturation.  At  about  3  sec  into 
the  maneuver,  aerodynamic  surface  rate  saturation  ceases,  and  thrust  vector 
commands  are  not  needed. 
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Fig.  5.37  Loaded  Roll  Maneuver  'Hme  Response 
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Fig.  5.38  (cont.)  High  Rate  Roll  Maneuver  Time  Response 

A  turn  reversal  is  also  performed  at  Mach  0.6  and  20,000  feet.  This 
maneuver  consists  of  a  bank  to  the  right,  a  baidc  to  the  left,  another  bank  to  the 
right,  and  a  bank  back  to  steady,  level  flight.  The  turn  reversal  maneuver  is  plotted 
in  Fig.  5.39.  The  roll  rate  responds  as  desired,  and  the  sideslip  angle  is  maintained 
at  less  than  4  deg  throughout  the  maneuver.  This  maneuver  is  also  extreme  enough 
to  warrant  thrust  vectoring  about  the  lateral/dircctional  axes  because  of 
aerodynamic  surface  saturation. 


Fig.  5.39  Turn  Reversal  Maneuver  Time  Re.spoiase 
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Fig.  5.39  (cent.)  Turn  Reversal  Maneuver  Time  Response 


Note  that  stability  ar  d  performance  of  the  aircraft  are  maintained  even  though  the 
high  rate  roll  and  turn  reversal  maneuvers  are  extreme  enough  to  saturate  the 
thrust  vectoring  nozzles. 
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5.6  CONCLUSIONS  AND  LESSONS  LEARNED 


A  controller  structure  is  presented  that  allows  the  separation  of  design  goals 
during  the  synthesis  process.  An  inner  loop  controller  is  used  to  equalixe  plant 
dynamics  across  the  flight  envelope.  An  outer  loop  controller  is  used  to  provide 
robust  flying  qualities.  A  control  selector  allows  the  inner  and  outer  loop  controllers 
to  command  body  axes  rotational  accelerations  thus  normalizing  the  control 
effectiveness  across  the  flight  envelope.  Daisy-chaining  is  used  within  the  control 
selector  to  select  thrust  vector  control  only  when  aerodynamic  control  is  not 
sufficient. 

It  is  found  that  the  prioritization  of  redundant  conti'ol  effectors  within  the 
control  selector  is  important.  Since  the  elevators  are  the  only  aerodynamic  pitch 
control,  the  differential  elevator  commands  must  be  reduced  to  allow  for  symmetric 
elevator  commands  if  needed.  The  choice  of  the  inner  loop  equalized  dynamics 
turns  out  to  be  crucial.  The  equalized  dynamics  should  bo  chosen  similar  to  the 
high  dynamic  pressure  models  to  avoid  deaugmentation  of  the  aircraft  at  high 
dynamic  pressure  conditions.  The  control  system  seems  to  make  the  aircraft 
sensitive  to  angle  of  attack  measurements.  This  problem  could  possibly  be  reduced 
by  including  angle  of  attack  measurement  dynamics  in  the  outer  loop  design  model. 
It  is  also  noted  that  the  methods  used  for  inner  and  outer  loop  control  design  give 
high  gain  controllers  for  the  longitudinal  axis.  It  is  suspected  that  additional 
feedforward  compensation  and  tuning  of  the  design  weights  could  have  solved  this 
problem. 
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APPENDIX  6 

This  section  gives  the  model  data  of  the  design  presented  in  the  previous 
sections.  First,  detailed  actuator  models  are  given,  then  linear  models  and 
information  about  the  flight  conditions  at  which  the  linear  models  are  generated  are 
given.  Inner  loop  control  gain  schedules  are  given,  and  finally,  full  and  reduced 
order  oute.r  loop  dynamical  controllers  are  presented. 

ACTUATOR  MODELS 

High-order  detailed  aerodynamic  actuator  models  are  reduced  to  simple  4th" 
order  and  2nd-ordef  models,  shown  in  Table  A5.1,  to  avoid  very  small  time-step 
integration  during  nonlinear  simulation. 


Table  A5.1  Actuator  Models 


Surface 

Actuator  Dynamics 

Rate 

Limits 

(deg/sec) 

Position 

Limits 

(deg) 

Elevator 

,  s  ,  2  2(0.0(38)  ,  ,, 

^82.9^  ^  82.9 

,,  s  ,2  2(0.41)  ,  ,,  s  -2 .  2(0.59) ,  , 

'^36.4^  3(3.4  *  ’  105.3^  105.3 

-60,0 

460.0 

-24.0 

410.5 

Leading 
Edge  Flap 

1.0 

'2(3.9  '■  '82.9  ^ 

-40.0 

4  40.0 

0.0 

435.0 

Trailing 
Edge  Flap 

LO 

s  2  2(0.71) 

Crr/-  +  ~  -iF  n  ^ 
o5.0  3.).0 

•100.0 

4 !  00,0 

-10.0 

445.0 

Aileron 

1.0 

,"r  ;2  '  2(0:59)  “ 

"  75.0 

-100,0 

4 100.0 

-25.0 

445.0 

Rudder 

1.0 

—  -  - 

720^  72,0 

-100,0 

4l00.() 

-30.0 

430.0 

Vectoring 

Nozzles 

1.0 

s  ’  2 ''  2(0.(30)'' 

20.()  20.0  ^  ^ 

-60.0 

460.0 

-30.0 

430.0 

I 
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Thrust- vectoring  nozzle  actuators  are  assumed  to  be  2nd-order  filters  and  are  also 
given  in  Table  A5.1.  These  reduced-order  models  are  used  in  the  nonlinear 
simulation  and  linear  robustness  analysis. 


LINEAR  DESIGN  MODELS  AND  I^’L  GHT  CONDITIONS 

The  flight  conditions  that  are  chosen  for  linear  control  design,  flying  qualities 
analysis,  and  robustness  analysis  a»-e  given  in  Table  A5.2.  The  conditions  chosen 
represent  a  broad  range  of  dynamic  pressures  throughout  the  envelope  since  inner- 
loop  gains  are  scheduled  versus  dynamic  pressure. 

Table  A5.2  Flight  Conditions 


a.  Longitudinal  Design  Conditions 


Mach 

Number 

Altitude 

(ft) 

q 

(psf) 

a 

(deg) 

Flying 

Qualities 

Analysis 

Robustness 

Analysis 

0.3 

26000 

47.4 

25.2 

V 

0.5 

40000 

68.5 

16.8 

V 

V 

0.6 

30000 

158.4 

5.2 

V 

0.4 

6000 

189.9 

6.0 

V 

V 

0.7 

14000 

426.4 

2.6 

V 

V 

0.8 

12000 

603.0 

1.9 

V 

V 

0.95 

20000 

614.4 

1.6 

V 

0.8 

10000 

652.0 

_  .  1,7 

V 

0.8 

5000 

789.1 

1.5 

V 

0.9 

10000 

825.2 

1,4 

V 

0.85 

5000 

890.8 

1.4 

V 

0.9 

5000 

998.7 

1.3 

V 

V 

Table  A5. 2  (cont.) 


b.  Lateral/Direc 

Lional  Design  Condition.s 

Mach 

Number 

Altitude 

(ft) 

9 

(psD 

a 

(deg) 

Flying 

Qualities 

Analysis 

Robustness 

Analysis 

0.2 

40.75 

29.7 

0.3 

91.69 

10.0 

0.4 

BWHI 

163.00 

5.6 

V 

0.5 

BBHH 

254.68 

3.6 

0.7 

BEISIffllli 

499.18 

1.8 

V 

0.9 

825.17 

1.2 

v 

V 

0.3 

BMWi 

61.27 

12.6 

< 

0.4 

20,000 

108.92 

7.8 

0.5 

20,000 

170.19 

5.4 

V 

V 

0.6 

245.07 

3.9 

0.75 

382.92 

2.8 

0.9 

551.40 

2.1 

n' 

0.4 

30,000 

70.39 

17.4 

\I 

V 

0.5 

30,000 

109.98 

8.1 

0.6 

30,000 

158.37 

5.2 

V 

0.7 

fc!8R!8M 

215.55 

3.6 

0.8 

■iMiM 

281.54 

2.2 

V 

1.5 

V 

The  linear  models  for  each  flight  condition  of  Table  Ab.2  are  listed  below,  whore  the 
nomenclature,  for  example  ,  is  the  longitudinal  state  matrix  at  Mach  .9  and 

lOkft. 

The  longitudinal  linear  models  have  the  following  form: 


r  •  ^ 
a 

r  1 

r  (X ' 

6e 

-  q  - 

M  M 

<x  q  j 

1 

L  q  J 

+ 

-  M5E  Mgpi-y  . 

.5p'"v. 

=  Ai 


ong 


l^long 


bE 

.5fi'V- 


(A5.1) 
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.  m3h26 

■  -0.2296 

0.9931 

1  r.ni3h26 

J  “long 

-0.04034 

-0.01145 

^long 

” 

.  0.02436 

-0.2046 

-1.73 

-0.517 

.  m5h40 

r  -0.2423 

0.9964  1  „m5h40 

-0.1737  J  ^Oong 

-0.0416 

-0.01141  ■ 

•^long 

"1 

-2.342 

-2.595 

-0.8161 

.  m6h30 

■  -0.5088 

0.994  ■ 

pinGhSO 

“long 

r 

-0.09277 

-0.01787  ■ 

■^long 

-1.131 

-0.2804  . 

L 

-6,573 

-1.525 

*  m4h6 

■  -0.8018 

0.9847  ■ 

pm4ti6 

-0.1508 

-0.02776  ■ 

^long 

-1.521 

-0.5944  . 

“long 

-7.926 

-1.751 

.  m7hl4 

■  -1.175 

0.9871  1 

pm7hl4 

“long 

-0.194 

-0.03593  1 

^long 

,  -8.458 

-0.8776  J 

!_ 

-19.29 

-3.803  J 

•  m8hl2 
^long 

■  -1.562 

0.9862  1 

pm8hl2 

“long 

“ 

-0.2316 

-0.04349  1 

•M 

-14.94 

-1.132  J 

. 

-26.48 

-5.323  J 

*  m95h20 

■  -1.905 

0.9895  ] 

pni95h20 

“long 

* 

-0.1867 

-0.03287  1 

^long 

““ 

-33.88 

-0.9872  J 

“ 

. 

-27.22 

-4.573  J 

.  inShlO 
“long 

■  -1.675 

0.9853  ] 

pinShlO 

“long 

r 

-0.2449 

-0.04649  1 

-16.16 

-1.212  J 

““ 

L 

-28.34 

-5.742  J 

*  m8h5 
“long 

■  -1.994 

0.9828  I 

pm8h5 

“long 

“ 

-0.2852 

-0.05567  ] 

““ 

-19.44 

-1.427  J 

-33.44 

-6.931  J 

.  ra9hl0 

■  -2.452 

0.9856  ] 

pHiOhlO 

“long 

* 

-0.2757 

-0.05226  1 

“long 

-38.61 

-1.34  J 

-37.36 

-7.247  J 

» m85h5 
“long 

■  -2.328 

0.9831  ] 

pm85h5 

^long 

- 

-0.3012 

-0.05866  ■ 

-30.44 

-1.493  J 

-38.43 

-7.815 

*  m9h5 

■  -2.911 

0.9835  1 

“long 

r 

1 

-0.3161 

-0.06231  ■ 

“long 

-46.47 

-1.553  J 

s 

L 

-43.65 

-8.752 
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The  iaterai/directional  linear  models  have  the  following  form: 


Yp 

sin  a 

-cos  (X 

P 

Lp 

Lp 

Lr 

Lnp 

Np 

Nr  _ 

r 

5dt 


^8DT  ^  5A  ^SR'l’V  "‘^SY’l'V  ' 
^SDT  Li5a  L'5R  ^^5RTV  ^■‘SY’rV 
.  ^5DT  ^5A  ^5R  ^5RTV  ^6YTV . 


5a 

^^R 

6rtv 


LSyi'vJ 


~  Ajat/dir 


I^lat/dir 


8ot 

5a 

5r 

6rtv 


L  <‘»YT  j 


.  m2hl0 
■^lat/dir 


Din2hlO 

“lat/dir 


.  mShlO 
^lat/dir 


rsmShlO 

“lat/dir 


.  m4hlO 
^lat/dir 


T3m4hlO 

at/dir 


.  m5hl0 
■^lat/dir 


QinDhlO 

^lat/dir 


■  -0.05904  0.4959  -0.8703  1 

-5.513  -0.9391  0.6655 

.  0.06838  0.02632  -0.1038  . 

■  0.005629  0.005764  0.003685  0  0,0904  1 

1.879  1.328  0.02922  0.6754  0.217 

.  -0.1092  -0.09645  -0.08404  0.006811  -2.974  . 

‘  -0.1292  0.1738  -0.9833  I 

-8.643  -1.129  0.5986 

1.519  -0.01327  -0.1105 


'  -0.006987  -0,005249  0.01285  0 

5.096  6.075  0.51  0.1781 

.  0.1908  -0.1522  -0.3872  0.0008849 


0.006894 

0.02478 

-0.3397 


-0.1544  0.09691  -0.9939  ' 
-9.965  -1.721  0.599 

2.169  -0.01995  -0.1447 


-0.01187  -0.006276  0.01785  0  0.002272  ‘ 

9.643  12.16  0.9326  0.1495  0.01088 

.  0.2768  -0.2727  -0.7155  -0.000305  -0.1492 

'  -0.1932  0.06234  -0.9968  1 
-12.37  -2.164  0.6034 

3.199  -0.0211  -0.1802 


-0.01652 

15.18 

0.3458 


-0.007459 

18.04 

-0.3975 


0.02203 

1.412 

-1.099 


0 

0.3035 

-0.002071 


0.002821  ■ 
0.01689 
-0.2315 


A 


inVhlO 
lat/dir  " 


B 


B 


B 


B 


.  in6h20 
^lat/dir 


B 


m7hlO 
lat/dir  ■" 


m9hlO 
Mat/dir  “ 


m9hlO 
lat/dir  " 


m3h20 
Nat/dir  “ 


Qm3h20 

^lat/dir 


m4h20 
Mat/dif  “ 


m4h20 
lat/dir  = 

in5h20 
Mat/dir  “ 

m5h20  _ 

lal/dir  “ 


m6h20  _ 

lal/dir  “ 

m75h20 
^at/dir  “ 


•rpm75h20 

“lat/dir 


.  m9h20  _ 

■“lat/dir  " 


■  -0.2701 

0,03162  -0.9984  1 

* 

-17.77 

-3.177  0.5446 

5.987 

-0.0205  -0.2555  . 

■  -0.02472 

-0.00764  0.0285r 

> 

0 

0.01074 

29.4 

25.95  2.503 

1.739 

0.09002 

0,4006 

-0.3672  -1.982 

-0.02069  -1.234 

■  -0.321 

0.02008  -0.9987  ^ 

1 

= 

-17.6 

-5.716  0.5193 

9,433 

-0.02149  0.3391  . 

“  -0.0301 

0  0.03051 

0 

0.02583  ■ 

= 

46.68 

19.41  4.054 

5.024 

0.2782 

0.2385 

-0.3226  -2.951 

-0.06746 

-3.813 

■  -0.0714 

0.2991  -0.9539 

1 

=: 

-6.746 

-0.5918  0.4968 

0.4099 

0.003356  -0.08297  _ 

-0.00274  -0.002723  0.006669 
2.784  2.627  0.2294 

0.04708  -0.1018  -0.2059 

-0.112  0.1408  -0.9889 

-8.538  -1.171  0.5146 

1.619  -0.01304  -0.103 

-0.007388  -0.004613  0.01238 
6.245  7.673  0.611 

0.1998  -0.1822  -0.4722 

-0.1354  0.09036  -0.9949 
-10.37  -1.469  0.5126 

2.281  -0.01482  -0.1277 

-0.01091  -0.005695  0.01555 
9.93  12.12  0.9416 

0.2757  -0.2797  -0.7419 

-0.166  0.0629  -0.9971 

-12.97  -1.761  0.5083 

3.191  -0,01417  -0.1529 

-0.0142  -0,00686  0.01851 
14.38  16.76  1.316 

0.3389  -0.385  -1.051 

-0.1982  0.03905  -0.9984 
-17.3  -2.505  0.4624 

4.688  -0.01064  -0.1942 

-0.01816  -0.007716  0.02168 
22.53  23.71  1.905 

0.4211  -0.238  -1.555 

-0.2257  0.02638  -0.9989 
-19.08  -3.708  0.4264 

6.586  -0.01925  -0.2379 


0 

0.5298 

0.005966 


0,03367 
0.1 165 
-1.597 


0 

0,3615 

0,0008652 


0.008751 

0.04034 

-0.553 


(,) 

0.3977 

-0.001175 


0.00489 

0.02817 

-0.3861 


0 

0.7007 

0.00475 


0.005817 

0,0402 

-0.5511 


0 

1.697 

-0.01765 


0.01018 

0.08794 

-1.206 
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■rjm9h20 

“lat^dir 


»  m4h30 
■^lat/dir 


I^ra4h30 

“lat/dir 


.  m5h30 
■^lat/dir 


T3ni5h30 
“1  at/dir 


.  m6h30 
■^lat/dir 


nin6h30 

^lajydir 


.  m7h30 
■^1  at/dir 


Qm7h30 

“lat/dir 


*m8h30 

■^'iat/dir 


„m8h30 

“lat/dir 


.  m9h30 
■^1  at/dir 


^m9h30 

^lat/dir 


0 

0.606 

0.004285 


-0.02089  -0.005598  0.02814  0 

31.29  25.88  2.879  8.359 

0.3459  -0.1457  -2.152  -0.04183 

-0.07258  0.2183  -0.9751 

-6.216  -0.6013  0.427 

0.9358  -0.005054  -0.07243 

-0.003475  -0.008869  0.007457 
3.491  3.867  0.3446 

0.09003  -0.115  -0.2766 

-0.09336  0.1351  -0.99 

-8.39  -0.9786  0.487 

1.629  -0.009743  -0.0886 

-0.006711  -0.003948  0.01047 

6.183  7.5  0.5958 

0.2107  -0.1819  -0.475 

-0.1118  0.09868  -0.9949 

-10.22  -1.169  0.4318 

0.1056  _ 

0.01262 
11.06  0.842 

-0.2637  0.6841 

0.06791  -0.997 


0.01809 

0.1875 

-2.571 


0.02451 

0.1084 

-1.486 


2.201  -0.009853 

-0.009231  -0.004885 
9.02 
0.2957 
-0.1275 


0 

0.6251 

0.001079 


0 

0.7304 

-0.001928 


0.01242 

0.06865 

-0.9411 


0.008328 

0.05522 

-0.757 


-11.92 

2.776 

-0.01108 

12.54 

0.3656 

-0.1397 

■14.78 

3.437 

-0.01269 

16.56 

0.3773 

-0.1542 

-16.19 

4.42 

-0.01401 

20.34 

0.3122 


-1.404  0.4283 

•0.01099  -0.124 

-0.005726  0.01428 
14.81 


-0.2088 

0.04954 

-1,859 

-0.003867 

-0.006594 

18.24 

-0.0H701 

0.03692 


1.085 
-0.9165 
-0.9981  ■ 
0.3924 
-0.1429  . 

0.0 1585 

i.:r/2 

0.998V 


0 

0.9813 

-0,005925 


0.007537 

U.0583 

•0.7992 


-2.38-1  0.3588 

-0.01436  -0.1624 
-0.006381 

20.8  1 .879 

-0,0334  1  -1.472 


0 

1  527 
-0.01347 


0 

2.187 

-0.0236'; 


0.009622 

0.08506 

-1.166 


w.v;  1 

0.1249 

-1.712 


INNER  LOOP  (IAIN  SCHEDULES 


The  longitudinal  inner  loop  equalization  gain  schedules  are  given  by: 


F  =  -40 


Kf=ll  Ij  T  =  [0  .0247  Oj 

N=N(q)  M  =  [Mi(5)  M^(q)). 

N(q)  =  -.312  q  + 

MjCq)  = -.058  q +  50.5  MgCq)  = -.006  q  +  8.11  (A5.3) 

where  q  is  the  dynamic  pre.ssure.  The  lateral/directional  inner  loop  equalization 
schedules  are  given  by 

gain  =  C2  q^  +  Cl  q  +  CO  ,  {A5.4) 

where  the  coefficients,  Cj,  are  given  in  Table  A5.3. 

OUTER  LOOP  CONTROLLERS 

The  full  Oirder  and  reduced-order  outer  loop  robust  performance  controllers 
have  the  folio wiVij^  focTm: 

[  numll  numi2 
I  r  11021  num22 

. .  '"“derr  ' 

nuin  --  X(s  +  zi)(s  +  -is  +  z,i.iHs  +  Zu) 

den  ~  (s  •!  1- iXs -t- P2)  -  0s pn-iKs  +  Pi,) .  (A5.5) 

The  poles,  zeros,  and  high  iioquency  gains  of  eq.  (Ab.b)  arc  given  in  Tables  A5.4, 
/\5.5,  A5.6,  and  A5.7. 
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Table  A5.3  Inner  Loop  Lateral/Directional  Gain  Functions 


Gain 

Altitude 

C2 

CO 

p  ->  Pc 

all 

-3.7683e05 

4.8235e-02 

3.5702e+00 

P  -»Pc 

all 

2.2130e-06 

4.3849e-03 

-1.1015e+00 

r->Pc 

10  kft 

-1.7065e-07 

3.0285e-04 

-6.5368e-01 

r-^Pc 

20  kft 

4.7299e-07 

-1.1692e-04 

-5.0102e-0i 

r->Pc 

30  kft 

3.4239e-07 

-5.3058e-05 

-4.3078e-01 

P 

all 

4.390  le-06 

-1.6425e-02 

9.3445e+00 

r  ->fc 

10  kft 

6.2272C-07 

-1.5514e-04 

-3.8156e+00 

r  ->  fc 

20  kft 

2.3004C-07 

3.1432e-04 

-3.973{)e+00 

r  fc 

30  kft 

-1.3339e-07 

3.252  le-04 

-4.0197e+00 

U14 


Table  A5.4  Full  Order  Outer  Loop  Lontritudinal  Controller  Tcrn’is 


Terms 

High  q  Controller 

K 

3.2460e+01 

zi 

-l.OOOOe+04 

Z2 

-5.6966e+01 

Z3 

-3.9957e+01 

Z4 

-3.2275e+01 

Z5 

ZG 

-l.l415e+0l- 

J9.8971e+0() 

Z7 

-9.10260+00+ 

j7.2716o+0() 

Z8 

-9.1026o+0()- 

j7.2716e+0() 

ZQ 

gKfgfUif 

ZlO 

-5.5028e+()0- 

j5.3695o+00 

Zll 

-8.9949e-01 

Z12 

-1.9290e+0() 

PI 

-1.0090c  t-U2 

P2 

-r?.6976o+01 

.  -P3 

-6.4479e+01 

P^ 

-3.8398c+()l 

Pf) 

-3.2i;ne+01 

P6 

-2.8024(‘+UI 

P7 

-2.1437C+01 

P8 

-8.640He+00+ 

,ja.0166c+00 

P9 

■imgEngfQijH 

8.2418C+00 

-7.2782e-02 

P12 

■8.8269e-01 

P13 

-1.9109O+00 

Low  q  Ciontrollcr 


1.16620+01 


-l.OOOOe+04 
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Table  A5.5  Reduced-Order  Outer  Loop  Lon^tudinal  C'ontrolier  Terms 


Coefficients 

High  q  Controller 

Low  q  Controller 

K 

1.1995e+02 

1.2951e+02 

-2.4180e+03 

-6.0168e+02 

Z2 

-3.8407e+00+ 

.i6.2507e+00 

Z3 

-3.8407e+00- 

6.2507e+00.i 

PI 

■■QbEBBSHI 

-5.9135e+01- 

.j5.l963e+01 

P2 

-5.9135e+0U 

.i5.1963e4-01 

P3 

-7.2319e-02 

-3.6321e-02 

P4 

-2.0509e+01 

-6.6426e+00 
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Table  A5.6  Full-Order  Outer  Loop  Lateral/Directiona 

Controller  'iV'i  ei;- 

Terms 

numll 

numl2 

num21 

nun)22  | 

K(pi) 

1.5684e+01 

4.9687e-02 

-1.2148e+02 

^'39360  {)i  1  ■' 

•i  '  /  2c-f03 

zi  (P2) 

-5.0000e+04 

-5.0000e+04 

-6.0000e+04 

:o:.-'.;02 

22  (P3) 

-4.6906e+02 

-3.8500e+03 

-4.6906e+02 

-4.6yoric+U3  ■;  ■ 

0'V2e:02 

22  (P4) 

-8.1594e+01 

-4.6906e+02 

-8.0915e+01 

-1. 14220+02  !  3 

1' 

^53  'c+Ol  ' 

24  (P5) 

-8.0915e+01 

-8.0915e+01 

-8.1594e+01 

-8.09150  1-01  1  3 

'5940+01 

25  (P6) 

-4.7416e+01+ 

.il.9234e+01 

-3.1594e+01 

-4.733  le+01+ 
jl.9169e+01 

-8.1394i!+0  i  '  -4 
‘  !’ 

o.'7i';rbi-r 

'■0  .  'C-IO] 

26  (P?) 

-4.7416e+01- 

.11.92346+01 

-4.727  le+01+ 
.jl.8037o+01 

-4.7331e+0j.- 

jl.9169e+01 

-4.7'^21e  :-()!:•  , 
i2.O2;;S(.;0l  L 

I  n  ).()]- 

:  '(H  0  i 

27  (ps) 

-4.7271e+01+ 

.jl.8037e+01 

-4.7271e+01- 

.jl.8037e+01 

-4.7271e+01+ 

.il.8()37e+01 

-4.742  ' 

j2.()2:TtM ')'  ;■ 

'.].()  0^. 

;  01 

28  (pg) 

-4.7271e+01- 

.jl.8037e+''l 

-4.7849e+01+ 

.j2.0195e+01 

-4.7271e+01- 

jl.8037e+01 

-4.727  te+()1+  ^  ‘i-. 

.  ‘  v  M:+()i- 
■  :■  '■  '-Oi 

29  (pio) 

-1.8854e+01 

-4.7849e+01- 

.j2.0195e+01 

-1.5469e+01 

-4,727 .lo+(;  :■  3 

jl.803>7i.MO  1  i  i/ 

CO');  l-OO-i- 
0  1!)  Ye  {-00 

210  (Pll) 

-1.1416e+01 

-1.9301e+01 

-5.5529e+00+ 

j7.1005e+00 

-L4948y-!-0i  !  -fi. 

f  i7 

4-U'  )+r()b- 
C  O.'V-  M)() 

211  (P12) 

-4.2560e+00+ 

j5.7923e+00 

-5.0577e+00+ 

j6.5979e+00 

-5.5529e+00- 

j7.1005e+00 

-3l8b2o-A)0  :•  ■  -3. 
.)7,049.'k.' '  00  .0 

i.’O  ■rOO+ 

0  (;0 

Z12  (P13) 

-4.2560e+00- 

.i5.7923e+00 

-5.05776+00- 

.j6.5979e+00 

-4.7010e+00+ 

j3.5892e+00 

-3.l302rC)v3  .  •. 

j7.0493oiOO  : 

'  {00 
.■■■  .■/o.i,.o() 

213  (P14) 

-4.8355e+0C+ 

.i3.6809e+00 

-5.0731e+00+ 

.il.0852e+00 

-4,7010e+00- 

,i3.5892e+00 

■4.4362C-I  110;  '  3 

.i2,023oo:on  i 

O';  '■■^.•l•oo 

214  (Plb) 

-4.8355e+00- 

j3.6809e+00 

-5.073  le+00- 
.i.l.0852e+00 

-6.8363e+00 

■4,430‘>''-,<.()  ' 

j2.()233..:.'0  '  ,■ 

'  ''ibo: 

I'.'O 

215  (pie) 

-5.7669e+00 

-4.8788e+00 

-2.3675C+00+ 

j2.5835e+00 

-5.9:}!,)e'  71 

; .  ;-(J0-. 
'-.(.'O 

216  (pi?) 

-3.1217e+00 

-3.1217e+00 

-2.3675e+00- 

j2.5835e+00 

-3.8oiV/i  1  1;  1  ,  - 

:  l  OO'' 

217  (pis) 

-3.1247e+0U 

-3.1160C+00 

-3.1222e+00 

-3.12!)«eiOO 

'  1  ()( 

218  (P19) 

3.308ie-01 

-1.7925C+00 

-3.1217e+00 

-3,121  Ve  :  O'-  ■  •  ^ 

I 

;r;;-03 

zi:;  (o2o) 

-G.0353e-04 

-1.8030e-02 

-7. 658  le -03 

-6.63120  03  j 

.V03’e-()2  "■ 

220  (P21) 

-1.8030e-(2 

-1.8279e-02 

-1.80646-02 

-f).7550c' 02  5  -7 

.  07 100-03 

Z21  (P22) 

-1.7979e-02 

-4.4855e-03 

-5.7726e-02 

-1.7970e-(i2 

1 

.7  /  -.llo  02 

Z22  (P23) 

-5.7737e-02 

-5.7682e-02 

-5.7726e02 

-r).7726(>  02  : 

I 

'  -  .'  'Oo  02 

Z23  (P24) 

-5.7726e-02 

-5.7726O-02 

-1.8030e-02 

-l.H03o>- : 

.  o 
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Table  A5.7  Reduced-Order  Outer  Loop  Lateral/Directional  Ccntroller  'i'enas 


Terms 


K 

(-) 

21 

(pl) 

22 

(P2) 

Z3 

(P3) 

24 

(P4) 

25 

(P5) 

26  (P6) 


28  (P8) 


3.3817e-01 

-1.6278e-h00 

-4,4999e-03 

-5.4958e-02 

-5.2459e-02 

-5.2431e-02 

-5.3326e-02 

-7,64i57e-03 

-6.0643e-04 

-2.5573e-02 

-2.5330e-02 

-2.4672e-02 

-2,5310e-02 

-2.4413e-02 

-4.484  le-03 

-7.6325O-03 

-8,64[)3e-03 

-5,2293e-02 

6  CONCLUSIONS 


This  document  provides  theory  and  examples  for  a  design  methodology  for  manual 
flight  control  system  design.  The  inner/outer  loop  structure  gives  a  versatile 
framework  for  the  application  of  advanced  multivariable  control  theory  to  aircraft 
control  problems.  The  engineer  has  the  freedom  to  choose  any  method  which 
achieves  the  design  goals.  This  freedom  is  demonstrated  in  the  successftal  application 
of  three  different  inner  loop  equalization  approaches  in  the  VISTA  F-16  and 
supermaneuverable  vehicle  design  examples.  This  framework  is  also  versatile  enough 
to  include  a  broad  range  of  design  specific  issues  such  as  the  management  of 
redundant  control  effectors,  flying  qualities  requirements,  modeling  uncertainties,  and 
nonlinearities. 

While  the  design  examples  presented  in  this  document  illustrate  a  broad  range  of 
design  issues,  they  still  do  not  represent  flight  ready  control  laws.  Issues  not 
addressed  in  this  study  that  are  worthy  of  further  attention  include:  turbulence 
responses,  sensor  noise,  changing  inertial  properties  due  to  different  loadings  and  fuel 
bum,  redundancy,  and  the  details  of  pilot  interfaces.  Continuing  work  wall  focus  on  a 
dynamic  inversion  approach  to  inner  loop  design  which  includes  inertial  coupling  and 
high  order  terms.  High  angle  of  attack  and  post-stall  control  will  be  studied  for  the 
supermaneuverable  vehicle  with  thrust  vectoring  and  the  F-16  with  vortex  control. 
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